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CHAPTER  1 
PROGRAM  SUMMARY 

1.1  INTRODUCTION 

Current  pavement  design  approaches  are  mostly  empirical;  based 
primarily  on  field  experience  and  engineering  judgement.  These 
approaches  involve  correlating  the  required  thicknesses  of  pavement 
components  to  certain  traffic  characteristics,  soil  strength  prop- 
erties obtained  from  specified  included  soil  tests,  and  the  properties 
of  the  layer  materials.  These  methods  have  been  somewhat  successful 
in  the  past;  however,  their  applicability  to  the  initial  design 
or  rehabilitation  of  current  pavement  structures  is  severely  limited. 
This  deficiency  is  vividly  illustrated  by  the  more  rapid  deterioration 
of  new  and  rehabilitated  pavement  structures  throughout  the  country 
than  anticipated  by  the  designers  of  these  structures.  That  is,  the 
experience  gained  from  one  set  of  materials,  traffic  loading  conditions, 
environmental  conditions,  etc.  cannot,  in  general  be  used  to  predict 
pavement  performance  for  a  different  set  of  materials,  traffic  loading, 
environment,  etc.  due  to  the  complicated  nonlinear  relations  which 
exist  among  these  parameters  and  pavement  performance. 

The  above  shortcoming  takes  on  additional  importance  with  the 
increasing  scarcity  and  cost  of  petroleum.  Accordingly,  a  need  exists 


for  the  development  of  realistic  pavement  design  and  analysis  procedures 
which  have  a  firm  theoretical  foundation.  These  procedures  should  be 
aimed  at  optimizing  pavement  performance  from  a  highway  user  and 
maintenance  point  of  view.  Overall  costs  should  be  minimized  through 
trade-offs  between  initial  costs,  recurring  costs,  maintenance 
strategies  and  pavement  performance.  Accomplishment  of  these  objectives 
requires,  in  particular: 

1.  Accurate  descriptions  of  loads  and  environments 
throughout  the  anticipated  useful  life  of  the 
pavement  system. 

2.  Accurate  determination  of  material  response 
and  failure  properties  for  each  material  in 
the  pavement  system  under  laboratory  conditions 
representative  of  field  use  conditions. 

3.  Accurate  predictions  of  the  stresses,  strains 
and  deflections  throughout  the  pavement  structure 
due  to  the  prescribed  loads  and  environments. 

4.  Accurate  predictions  of  distress  and  damage  due 
to  the  prescribed  loads  and  environments. 

5.  Accurate  relations  between  pavement  system  distress 
and  pavement  performance. 

6.  Verification  and  evaluation  of  pavement  performance 
predictive  models  under  field  use  conditions. 

7.  Modification  and  re-evaluation  of  predictive  models 
as  required. 

In  addition,  design  of  overlays  requires  in  site  methods  for  determin- 
ing reduced  or  equivalent  material  response  properties  (i.e.,  moduli 
or  compliances)  for  the  layers  of  the  pavement  system,  and  for  the 
evaluation  of  the  load  bearing  capability  of  the  existing  pavement 
structure  at  the  time  the  pavement  is  overlaid. 

The  Federal  Highway  Administration  (FHWA)  has  addressed  various 
aspects  of  the  above  requirements  under  Research  Project  5C  entitled 


"New  Methodology  for  Flexible  Pavements".  Significant  and  substantial 
progress  has  been  made  in  this  direction  through  the  development  of 
an  overall  pavement  design  framework,  known  as  VESYS  II,  which  is 
based  on  mechanistic  concepts  [l,  2]. 

The  present  report  summarizes  University  of  Utah  work  related 
to  the  modification  and  extension  of  VESYS  II.  These  extensions/ 
modifications  have  led  to  an  improved  pavement  design  framework  which 
is  referred  to  as  VESYS  TIM  herein. 

The  primary  improvements  to  VESYS  II  consist  of: 

1.  Extension  of  Three-Layered  Viscoelastic  solution 
to  an  N-layered  viscoelastic  solution. 

2.  Development  of  approximate  closed  form  probabilistic 
solution  for  stresses,  strains  and  deflections  for 
and  N-layered  viscoelastic  system. 

3.  Development  of  an  improved  rut  depth  model  which 
utilizes  simplified  laboratory  characterization 
tests. 

Additional,  minor  modifications  have  also  been  made  to  VESYS  II 

which  are  described  in  subsequent  sections  of  this  report. 

Initial  efforts  on  the  conceptual  development  of  an  optimizing 
and  decision  making  program  which  interrelates  pavement  structural 
make-up,  system  costs,  maintenance  strategies  and  service  life  were 
carried  out  under  subcontract  by  Austin  Research  Engineers,  Inc., 
Austin,  Texas.  The  results  of  their  study,  as  submitted  in  their 
final  report  to  the  University  of  Utah,  are  presented  in  Chapter  7. 

Typical  input  data  and  examples  are  included  in  the  discussions 
of  the  various  system  components.  Analytical  derivations  and 
supporting  documentation  for  certain  models  are  included  as  appendix 
material  to  avoid  disruption  of  the  continuity  of  the  main  report. 


Additional  program  documentation  and  a  user's  manual  are  contained 
in  Volume  II  of  this  final  report. 

The  computer  code  has  been  submitted  under  a  separate  cover. 

1.2  PROJECT  OBJECTIVES  AND  ACCOMPLISHMENTS 

The  original  project  work  statement  consisted  of  four  specific 
tasks,  in  addition  to  a  literature  review: 

1.  Develop  a  VESYS  II  viscoelastic  five  layer  solution. 

2.  Develop  a  closed  form  probabilistic  solution  for  stresses, 
strains  and  deflections  due  to  statistical  variations  of 
creep  compliance  at  various  times  and  rates  of  loading. 

3.  Develop  modular  capability  of  the  design  framework 
for  substitution  of  improved  primary  response  or 
damage  models  as  they  are  developed.  Include  simple 
models  for: 

a.  thermal  fatigue 

b.  fatigue  fracture  mechanics 

c.  permanent  deformation 

d.  power  spectral  representations  of  pavement  roughness. 

4.  Develop  concepts  to  interrelate  structural  make-up,  systems 
costs,  maintenance  strategies  and  service  life  in  an  optimiz- 
ing and  decision-making  program. 

All  technical  objectives  have  been  accomplished  or  exceeded.  In 
place  of  a  viscoelastic  five-layer  solution,  a  general  N-layer  visco- 
elastic solution  has  been  developed  using  the  quasi-viscoelastic 
method  of  solution.  This  approach  was  adopted  for  a  variety  of 
reasons. 

The  three-layer  VESYS  II  is  constrained  to  consideration  of 
incompressible  materials  and  re-derivation  of  the  solution  for  five- 
layers  within  the  standard  convolution  integral  approach  for  compressible 
as  well  as  incompressible  materials  (Poisson's  ratio  =  1/2)  represented 
a  considerable  undertaking  beyond  the  scope  of  the  present  contract. 


Furthermore,  introduction  of  a  closed  form  probabilistic  solution  for 
stresses,  strains,  and  deflections  within  the  convolution  integrals 
for  the  transient  viscoelastic  solution  is  complicated.  The  quasi- 
viscoelastic  method  of  solution  avoids  these  difficulties  as  a  series 
of  elastic  solutions  are  combined  using  instantaneous  values  of 
relaxation  modulus  or  creep  compliance  to  obtain  the  viscoelastic 
solution. 

The  original  intent  was  to  introduce  an  approximate  closed  form 
solution  into  the  CHEVRON  5-layer  program  [3].  However,  the  solution 
method  in  CHEVRON-5  is  specifically  applicable  to  a  layered  system 
with  no  more  than  five-layers,  whereas  the  approximate  closed  formed 
probabilistic  solution  was  obtained  by  a  Taylor's  series  expansion 
of  the  general  analytic  solution  for  N-layers  about  the  mean  values 
of  the  material  properties.  It  appeared  that  approximately  the  same 
effort  was  required  to  either  program  the  solution  to  the  probabilistic 
equations  using  the  method  employed  in  CHEVRON-5  or  to  completely 
re-program  the  solution  to  the  N-layered  problem.  Since  CHEVRON-5 
also  has  other  limitations  and  in  the  interest  of  greater  program 
efficiency  and  flexibility  it  was  decided  to  develop  a  new  computer 
code  for  the  deterministic  and  probabilistic  solution  for  a  N-layered 
system. 

An  approximate  closed  form  probabilistic  solution,  as  noted 
above,  was  obtained  using  Cornell's  first  order,  second  moment 
theory  [4].  A  Taylor's  series  expansion  is  taken  about  the  mean 
values  of  layer  creep  complicances.  The  deterministic  solution 
is  used  to  represent  the  mean  solution  and  the  variance  of  stresses, 
strains  and  deflections  is  expressed  in  terms  of  the  first  partial 
derivatives. 
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The  various  subsystems  are  assembled  in  a  modular,  overlay 
fashion;  however,  rather  than  relay  on  machine  dependent  system 
commands  to  control  execution  of  the  various  program  elements,  a 
simple  command  program,  written  in  Fortran,  controls  execution  of 
the  overlays.  This  procedure,  while  not  as  efficient  for  a  specific 
computer,  increases  program  portability. 

The  executive  program  controls  execution  of  all  program  overlays. 
The  three  main  elements  of  the  framework  consist  of  the  Primary 
Response  Model,  the  Distress  Model  and  the  Serviceability  Model, 
which  are  not  connected  by  any  common  statements.  Internal  logic 
in  the  executive  program  determines  which  subsystem  is  to  be  executed 
and  calls  a  Read  program  which  reads  the  necessary  input  data  and 
creates  a  data  file  for  access  by  the  subsystems.  In  this  way,  the 
various  subsystems  can  be  executed  separately  and  readily  replaced 
as  improved  mechanistic  models  are  developed. 

1.3  DESIGN  AND  ANALYSIS  METHODOLOGY 

The  recommended  methodology  for  design  and  analysis  of  pavement 
structures  consists  of  the  four  phases  shown  schematically  in  Figure  1-1. 

1.  Problem  Definition 

2.  Primary  Response  Analysis 

3.  Distress  Analysis 

4.  Verification 

Additionally,  a  fifth,  surveillance  phase,  consists  of  periodically 
updating  analyses  and  pavement  structural  integrity  assessment  during 
the  service  life  of  the  pavement  system.  The  approach  depicted  by 
the  diagram  is  typical,  but  may  be  modified  for  special  types  of 
analyses. 
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FIGURE  1.1   PAVEMENT  DESIGN  METHODOLOGY 


The  scope  of  analysis  for  a  particular  situation  must  be  determined 
in  early  phases  of  the  analysis  so  that  data  requirements  can  be  clearly 
defined.  Although  the  nature  and  precision  of  the  data  required  will 
depend  upon  the  type  of  analysis,  minimum  requirements  generally 
consist  of: 

1.  Description  of  the  Pavement  Structure  and  System 
Components. 

2.  Layer  Material  Response  and  Failure  Properties. 

3.  Definition  of  Loads  and  Environments. 

4.  Required  Service  Life. 

5.  Maintenance  Requirements. 


CHAPTER  2 


OVERVIEW  OF  PAVEMENT  DESIGN 
&  ANALYSIS  METHODS 


2.1   INTRODUCTION 

Much  effort  during  the  past  twenty  years  has  gone  into  the 
development  of  pavement  design  procedures.  Some  of  these  methods  are 
widely  used  and  firmly  established  with  some  design  agencies.  These 
design  methods  can  be  classified  into  three  groups:  Elasticity 
methods,  Semi-Empirical  and  Statistical  methods,  and  Empirical  methods. 

Elacticity  methods  consider  the  behavior  of  pavement  under  traffic 
loading  assuming  the  pavement  system  remains  proportional  to  the  applied 
loads.  The  principal  design  criterion  consists  of  limiting  stresses 
or  strain  as  calculated  by  an  elastic-layered  theory.  The  methods 
belonging  to  this  group  include  those  used  by  the  Kansas  Highway 
Department  [5],  the  U.  S.  Corps  of  Engineers  (CBR  Method)  [6],  Shell 
Company  [7]  and  the  Texas  Highway  Department  [8]. 

Semi-Empirical  and  Statistical  methods  are  based  mainly  on  the 
correlation  of  field  performance  with  certain  "indirect"  strength 
properties  of  pavement  materials.  Methods  belonging  to  this  group 
include  those  of  the  state  of  California  [9,  10],  the  McLeod  Method  [11], 
the  AASHO  Intermin  Guidelines  for  Flexible  Pavement  design  [12],  and 
the  Asphalt  Institute  Method  [13]. 
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FIGURE  2.1  Partial  Listing  of  Categories  of  Pavement  Distress 
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Empirical  methods  relate  the  pavement  thickness  to  some  particular 
soil  conditions.  Only  the  Michigan  Highway  Department  Method  [14]  and 
the  Canadian  Good  Road  Association  Method  belong  to  this  group. 

Though  some  of  these  methods  have  been  used  widely  and  in  some 
sense  successfully,  none  of  these  methods  have  the  ability  to  fully 
take  into  account  all  of  the  pavement  distress  associated  with  traffic 
and  non-traffic  loading,  as  illustrated  in  Figure  2-1.  In  addition, 
the  rapid  changes  in  traffic  volumes,  the  steadily  increasing  wheel 
loads,  the  use  of  new  and  sometimes  poorer  quality  materials,  more 
strength  maintenance  requirements  and  longer  required  service  life 
and  the  variability  of  environmental  conditions  make  these  methods 
unsatisfactory.  Accordingly,  more  rational  pavement  design  procedures 
capable  of  considering  each  distress  mechanism  have  been  developed 
or  are  under  development  by  several  agencies.  One  of  the  first 
rational  pavement  design  frameworks  was  developed  by  the  Massachusetts 
Institute  of  Technology  under  contract  from  the  Federal  Highway 
Administration  and  is  known  as  the  VESYS  II  viscoelastic  pavement 
design  framework  [1]. 

2.2  REVIEW  OF  VESYS  II  PAVEMENT  DESIGN  FRAMEWORK 

The  VESYS  II  pavement  design  system  framework  developed  by  MIT 
and  shown  schematically  in  Figure  2-2,  consists  of  five  components 
which,  when  integrated  together,  result  in  a  predictive  pavement 
design  system.  Each  component  consists  of  several  elements  of  sub- 
models which  form  the  rational  basis  for  the  design  procedure.  A 
brief  review  of  the  VESYS-II  is  given  in  this  chapter.  More  detailed 
information  is  provided  in  references  1,  2  and  15  to  18. 
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FIGURE  2-2  VESYS  II  Pavement  Design  System 
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In  the  VESYS  II  pavement  design  framework  the  pavement  system  is 
represented  by  a  three-layer  model  (Figure  2-3);  the  upper  two  layers 
have  finite  thicknesses  and  the  third  layer  is  infinite.  Each  layer 
is  infinite  in  the  horizontal  directions  and  may  have  elastic  or 
viscoelastic  behavior. 

The  pavement  system  is  subjected  to  axisymmetric  loading  with 
varying  intensities  and  durations  represented  by  statistical  distributions, 
Montly  mean  temperatures  are  used  to  represent  the  prevailing  regional 
climatic  conditions.  The  effects  of  moisture  on  the  subgrade  response 
is  not  considered. 

Much  of  the  input  and  output  is  described  in  terms  of  statistical 
distributions  instead  of  single-valued,  deterministic  estimates.  These 
probabilistic  formulations  are  of  two  types:  the  method  of  statistical 
trials  (monte  carlo  simulation)  and  closed  form  probabilistic  solutions. 
For  the  most  part,  simulation  techniques  are  used  in  the  structural 
subsystem,  although  some  approximate  closed  form  solutions  are  employed. 

A  block  diagram  of  the  structural  model  is  presented  in  Figure  2-4. 
The  structural  subsystem  consists  of  a  primary  response  model,  a  damage 
model  and  a  serviceability  model.  In  the  primary  response  model,  the 
response  of  the  layered  system  to  a  stationary,  moving  or  repeated 
load  applied  at  the  surface  is  obtained.  Using  the  viscoelastic  three- 
layered  system  program,  materials  properties  variations  are  accounted 
for  using  Monte  Carlo  techniques  and  approximate  closed  form  probabilistic 
solutions.  The  statistical  estimates  of  stresses,  strains  and  deflections 
so  determined  along  with  the  loading  characteristics  and  temperature 
history  are  used  as  inputs  to  the  Damage  Model  to  obtain  the  extent 
of  rut  depth,  roughness  or  slope  variance,  and  the  extent  of  cracking. 
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FIGURE  2-3  Idealization  of  Pavement  Structure  in  VESYS  II 
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The  distress  indicators  obtained  from  the  distress  model  are  combined 
in  the  form  of  a  equation  suggested  by  AASHO  to  provide  a  subjective 
measure,  referred  to  as  the  Present  Serviceability  Index  (PSI),  of 
the  serviceability  level  of  the  system  at  any  time  period. 

The  different  elements  of  the  VESYS  II  design  framework  are 
discussed  briefly  in  the  following  subsections. 

2.2.1  INPUT  VARIABLES 

The  inputs  to  the  structural  model  are  divided  into  four 
categories:  system  geometry,  materials  properties,  loads  character- 
istics and  temperature  history. 

The  inputs  for  the  SYSTEM  GEOMERTRY  consist  of  layer  thicknesses 
for  the  first  two  layers.  The  thicknesses  are  input  as  deterministic 
quantities. 

The  MATERIALS  PROPERTIES  input  for  each  layer  consist  of  creep 
or  elastic  compliances,  time-temperature  shift  factors,  fatigue  data 
in  terms  of  strain  level  versus  number  of  cycles  to  failure,  and  the 
spatial  variation  of  materials  properties.  Each  layer  is  assumed  to 
be  incompressible  (v  =  1/2). 

The  Creep  compliance  data  are  fitted  to  an  exponential  (Dirichlet 

or  Prony)  series  of  the  form 

N 
J(t)  =  n  L     J.exp(-t/_  )  (2-1) 

i=l  1  Ti 

using  a  least  squares  curve  fitting  subroutine.  In  Eq.  (1),  J(t) 

represents  the  creep  compliance  at  time  t,  and 

n  =  random  variable  with  a  mean  of  unity  and  a  coefficient 

of  variation  determined  by  the  scatter  of  the  creep 

test  data. 
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J.  =  coefficients  for  the  mean  creep  compliance 
t.  =  retardation  time  constants 
N  =  number  of  terms  in  the  series 
When  considering  the  effects  of  environmental  factors  such  as 

temperature,  Eq.  (1)  is  written: 

N 
J(t,J)  =  ay  +  BTri  I     J1exp(-t/T|<aT)  (2-2) 

where  ou-  represents  a  vertical  shift  on  an  arithmetic  scale,  6T 

represents  a  vertical  shift  on  a  logarithmic  scale  and  aT  represents 

a  horizontal  shift  on  a  logarithmic  scale  of  time.  Assuming 

thermorheologically  simple  material  behavior  ou-  =  0,  3T  =  p  T  -  1 

o  o 
where  p  denotes  the  density  and  the  subscript  zeros  indicate  reference 

values.  The  term  aT  is  referred  to  as  the  time-temperature  shift 

factor  and  is  a  function  of  temperature  (Rankine  or  Kelvin)  only. 

Fatigue  test  data  are  represented  by  the  equation 

N  =  k1Ae"k2  (2-3) 

where  N  represents  a  number  of  cycles  to  failure  at  a  strain  amplitude 
of  Ae.  k,  and  k?  are  material  property  coefficients,  k,  and  k?  are 
statistically  correlated  and  their  temperature  dependence  determined 
experimentally.  The  coefficients  k,  and  k^   are  given  by  their  means, 
variances  and  coefficients  of  correlation. 

The  spatial  variation  of  materials  properties  is  given  by  a 
spatial  auto-correlation  function  which  measures  the  statistical 
correlation  of  two  points  a  given  distance  apart.  This  coefficient 
denoted  by  c,  is  represented  by  the  equation 

C  =  Aexp(-Ax2/B2)  (2-4) 

where  Ax  is  the  distance  between  the  two  points,  and  A  and  B  are  constants, 
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The  traffic  LOAD  HISTORY  consists  of  the  duration,  intensity 
and  frequency.  A  circular  contact  area  is  assumed  for  the  wheel 
loads  with  the  contact  area  assumed  to  be  constant.  Hence,  the 
variability  of  the  wheel  load  is  expressed  in  terms  of  mean  and 
variance  of  the  tire  pressure.  The  mean  and  the  coefficient  of 
variation  of  load  duration  are  used  to  represent  the  statistical 
variation  of  wheel  loads  duration.  Frequency  of  loading  is  expressed 
in  terms  of  a  mean  number  of  load  applications  per  month.  In 
addition,  a  traffic  load  channelization  factor  is  also  input. 

The  ENVIRONMENTAL  VARIABLES  are  assumed  to  be  temperature  and 
moisture;  although  moisture  effects  are  yet  explicitly  included. 
Temperature  is  given  by  its  mean  and  variance  for  each  basic  time 
period,  which  is  usually  taken  to  be  a  month. 

2.2.2  PRIMARY  RESPONSE  MODEL 

The  method  of  solving  the  three-layer  (linear)  viscoelastic 
boundary  value  problem  consists  of  three  steps: 

1)  the  elastic  solution  due  to  a  stationary  load  is  obtained; 

2)  the  elastic-viscoelastic  corresponding  principle  is  applied 
to  the  elastic  solution  to  obtain  the  corresponding 
viscoelastic  solution; 

3)  the  stress-strain  and  deflection  due  to  a  stationary  load 
together  with  the  time  functions  representing  the  repeated 
and  moving  loads  are  used  to  generate  the  stress,  strain 
and  deflection  due  to  repeated  and  moving  loads  using 
linear  superposition. 
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In  the  formulation  of  the  three-layer  elastic  solution,  the  loading 

and  geometry  of  the  system  are  assumed  to  be  axial ly  symmetric.  Therefore: 

the  equilibrium  equations  and  compatibility  equations  expressed  in 

polar  coordinates,  and  the  stress-strain  and  strain-displacement 

relations  for  an  incompressible  body  form  a  complete  set  of  equations 

for  solution  of  the  stress,  strain  and  deflection  in  each  layer.  By 

using  a  proper  set  of  exterior  boundary  conditions  (on  the  free 

surface),  interior  boundary  conditions  (continuity  on  the  interfaces) 

and  the  boundary  condition  at  infinity,  together  with  the  governing 

equations  in  each  layer,  the  stress,  strain  and  deflection  anywhere 

in  the  three-layer  system  can  be  uniquely  determined.  For  example, 

the  deformation  at  the  free  surface,  w,  is  expressed  in  the  form: 

18 

w  =  2(1.wi)2¥^iw -|!i!i  (2.5) 

A0* 

In  this  expression  the  <b. 's  and  8.;'s  depend  on  the  geometry  (thickness 

of  each  layer),  and  the  a.'s  and  g.'s  depend  on  the  products  of  the 

elastic  properties  of  each  layer.  For  example, 

2  2       2  3 

a,  =  E-. E? ,  cu  =  E^TEpE-,  and  ou=  E-,EpE^,  etc. 

In  the  viscoelastic  solution  of  a  three-layer  system  under  a 
stationary  load,  it  is  only  necessary  to  replace  the  modulus  for 
each  layer  in  the  previous  solution  by  the  corresponding  viscoelastic 
properties  (relaxation  modulus  or  creep  compliance)  using  the  elastic- 
viscoelastic  correspondence  principle.  For  example,  for  the  deflection 
at  the  free  surface,  Eq.  (5)  becomes 
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9  t  93, (t) 

I     8.  S    w(t-x)  — 5 dx  =  R  E  (J). a.  (2-6) 

i=l     n     o  dT  V1 

2 
where  3-  and  a.  are  integral  operators.     For  example,  when  ou  =  E,E?Er 


then 


«4(t)  ■  t\  E^t-x,)  g-    £   Efy   -  x2)  £-  [/J2  E2(x2  -  x3) 

3E3<T3>  .  , 
-SrT- dT3] 


and 


t       9B4(x) 


(2-7) 


=  /0tw(t-x)|T[VE1(t-x1)^T/(;iE1(Tl  -x2) 

9     T2  T3  8E3^T4^ 

^-  /0  El;%  -  t3)  X  /o  E2(t3  -  t4)  -^—  dx4] 

Equation  (6)  is  an  integral  equation  for  the  deflection  w.  E, (t),  Ep(t) 
and  E^(t)  are  the  relaxation  moduli  for  the  three  layers.  Numerical 
solution  of  this  equation  is  straightforward  although  it  requires 
very   lengthy  calculations,  particularly  if  reasonable  accuracy  is 
desired.  Once  the  primary  response  parameters  (stress,  strain  and 
deflection)  for  a  stationary  load  are  obtained,  then  the  corresponding 
parameters  for  repeated  and  moving  loads  are  obtained. 

The  approach  mentioned  above  for  the  solution  of  three-layer 
viscoelastic  system  is  sound.  However,  it  is  not  without  faults. 
In  particular,  the  computation  time  required  in  obtaining  the  numerical 
solution  of  Eq.  7  is  substantial.  It  is  well  known  that  the  accuracy 
of  the  solution  depends  upon  the  fineness  of  the  time  interval  chosen. 
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Greater  accuracy  can  be  achieved  using  smaller  time  intervals  in  the 
numerical  solution.  However,  the  computation  time  is  substantially 
increased  as  the  number  of  time  intervals  increases. 

The  stationary  loading  program  yields  the  primary  responses  of 
the  system  which  include  the  deflection  at  the  surface  and  the  radial 
strain  at  the  first  interface  beneath  the  center  of  the  loaded  area. 

The  curve  fitting  subroutine  is  used  to  obtain  an  exponential 
series  representation  of  the  primary  responses.  This  series  along 
with  the  traffic  and  temperature  histories  and  the  fatigue  properties 
of  the  materials,  are  then  used  as  input  to  the  repeated  loading 
subroutine. 

Assuming  the  loads  to  be  functions  of  time,  the  response  after 
repeated  loading  is  determined  using  linear  superposition.  This 
deterministric  solution  is  then  converted  into  an  approximate 
probabilistic  closed  form  solution.  This  probabilistic  solution  for 
the  primary  responses  (i.e.,  stresses,  strains,  and  deflections) 
is  input  to  the  damage  model  to  predict  pavement  rutting,  cracking 
and  slope  variance. 

The  damage  model  yields  the  history  of  the  statistical  characteristics 
of  the  damage  indicators  which  are  input  to  the  serviceability  model 
as  described  subsequently. 

2.2.3  OUTPUT 

The  outputs  of  the  structural  model  include  the  statistical 
estimates  of  AASHO's  damage  indicators  (i.e.,  rut  depth,  roughness 
and  cracking)  determined  at  different  times. 
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Rut  depth  is  considered  to  be  mainly  due  to  channelization  of 
traffic  which  causes  a  differential  deformation  in  the  transverse 
direction,  and  is  given  as  a  differential  settlement  between  the 
most  travelled  channel  and  its  surroundings. 

Slope  variance,  or  roughness,  is  an  indication  of  the  state  of 
the  pavement  in  the  longitudinal  direction,  and  is  caused  mainly 
by  variation  in  the  properties  of  the  materials.  Layers  which  are 
assumed  to  be  elastic  contribute  only  to  the  resilient  deformations 
of  the  pavement  system. 

A  linear  Miner's  law  is  used  to  account  for  crack  accumulation 
in  the  top  layer  assuming  that  a  fatigue  mechanism  is  responsible 
for  progression  of  cracking  of  pavements.  The  strains  at  the  first 
interface  are  used  to  determine  the  fatigue  life  employing  Eq.  (3). 

2.2.4  SERVICEABILITY  MODEL 

The  distress  indicators  obtained  from  the  damage  model  (i.e., 
rut  depth,  slope  variance  and  cracking)  are  input  to  the  serviceability 
model  shown  in  Figure  2-5  to  define  the  life  expectancy  of  the  pavement 
system.  Figure  2-5  also  includes  the  maintenance  model  and  highway 
cost  model . 

The  distress  indicators  obtained  from  the  damage  model  are 
combined  in  a  linear  regression  form  suggested  by  AASHO  to  obtain 
the  Present  Serviceability  Index  (PSI). 

PSI  =  aQ  +  a]X1  +  a£X2  +  a3X3  (2-8) 

where  X,,  X?  and  X-  are  statistical  measures  of  the  damage  caused 
by  rutting,  slope  variance  and  cracking.  The  probabilities  of  having 
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any  value  of  PSI  at  any  time  are  referred  to  as  the  state  probabilities, 
and  the  probability,  at  any  time,  of  being  above  some  unacceptable 
value  of  PSI  is  defined  as  the  reliability  of  the  system  at  that  time. 
The  unacceptable  limit  of  serviceability  is  subjectively  defined  by 
the  user  and  defines  some  threshold  for  failure,  from  which  the  life 
expectancy  of  the  system  is  determined. 

The  outputs  of  the  serviceability  model  consist  of  the  following: 

1)  Means  and  variances  of  AASHO's  PSI  at  different  times 
using  the  following  expression: 

PSI  =  aQ  +  a]  log[l  +  (SV)]  +  a2/[GAT  +  a3(RD)2  (2-9) 

where  (SV)  represents  the  slope  variance,  (CA)  represents 
the  cracked  area  and  (RD)  represents  the  rut  depth.  If 
the  damage  indicators  are  statistically  independent,  the 
following  equations  for  the  expected  value  of  the  PSI 
is  used: 

E[PSI]  =  ^   +  a2{log[l+(SV)]  -  ^£l+(SV)2  a$v2]} 

3/2     2 

+  *3{/m  °Am}  +  34{0RD  +  ^] 

(2-10) 

2 
The  upper  bars  represent  an  expected  value  and  the  a 

represent  a  variance.  The  value  of  the  variance  of  the 

PSI  is  similarly  derived. 

2)  Assuming  that  the  serviceability  index  has  a  normal  distri- 
bution and  that  an  unacceptable  level  of  serviceability  is 
defined,  the  reliability  of  the  system  at  a  given  time  is 
computed  as  the  probability  of  having  the  system  above  the 
failure  level . 
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3)  The  expected  life  of  the  system  is  given  as  the  mean  time 
for  the  system  to  reach  the  failure  level  in  the  absence 
of  maintenance  activities.  Probabilistic  distributions 
of  the  life  are  also  given. 

4)  The  probability  of  the  system  being  between  specified 
levels  of  the  serviceability  index  (i.e.,  the  marginal 
state  probability). 

5)  The  state  transition  probability,  i.e.,  the  probability 
of  going  from  one  level  to  another  at  a  specified  time. 

2.2.5  MAINTENANCE  SUBMODEL 

The  maintenance  submodel  shown  in  Figure  2.5  generates  various 
maintenance  strategies  over  the  operational  life  of  the  pavement 
system  and  evaluates  the  consequences  of  these  strategies  in  terms 
of  serviceability,  reliability  life  and  economic  attributes.  A 
decision-making  algorithm  selects  those  strategies  which  provide 
optimal  costs  for  a  given  design  configuration. 

2.2.6  COST  SUBMODEL 

The  cost  model,  also  shown  in  Figure  2-5,  addresses  itself  to 
the  determination  of  the  total  costs  of  construction,  operation  and 
maintenance  of  highway  systems,  and  incorporates  three  components: 
construction  costs,  roadway  maintenance  costs  and  vehicle  operating 
costs.  Each  component  in  this  model  provides  estimates  of  resource 
consumption,  and  yields  the  cost  estimates  of  these  resources  using 
separately  defined  unit  prices. 
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Input  variables  define  the  project  to  be  analyzed  and  the 
structural  design-maintenance  strategy  to  be  evaluated.  Grade, 
alignment,  width  and  depth  of  surfacing,  maintenance  and  reconstruction 
policy  are  specified  along  with  prevailing  unit  labor,  materials  and 
equipment  costs,  the  submodel  estimates  the  construction,  maintenance 
and  user  costs  for  each  analysis  period  which  are  discounted  to  find 
the  present  value  of  the  structural  design  -  maintenance  strategy 
that  is  being  evaluated.  The  result  is  the  ability  of  the  highway 
designer  to  compare  a  number  of  strategies  on  the  basis  of  their 
total  predicted  costs. 

These  latter  two  submodels  have  not  been  fully  implemented 
with  the  other  models  of  the  structural  subsystem. 
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CHAPTER  3 

SYSTEM  INPUTS  -  LOADS  AND  ENVIRONMENTS 

3.1   INTRODUCTION 

The  system  inputs  for  VESYS-IIM  consist  of  the  following  six 
categories: 

(1)  Geometry  of  the  pavement  system 

(2)  Traffic  Loading 

(3)  Temperature 

(4)  Moisture 

(5)  Material  Response  Properties 

(6)  Material  Damage  Properties 

All  the  inputs,  except  the  geometry  of  the  pavement  system, 
are  in  probablistic  format  such  that  the  variability  of  the  various 
parameters  can  be  properly  taken  into  consideration  in  evaluating 
the  response  and  performance  of  the  pavement  system.  Thus,  the 
variability  of  the  pavement  performance  as  well  as  the  mean  performance 
of  the  pavement  system  can  be  estimated  in  a  more  realistic  manner. 

From  a  practical  point  of  view,  in  the  course  of  design,  the 
inputs  for  the  system  are  kept  as  simple  and  straight  forward  as 
possible.  Ultimately,  it  is  hoped  that  all  intermediate  data  reduction 
can  be  handled  by  a  computer.  Some  of  these  intermediate  calculations 
reguire  additional  knowledge  and  experience  from  the  one  doing  the 
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data  reduction. 

In  the  chapter  the  first  four  categories  will  be  included.  The 
last  two  categories  with  regard  to  the  material  characterization  are 
discussed  in  the  following  chapters. 

3.2  GEOMETRY  OF  THE  PAVEMENT  SYSTEM 

The  geometry  requirement  for  the  pavement  system  is  the  heights 
of  the  first  n-1  layers  with  the  thickness  of  the  nth  layer  being 
infinite.  Therefore,  for  a  5  layer  system  it  is  necessary  to  specify 
the  heights  of  the  first  four  layers. 

3.3  TRAFFIC  LOADS 

The  primary  responses  (stresses,  strains  and  deflections) 
depend  upon  the  total  weight,  the  contact  area,  the  contact  pressure 
and  wheel  configurations.  The  manifestation  of  pavement  distress 
in  the  form  of  fatigue  cracking,  rutting  and  roughness,  on  the  other 
hand,  depends  on  the  stresses,  strains  and  deflections  as  well  as 
total  number  of  load  applications.  Hence  a  more  realistic  way  of 
handling  the  traffic  loading  is  needed.  Instead  of  just  reading  in 
a  certain  fixed  intensity  (such  as  equivalent  18  )  it  is  proposed  to 
input  (1)  the  actual  wheel  load  distribution  such  as  the  distribution 
from  the  standard  W-4  form  which  provides  the  number  of  axles  per 
thousand  trucks  (or  vehicles)  in  each  axle  load  group  (10  groups  for 
single  axles,  and  converting  the  tandem  axles  into  equivalent  single 
axles);  (2)  Average  Daily  Traffic,  ADT  (including  mean  and  variance); 
(3)  monthly  variation  of  traffic;  (4)  average  growth  of  traffic  per 
year;  (5)  speed  of  trucks  (mean  and  variation);  (6)  average  radius 
of  the  contact  area;  (7)  channelization  distribution  factor;  (8)  percent 
of  truck  in  traffic  stream;  (9)  percent  trucks  in  the  design  lane. 
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Input  (1)  provides  sufficient  information  with  regard  to  the  wheel 
load  distribution.  Input  (2)  and  (3)  allows  generation  of  the  monthly 
traffic.  Input  (4)  is  used  to  estimate  future  traffic  growth. 
Inputs  (5)  and  (6)  are  used  to  calculate  the  duration  and  (7),  (8) 
and  (9)  are  used  for  calculation  of  the  percent  of  traffic  under 
the  wheel  paths  in  the  design  lane. 

In  the  present  VESYS  II  primary  response  model,  the  contact 
area  (radius  =  R)  is  considered  to  be  a  fixed  value.  Hence,  the 
load  intensity  (or  pressure  p)  is  directly  proportional  to  the  wheel 
load  P. 

P=   P 


ttR' 


Since  VESYS  II  is  a  linear  system,  that  is  the  primary  responses, 
strains  and  deflections  at  any  position  in  the  pavement  system  are 
proportional  to  the  load  intensity  (pressure),  the  primary  responses 
are  thus  linearly  proportional  to  the  wheel  load.*  Therefore,  in 
the  calculations  of  the  primary  responses  [o(l),  e(1)  and  5(1)]  only 
the  response  from  a  unity  pressure  need  to  be  calculated.  The 
primary  response  of  a(l),  c(l)  and  5(1)  correspond  to  a  wheel  load 
of  1  x  (ttR2).  The  corresponding  primary  responses  for  different 
axle  load  groups  can  be  obtained  straight  forwardly  by  multiplying 
the  factor  which  equals  the  wheel  load  in  each  axle  group  divided 
by  (2ttR2).  Thus  a  table  which  consists  of  frequency  versus  stresses, 
strains  and  deflections  per  thousand  vehicles  can  be  obtained  from 


*Notice  that  this  argument  of  proportionality  between  the  primary 
responses  to  the  wheel  load  is  not  valid  if  the  contact  area  is  varied. 
For  example,  for  a  constant  pressure,  the  contact  radius  is  proportional 
to  VP~.  Non linearity  is  introduced. 
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input  (1).  Input  (2),  (3)  and  (4)  is  then  used  to  calculate  the 
total  number  of  stress,  strain  and  deflection  intensities  in  each 
category  in  each  month,  each  year  and  total  accumulated  values  for 
any  period.  This  information  is  then  input  to  the  damage  models  to 
calculate  the  total  number  of  stress,  strain  and  deflection  intensities 
in  each  category  in  each  month,  each  year  and  total  accumulated 
values  for  any  period.  This  inforamtion  is  then  input  to  the  damage 
models  to  calculate  the  distresses.  A  flow  diaqram  is  shown  in 
Figure  3-1 . 

At  the  present  time  the  effects  of  wheel  load  configurations, 
and  shear  force  on  the  stresses  induced  in  the  pavement  system  are 
not  included  in  the  primary  response  model  of  the  VESYS  IIM  system. 
This  is  mainly  due  to  the  assumption  of  axial  symmetry  for  loading 
as  well  as  geometry  on  the  solution  of  layered  system  which  excluded 
the  multiple  wheel  loads  and  shear  forces  automatically.  Although 
it  should  not  be  too  difficult  a  incorporate  non-symmetric  loading 
into  the  layered  system.  The  computer  program  developed  recently 
by  Harrison,  Wardle,  and  Gerrard  [19]  in  Australia  can  handle  the 
horizontal  shear  force. 
3.4  TEMPERATURE 

The  need  for  temperature  is  two  fold:  one  is  the  calculation  of 
time- temperature  superposition  used  to  calculate  the  effect  of  temperature 
on  the  creep  compliances  of  the  materials  and  the  other  is  to  calculate 
thermal  fracture.  The  temperature  is  given  by  the  mean  and  variance 
of  air  temperature  for  each  basic  period.  The  basic  period  is  usually 
taken  to  be  a  month.  This  data  may  be  obtained  from  the  weather 
bureau  in  the  region. 
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The  temperature  in  the  pavement  can  be  obtained  from  a  heat 
transfer  analysis. 
3.5  MOISTURE 

Subgrade  moisture  is  read  in  to  calculate  the  moisture- 
shift  factor  to  be  used  to  relate  the  effect  of  moisture  on  the 
creep  compliance  of  subgrade  materials.  This  information  is  also 
needed  in  order  to  calculate  the  freezing  and  thawing  effect  of  the 
subgrade  in  each  basic  period. 

In  the  VESYS  I IM  it  is  assumed  that  surface  and  base  materials 
are  relatively  insensitive  to  moisture  change.  Hence  moisture  effect 
is  considered  only  in  the  subgrade  materials. 
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CHAPTER  4 

SYSTEM  INPUTS  -  MATERIALS  PROPERTIES 

4.1   INTRODUCTION 

The  properties  of  the  materials  in  each  layer  in  pavement  system 
are  dependent  upon  the  constituents  of  the  materials  and  the  manner 
in  which  they  are  prepared  and  constructed.  Since  these  factors  are 
statistical  quantities,  the  material  properties  in  each  layer  are 
expected  to  be  statistical  quantities.  The  material  properties  which 
are  pertinent  to  pavement  performance  are  divided  into  the  following 
two  categories: 

(1)  Material  Response  Properties:  Creep  compliances  (for 
viscoelastic  materials)or  elastic  modulus  (for  elastic  response) 

and  Poisson's  ratio.  These  properties  are  needed  in  order  to  calculate 
the  stresses,  strains  and  deflection  responses  in  a  pavement  system 
under  the  application  of  external  loadings. 

(2)  Material  Damage  Properties:  Fatigue  properties  and  permanent 
deformation.  These  properties  together  with  the  stresses,  strains  and 
deflection  responses  in  the  pavement  system  are  used  to  estimate  the 
pavement  damage  of  cracking,  rutting  and  roughness  under  various 
stages  of  its  life. 
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4.2  CREEP  COMPLIANCES  AND  ELASTIC  MODULUS 

In  the  VESYS  IIM  framework,  the  pavement  system  is  considered 
to  be  an  n  layer  linear  viscoelastic  system.  For  example,  in  a  five- 
layered  system,  the  first  and  second  layers  represent  the  bituminous 
surface  course  and  base  course  respectively,  the  third  layer  represents 
the  compacted  s  bgrade,  the  fourth  and  fifth  layers  can  be  considered 
to  be  natural  ground.  The  materials  in  each  layer  are  considered 
to  be  linear  viscoelastic  or  linear  elastic  if  the  time-dependency 
of  the  material  response  under  a  constant  stress  can  be  neglected. 

4.3  LINEAR  VISCOELASTIC  CONSTITUTIVE  EQUATION 

For  a  linear  viscoelastic  material,  the  creep  strains  e.  .  .  (t) 
under  a  given  stress  history  a.,  (t)  can  be  represented  by  the  following 
constitutive  equations 

ft  3a.  .(0  ft  3a.  At) 

eij   (t)  =  6ij     I      *1    {t'°  ~lt dK  +  J         h  {t<)       dt~    dE>     (4_1) 

o  ° 

In  this  equation,  e..,  a.,   represent  the  strain  and  stress  tensors 
respectively  and  each  consists  of  six  independent  components.  Three 
normal  strains  and  stresses  and  three  shear  strains  and  stresses:  akk 
represents  the  dilational  (confining)  stress;  and  5..  is  a  kronecker 
tensor,  which  equals  1  when  i=j?  and  o  when  ij'j.  Thus,  under  a 
uniaxial  stress  (unconfined)state,  a.,  =  a,,  =  a  and  a..    =  °/3\   the 
axial  creep  strain  e,,  (t)  becomes 

l(t)  =  I  J        [(J)1  (t"?)  +  ^     (t-')]  ^f~^  (4"2) 


el 
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or 

-   ut, 

(4-3) 


Jo 


e„(t)-  I   ^  (t-O^-dC 


J^t)  =  j  [<D1  (t)  +  3  d)2(t)] 


The  radial  creep  strain  (transverse  strain)  e??  =  e,~  (t)  becomes 


«o(t)  =  1  /t  ^(t-Fj  ^-dF  (4-4) 

^0 


'22v  '   3  /   Hv  w  9F, 


Under  a  confining  test  in  which  the  radial  stress  equals  P  and  axial 
stress  a,,  =  a'  +  P,  the  creep  strain  components  become 


Ell 


£22 


(t)  =  1  J      J^t-0  -^-  dC  +  y1  ^(t-0  ^f1  d?   (4-5) 

(t)  =  E33(t)  =  i  j    ^(t-o  g[gtat3Eia]  dc  (4-6) 


In  both  cases,  <J),  (t)  and  (J>?(t)  or  <£-,(t)  and  J,(t)  can  be  determined. 
When  the  axial  stress  is  a  constant  o(t)  =  a  ,   (4-3)  and  (4-4) 
become 

en(t)  =  J1(t)aQ  (4-7a) 

e22  =  £33(t)  =  *l(t)a0  (4"7b) 

Thus  the  axial  creep  compliance  J-(t)  can  be  determined  by  dividing 

the  axial  creep  strain  by  the  applied  stress. 

eii(t) 
J^t)  =-4 (4-8) 
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Using  the  definitions  of  Poisson's  ratio,   the  Poisson  ratio  v 
of  a  linear  viscoelastic  material  can  be  determined  as  follows 


E22     ¥t)qo      *l(t)  ,,  ., 


Equation  (4-9)  implies  that  the  Poisson's  ratio  for  a  linear  viscoelastic 
material,  in  general,  is  time  dependent.  However,  in  the  VESYS  I IM 
program,  v  is  considered  to  be  time-independent,  though  it  is  not 
necessarily  equal  to  one  half.* 

When  the  material  response  is  insensitive  to  time,  the  creep 
strain  e -, -,  ( t )  under  a  constant  stress  is  independent  of  time  or  nearly 
so,  then  the  creep  compliance  in  (4-8)  becomes  a  constant,  and  the 
reciprocal  is  the  elastic  modulus.  Further,  if  only  the  resilient 
strain  instead  of  total  strain  is  used,  then  the  modulus  is  called 
resil ient  modulus. 
4.4  CREEP  COMPLIANCE  OF  ASPHALT  CONCRETE 

The  mix  design,  gradation  of  aggregates,  density  and  void  ratio 
of  the  asphalt  concrete  samples  made  in  the  laboratory  should  simulate 
the  in-place  asphalt  concrete  material.  The  test  samples  are  usually 
cylindrical  in  shape  with  the  height  of  the  specimen  about  1-1/2  to 
2  times  the  diameter.  In  creep  tests,  a  constant  stress  is  applied  and 
strain  is  measured  as  a  function  of  time.  Normally  this  test  is 
conducted  at  constant  temperature.   The  creep  compliance  can  be 
determined  from  (4.8).  A  typical  result  of  creep  compliance  versus 
time  is  shown  in  Figure  4-1. 


■The  original  VESYS  II  set  the  Poisson's  ratio  equal  to  one-half. 
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Laboratory  test  results  of  creep  compliances  versus  time  are  input 

to  the  VESYS  I IM  program.  The  program  curve  fits  the  creep  compliances 

to  an  exponential  (Dirichlet)  series  of  the  form: 

N 
J  (t)   ]T   G1  exp  (-t6.)  (4-10) 

i=l 

where  G.  are  the  coefficients  of  the  Dirichlet  series  for  the  creep 

compliance;  6.  are  the  exponents  of  the  series,  they  are  generally 

chosen  to  be  one  order  of  magnitude  apart  from  each  other;  and  N  is 

the  number  of  terms  in  the  series. 

In  the  present  VESYS  I IM  program,  the  coefficient  of  variation 

of  J(t)  is  used  to  take  into  consideration  the  variability  of  the 

creep  compliance.  It  can  be  estimated  partially  from  a  set  of  laboratory 

creep  tests  results  and  partially  from  the  in-place  quality  control 

and  the  construction  method.  The  above  discussion  is  directly  applicable 

to  any  pavement  material  exhibiting  time  dependency.  This  includes  the 

bituminous  treated  base,  clay  soils.  Creep  of  soil  and  the  effect 

of  moisture  are   discussed  in  a  later  section  of  this  chapter. 

4.5  TEMPERATURE  EFFECT  ON  CREEP  COMPLIANCE  OF  ASPHALT  CONCRETE  AND  TIME- 
TEMPERATURE  SUPERPOSITION 

At  different  temperatures,  the  creep  compliance  and  relaxation 

modulus  of  asphalt  concrete  will  be  different  from  those  determined  at 

a  reference  temperature.  This  dependency  on  temperature  as  well  as 

on  time  of  creep  compliance  can  be  expressed  by  the  following  equation: 

J  =  J  (T,t)  (4.11) 
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where  T  is  the  temperature  and  t  is  the  time.  This  creep  compliance 
function  can  be  determined  in  principle,  from  a  proper  set  of  tests 
which  require  several  creep  tests  under  different  isothermal  temperature 
levels.  This  poses  a  problem  if  the  number  of  tests  required  are  yery 
large.  Furthermore,  even  if  this  approach  is  workable  the  analysis  of 
pavement  layer  systems  involving  transient  temperatures  becomes  very 
involved. 

Fortunately,  theoretical  and  experimental  results  for  many 
viscoelastic  materials  [21,  22]  and  experimental  results  on  asphalt 
concrete  [20,  24,  25]  indicate  the  effects  due  to  time  and  temperature 
can  be  combined  into  a  single  parameter  though  the  concept  of  the 
"time-temperature  superposition  principle"  which  implies  that  the 
following  relations  exist 


J  (T,t)  =  J  (TQ,0  (4.12) 

and 

C  =  t/aT(T)  (4-12a) 

where  t  is  the  actual  time  of  observation  measured  from  first  application 
of  load,  T  is  the  temperature  and  £  is  the  "reduced  time",  and  T  is 
the  reference  temperature.  The  reduced  time  is  related  to  the  real 
time  t  by  the  temperature  shift  factor  AT(t),  which  is  function  of 
temperature.  The  time-temperature  superposition  principle  cited 
above  states  that  the  effect  of  temperature  on  the  time  dependent 
mechanical  behavior  is  equivalent  to  a  stretching  (or  shrinking)  of 
the  real  time  for  temperature  above  (or  below)  the  reference.   In 
other  words,  the  behavior  of  material  at  high  temperature  and  high 
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strain  rate  is  similar  to  that  at  low  temperature  and  low  strain  rate. 
Materials  exhibiting  this  property  are  called  "Thermoheologically 
simple".  Thus  determination  of  the  temperature  shift  factor  AT(t)  as 
a  function  of  temperature  will  provide  the  necessary  information  for 
determination  of  the  reduced  time  £.  Figure  4-2  and  4-3  are  the  master 
creep  curve  and  shift  factor  versus  temperature  from  the  creep  compliances 
shown  in  Figure  4-1 . 

Therefore  (4.10)  becomes 

N 
J9t,(T)  =  £  G.  exp  (-c6i}  (4.13) 

i=l 

for  creep  compliance  at  temperatures  other  than  reference  temperature. 
4.6  RESILIENT  MODULUS  OF  AGGREGATE  BASE 

Aggregate  base  materials  are  found  to  exhibit  no  creep  under 
sustained  load  and  therefore  the  response  of  aggregates  under  loading 
can  best  be  represented  by  the  resilient  modulus.   The  resilient 
modulus  of  aggregates  is  usually  obtained  through  a  triaxial -compression 
repeated-load  tests.  Under  the  application  of  repeated  stresses 
confining  pressure  plus  axial  deviatoric  stress,  the  axial  strain 
response  is  measured.  The  resilient  modulus  is  defined  as  the  ratio  of 
a  repeatedly  applied  deviatoric  stress  to  the  recoverable  axial  strain. 
The  resilient  modulus  ER  of  aggregate  base  material  is  found  to  be 
nonlinear  function  of  mean  normal  stress, 

ER  =  K(a)n  (4-14) 

Figure  4-4  shows  the  effect  of  confining  stress  on  the  resilient 
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modulus  of  a  granular  base  material. 

In  order  to  select  a  resilient  modulus  for  an  aggregate  base  material 
for  use  in  the  layerd  system,  the  approximate  mean  normal  stress 
induced  in  the  aggregate  base  of  the  layered  system  when  subjected 
to  the  anticipated  external  load  has  to  be  known  as  priori.  This 
approximate  mean  normal  stress  can  be  obtained,  through  estimation 
or  through  iteration  between  the  estimation  of  mean  normal  stress, 
the  resilient  modulus  and  the  actual  response  of  the  system. 
4.7  CREEP  COMPLIANCE  AND  MOISTURE  EFFECT  OF  SUBGRADE 

The  creep  compliance  of  subgrade  materials  can  be  determined  in 
a  manner  similar  to  the  determination  of  the  creep  compliance  of  asphalt 
concrete.  The  only  difference  is  that  confining  compression  creep 
tests  are  usually  used  for  subgrade  materials  instead  of  unconfined 
compression  creep  tests. 

The  significance  of  time-dependence  of  creep  compliance  of  subgrade 
materials  depends  on  the  type  of  subgrade  material  and  moisture  content. 
Usually,  temperature  has  little  effect  on  the  creep  of  subgrade 
materials.  For  silts  and  clays  of  high  liquid  limit  and  high  plasticity, 
the  materials  exhibit  more  time-dependency  than  those  with  low  liquid 
limit,  and  coarse  grained  soils.  Moisture  content  in  soils  has  very 
significant  effect  on  the  creep  behavior  of  the  material.  This  has 
been  demonstrated  in  [26]  and  Figure  4-5. 

From  a  phenomenological  standpoint,  the  effect  of  moisture  and 
time  on  the  creep  of  subgrade  materials  can  be  combined  into  a  single 
parameter  using  the  concept  of  "time-moisture  superposition  principle" 
similar  to  the  "time-temperature  superposition  principle"  combining  the 


44 


TgS. 


CM 


O 


00 


r-i 


o 


o 

IS 


ITS 

CVJ 

<-\ 

o 
2 


o 
o 
o 

o 


o 
o 


o 

lO 


CO 

o 
o 

w 

CO 


o 


H 
EH 


H 

d 


o 
o 
o 


o 
o 


o 


Lf\ 


o 

ld 


CO 
EH 

En 


S 

o 

o 


EH 

CO  K 

H  EH 
O 

S  En 

CO 
CO  H 
D  En 
O 

M    H 
K  EH 

K 
O 

P-H 


CO 


EH  a 


£ 


o 


Eh 

3g 
1-^   K 
PL,   W 
Oh 


O 
O 

P-. 


in 
> 


H 


9" 


oi  x  aoNYiidwoo  daano 


effect  of  time  and  moisture  into  a  single  parameter  called  reduced  time, 
Hence  the  creep  compliance  obtained  at  different  moisture  levels  can 
be  reduced  to  the  following  form 


J(M,t)  =  J(M0,C)  (4-15) 


where 

M  =  moisture  content 

M  =  reference  moisture  content 
o 

t  =  real  time 
£  =  moisture  reduced  time 
and  the  real  time  and  the  reduced  time  are  related  by 

5  =  t/Am(M)  (4-16) 

where  A (M)  is  called  moisture  shift  factor  which  is  a  function  of 
moisture  content.  Inserting  (4.15)  (4.16)  into  (4.10)  yields  the 
following  expression  for  creep  compliance  at  different  moisture  levels 

N 
(M,t)  =  2   G.  exp  (-&.)  (4-17) 

i=l 

As  mentioned  above,  the  creep  behavior  of  subgrade  material  is  not 
dependent  on  temperature  so  long  as  the  temperature  is  above 
freezing  temperature.  However,  if  the  temperature  does  become  a 
significant  factor;  that  is,  the  creep  compliance  depends  on  both 
temperature  and  moisture,  then 


J  (T,  M,  t)  =  J  (TQ,  MQ,  c)  (4-18) 
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where  the  temperature  moisture  reduced  time  z,   equals 


<=  =  AT(t)\(M)  <4'19> 


and  the  moisture  shift  factor  A  (M)  and  temperature  shift  factor 
AT(t)  can  be  determined  from  (4.10)  and  (4.12a). 
4.8  FATIGUE  PROPERTIES  OF  ASPHALT  CONCRETE 

In  the  present  VESYS  I IM  system  conventional  fatigue  properties 
of  asphalt  concrete  are  used  to  predict  the  fatigue  cracking  of 
pavement  surface;  though  eventually  a  fracture  mechanics  approach  to 
the  determination  of  pavement  cracking  should  be  used.  Results  of  the 
fatigue  properties  of  asphalt  concrete  have  been  collected  by  many 
investigators  in  the  past  ten  years,  see  for  example  [24].  Typically 
fatigue  test  results  follow  a  power-law  function  of  the  tensile 
strain  amplitude.  Hence,  in  the  present  VESYS  I IM  the  following 
relation  is  used  to  represent  the  fatigue  life  (N^)  under  a  constant 
strain  amplitude  (Ae) 


K 
Nf  =  K]  (Ae)  (4-20) 


K,  and  Kg  are  coefficients  related  to  the  materials  characteristics 
which  can  be  determined  from  a  set  of  fatigue  tests.  In  general, 
these  coefficients  are  generally  sensitive  to  temperature  and  are 
also  statistical  in  nature.  The  dependency  of  K,  and  K?  on  temperature 
is  to  be  determined  from  further  experimental  data. 

The  coefficients  K,  and  K~  may  be  statistically  correlated.  The 
coefficient  of  correlation  could  vary  between  -1  and  +1.  This 
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coefficient  of  correlation  is  an  important  effect  on  the  analysis.  If 
this  coefficient  is  (-1)  it  means  that  an  increase  of  K,  corresponds 
necessarily  to  a  decrease  of  K2  while  if  the  coefficient  is  zero  it 
means  that  K,  and  K2  are  statistically  independent  of  each  other. 
These  coefficients  are  given  by  their  means,  variances  and  coefficient 
of  correlation. 

The  cumulative  damage  law  based  on  Miner's  rule  as  shown  in  the 
following: 


m 
Damage  =  £  $-  (4-21) 

i=l    f  n 

where 

n  =  number  of  repetions  of  the  strain  designated  by  the 
range  i 


and 


Nf  =  number  of  repetions  in  the  range  i  required  to  produce 
failure. 
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CHAPTER  5 
PRIMARY  RESPONSE  MODEL 

5.1  INTRODUCTION 

The  VESYS  II  pavement  design  framework,  as  noted  previously, 
solves  for  the  stresses,  strains  and  deflections  in  a  three-layered, 
incompressible,  viscoelastic  system  using  a  convolution  integral 
formulation.  The  solution  is  only  valid  for  three  layers  which  are 
all  incompressible.  In  the  course  of  extending  the  solution  to  five 
layers  and  incorporating  an  approximate  closed  form  probabilistic 
solution,  it  was  decided  that  greater  efficiencies  could  be  obtained 
by  discarding  this  primary  response  model  and  developing  a  general 
layered  solution  which  is  valid  for  an  arbitrary  number  of  layers 
(the  limiting  factor  being  computer  storage  capabilities)  and  for  any 
value  of  Poisson's  ratio.  In  addition,  due  to  the  time  scale  of 
loading  for  pavement  systems  the  quasi-elastic  approximations  to  the 
viscoelastic  response  is  employed  to  further  reduce  computation  time. 
It  is  demonstrated  that  acceptable  solutions  are  obtained  with 
this  approximation. 

An  approximate  closed-form  probabilistic  solution  is  obtained 
using  the  first  order,  second  moment  theory. 
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5.2  QUASI-ELASTIC  ANALYSIS 

The  quasi-elastic  method  of  analysis  is  based  upon  the  work  of 
Schapery  [27]  and  is  applicable  to  viscoelastic  problems  in  which  the 
derivative  of  the  time-dependent  solution  with  respect  to  logarithmic 
time,  log  t,  is  a  slowly  varying  function  of  log  t.  The  method 
involves  replacing  elastic  moduli  (or  compliances)  by  instantaneous 
values  of  the  relaxation  moduli,  and,  as  necessary,  making  a  finite 
difference  approximation  to  the  Boltzmann  Superposition  integral 
of  the  loading  function. 

As  an  illustration  of  the  procedure  of  performing  a  quasi-elastic 
analysis,  the  steps  followed  in  obtaining  a  solution  to  a  problem  for 
a  single  layer  are  outlined  below.  The  same  procedure  is  followed 
for  any  number  of  layers. 

1)  Determine  a  value  for  the  layer  temperature 

2)  Subdivide  the  loading  curve  into  a  sequence  of  time 
intervals  in  such  a  manner  that  the  curve  is  reasonably 
well  approximated  by  a  series  of  step  functions 

P(t)  =0        ,  t  =  0 

P(t)  =  APj        ,  t,  <  t  <  T2 

P(t)  =  APi  +  AP2   ,  T2  <  t  <  T3 


P(t)  =  APx  +  ...  +  APk   ,   tk  <  t  <  xk+1 

3)   Perform  a  stress  analysis  for  a  unit  load  (P(t)  =  1  psi) 
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for  each  of  the  time  increments,  x.  ,  by  some  appropriate 
means  using  as  an  elastic  compliance  the  value  of  the  creep 
compliance  corresponding  to  the  reduced  times  T./aT. 
4)   Denoting  the  components  of  the  displacment  vector  at  the 
point  x.  (k  =  1,  2,  3)  and  time  obtained  from  the  analysis 
carried  out  in  Step  3  above  by  w. ,  the  components  w. 
of  displacement  corresponding  to  the  step-wise  loading  may 
be  determined  from  the  following  set  of  equations* 

W.(xk,t)  =  AP1W.(s|<,t-T1)  ,   X!  <  t  <  T2 

w.(xk,t)  =  AP1W.(xk,t-T2) 

+  PiW.(x|c,t-x2)  ,  x2  <  t  <  x3 


N+l 
w.(xk,t)  -I  APnWi(xk,t-xn) 


The  values,  W.,  can  be  found  from  computer  solutions  at  a  few  log  t 
points  plus  linear  interpolations  in  log  t. 

Although  the  above  illustration  was  for  a  single  layer,  and 
statistical  variations  of  the  loading  history  and  materials  properties 
were  excluded,  inclusion  of  these  effects  is  straight-forward,  and 
material  property  variability  is  included  in  the  computer  code 
developed. 


*  If  the  limits  N  -*  °°  and  AP.  -»•  0  x  <  t  <  t  +,  are  taken  (corresponding 

to  a  continuously  varying  pressure  P(t)),  this  set  of  equations  become 
the  Boltzman  Superposition  Integeral 
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Figures  5-1  and  5-2  on  Pages  53  and  54,  respectively  compare  MIT 
published  data  for  a  three-layered  viscoelastic  system  [15]  with  the 
solution  for  stress  and  deflection  obtained  using  discrete  moduli  values 
and  the  N-layer  program.  The  excellent  agreement  obtained  indicates 
the  applicability  of  the  quasi-elastic  approach. 
5.3  ELASTIC  N-LAYER  SOLUTION 

An  idealized  layered  system  can  be  represented  by  multilayers  of 
different  thicknesses  and  elastic  material  characteristics  as  shown 
in  Figure  5-3,  in  which  a  normal  load  of  intensity  q  is  distributed 
over  a  circular  area  of  radius  a. 

The  layer  materials  are  assumed  to  be  linearly  elastic,  isotropic, 
and  weightless.  The  layers  are  assumed  to  be  infinite  in  horizontal 
extent  and  the  lower  layer  (n-th  layer)  is  considered  to  be  semi- 
infinite  in  the  vertical  direction. 


»■       i 


E 


'21 


■n-1 


E1,v1 


E2,v2 


En-l*Vl 


i 


1 


En'vn 


z   ? 
Figure  5-3      An  n-Layer  Pavement  System 
52 


Q 

ui 

u 

m 

3 

Q 

ai 

> 

u 

to 

2 

h- 

02 

Q 

UJ 

o 

•J 

O 

O 

m 

q: 

p 

J 

«o 

Lf) 

o 

2 

b/'jo 


53 


o 
o 
o 


o 
o 


o 
o 

ro 


O 
O 


UJ 


o 

lT> 

a 

UJ 

u 

Q 

o 

U4 

ro 

in 

> 

o 

* 

t-H 

UJ 

H 

2 

h- 

O 

Q 

w 
a 

UJ 

j 

O 

< 

D 

u 

Q 

UJ 

i— i 

H 
OS 

m 

q: 

> 
Cvl 

ro 

1 

a 

s- 

3 

b-(*)Sl/M 


54 


In  order  to  solve  for  the  components  of  stresses  and  strains  at 
any  point  of  the  structure,  linear  small  deformation  elasticity  theory 
is  employed.  A  stress  function  4>(r,z)  which  satisfies  the  governing 
differential  equation  is  determined  for  each  layer. 


Vcj)  =  0 


(5-1) 


In  cylindrical  coordinates,  (Figure  5-4)  the  equations  of 
equilibrium  have  the  form 

3a    3t  _  a   -afl 
r  +   rz  +  ^L_e    0 

3r    3z      r 


(5-2) 


3a    3t    t 
z  +  _zr  +  ^r  =  0 


3z    3r 


(5-3) 


where  a  ,  a,,  t  ,  x.  xQ  .  and  x  ~  are  components  of  the  stress  tensor 
r   z   rz   zr   or      ro 

with 

Tre    =    Ter    =  ° 


*3fcs£±7 


Figure  5-4  Stress  State  of  an  Element  in  Pavement  System 
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The  equations  of  compatibility  in  terms  of  stresses  are  given 
by 

r  or 

r 

r 

where 


v2  ,  Ll»  +aL  (5-8) 

v4  =  vV 


0     ■    °r+     °Q+    °z 


The  stresses  and  displacements  are  given  by 


.  3 


p.-  A] 


ap  =  |j|vV^-  if 


P«-H*] 


ae  *  If  I 


=  fj  [(2-v)v2<j,  -i-|J  (5-n) 


°z 
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=  !f[0-v)v< 


xr2  =  -|(l-v)V2,-^| 


(5-12) 


u  -  JfL  |!|-  (5-13) 

w  =  i+v  r2{1.v)v2^  .  a!|j  (5_14) 


v  *  0  (5-15) 

where 

u  is  the  radial  displacement 
v  is  the  tangential  displacement,  and 
w  is  the  vertical  displacement 
The  continuity  conditions  at  the  i-th  interface  are  given  by 


(5-16) 
(5-17) 
(5-18) 
(5-19) 


Equations  (5-1 G  to  5-19)  state  that  each  layer  at  the  i-th  interface  is 
considered  to  be  "rough",  that  is  the  components  of  stress  a  and 
t   as  well  as  the  components  of  displacement  u  and  w  are  continuous 
at  z  =  hi  (i  -  1,  2,  ...,  n-1). 


1 

1+1 

az 

i 
Trz  " 

i+1 

T 

rz 

u1  - 

ui+1 

i 
w   * 

w1+1 
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Since  the  region  of  interest  is  infinite  in  the  r  direction,  the 
biharmonic  equation  is  readily  solved  using  Hankel  Transforms  [28] 

i(m,z)  =    r<fr(r,z)  JQ(mr)dr  (5-20) 

The  inverse  function  is 

<f>(r,z)  =  J  m$(m,z)  JQ(mr)dm  (5-21) 

For  partial  differential  V24>,  it  follows  that 
2 


and 


\     =     rm|i|-  m2$|j0(mr)dni  (5-22) 


00  o 

I    rv24>  JQ(mr)dr    =     (^  -  m2)   fr^  JQ(mr)dr 

J     r74d>  J  (mr)dr    =     (^  -  m2)2   f  r<J>  J  (mr)dr  (5-23) 

Jo  u  dz^  Jo  ° 

Thus  instead  of  solving  biharmonic  equation  (5-1 ),  one  is  led  to  the 
solution  of  the  ordinary  differential  equation 

[^T*m2|  *(m'z)  =  °  (5'24) 


where  $  is  defined  in  Equation  (5-20).  The  solution  of  Equation  (5-24) 
leads  to  the  expression 

$(n.z)  =  (A+Bz)emz  '  +  (C+Dz)e"mz  (5-25) 
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For  the  i-th  layer  of  the  system,  Equation  (5-25)  can  be 
written 

^(m.z)  *     (A.+BjZ)emz  +  (C.+D^e""12  (5-26) 


where  A.,  B.,  C  and  D.,  are  constants  of  integration,  to  be  determined 
in  such  a  way  that  the  boundary  conditions     are  satisfied. 

Now,  multiplying  <j>.   by  J   (mr)   in  Equation   (5-26),  and  defining 
4>.(m,r,z)  as  the  product  of  <j)(m,z)  J   (mr) 

^(m.r.z)     =    JQ(mr)     [(A1+B1z)emz  +  (C^D^e-"12]         (5-27) 

it  can  be  verified  that  Equation   (5-1)   is  satisfied  by  the  function 
(^(m.r.z). 

Making  use  of  Equation  (5-27),   Equations   (5-9)through  (5-14) 
become 

o\  =  m2J0(mr){(l-2v.)[Biemz  +  D^""12]  -  m[(A1+B.z)e'nz  -  (C^zfe-™2]}    (5-28) 

~1     _  _2,   ,„uo..  rD  jnz       n     -mzn   .  ml-/a   ,„  ,lflmz   .    (r   .n  ,\  -mz- 


Tpz  =  nrJ1(mr){2vi[Biemz  -  D^"™*]  +  m[(A.+B1z)e"M;  +  (C^zje  ""]}  (5-29) 

P  =  m2Jn(mr){(H2vi)[Bien,z+Die-mz]  +  m^B^e*2  -  (C.+D^Je^2]} 


2  ^2  {B.emz+Die-mz  ♦  m[(A.+BiZ)emz  -  (C^D^e-"12]}  (5-30) 


■m        mr 


5J  ■  2vim2J0(mr)[B.emz+Dje-mz]  +  n,2  ^^Of^M^)*™ 


-  (Ci+Diz)e-mz]}  (5-31) 

~f  .  ^mJl(mr){B1en,z+Die-mz+m[(Ai+BiZ)emz  -  (C^O^Je-"12]}  (5-32) 
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w1  »  -r1  J>r)m{2(l-2v.)[B.ernz-D.e"mz]  -mC(A,+B.z)emz+(C,+D.z)e'niz]}   (5-33) 

tjO  111  11  11 


Substituting  Equations  (5-28)  to  (5-31 )  into  Equations  (5-16)  and 
(5-19)  and  rewriting, 

5 


S^z) 


*z 

^rz 

-i 
u 

w1 


=  K(vi,Ei)  M(z,vi)  D(z) 


'i 


B 


where 


L°'J 


(5-34) 


KCv^E.)  = 


-m  JQ(mr) 


m  J, (mr) 


1+v, 


-  mJ, (mr) 


1+v, 


-  mJ  (mr) 


E.  -o 


(5-35) 


and 


M(z,vi) 


me 


1 
1 
1 

1 

mz 


mz+2^-1 
mz+2v. 


mz+1 


mzHv^-2 


-mz+Zv^l 
mz-2v. 
-mz+1 


mz-^^2 


mz 


D(z)  = 


me 


-mz 


(5-36) 


(5-37) 


-mz 
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Continuity  conditions  at  the  i-th  interface  can  be  expressed  by 
Sjthj-O)  =  S1+l  (h^  (5.38) 


or 


Ai 


D(hi)XiD(hi) 
4(vrl) 


"1+1 
Bi+1 
CH1 
Di+1 


(5-39) 


and 


Ai 


n-l     Dfh.JX.DthJ 

j=1  J 


"o  ' 

0 

Cn 

(5-40) 


Xi 


(4V.-3J-L         "^"WW       (2mh.+4vrl)  p(h.  ,v.  ,v1+l )  ' 

+Li^"hi'VrV     V+i)      -hP«hr'vi+Tvi) 

2(1,-1)       (2"V4Vl+1) 
(Lrl) 


0  L.(4vi+1-3)-l 

-p(-h-,v     v       )  q(h     v     v       ) 

;-r-i-)  ^i.^%)  ^v3)-Li  .h;(h;,v;;;,Vi 

[zo-l.)         I^i^i+r1)  °        Li<4vr3H  J 


(2mh.+4v      -1) 
(1-L|)     1+1 


where      p(h,v,u)  =  2(mh)     +  4mh(v-y)  +1   -  8vu  +  2y 
q(h,v,y)  =  2mh(2v-l)  ■     (1   +  8uv  -  6y) 

L 


E1       1+Vl 


'i 


Ei+1     1+vi 


(5-41) 
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So  far  the  boundary  .conditions  at  the  interface  have  been  con- 
sidered. At  the  loading  surface  a  vertical  uniform  loading  may  be 
assumed  in  the  form 


p(m)Jo(mr) 


(5-42) 


By  imposing  a 


=  -p(m)J  (mr),  Equation  (5-28)  reduces  to 


mA, 


+  (2v1-1)B1hbC1+  (2v1-l)"D1  =  ^=  K(m)    (5"43) 

m 


By  imposing  x 


1 
rz 


=  0,  Equation  (5-29)  reduces  to 

z=o  ^  x         ' 


mA]  +  2v1B1  +  mC1   -  2\>}Q    =  o 


(5-44) 


Equations   (5-43)  and  (5-44)  can  be  written  as 


K(m) 


m        2v,  m      -2v1 


m        2v,-l       -m       2v1-l 


Al 

Bl 

cl 

[Dl 

(5-45) 


or 


K(m) 


.[.] 


n_1   D(h.r1X.D(h.) 


j-1 


4^ry 


"o 

0 

Cn 

(5-46) 


where 


Q    = 


m     2v  m      -2v( 


m     2v-l     -m      2v-l 


(5-47) 
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The  constants  C  and  D  are  therefore  determined.  Knowing  C  and 
n     n  n 

D  ,  the  values  of  A.,B.,C.  and  D.  are  obtainable  from  Equation  (5-40) 
since  Equations  (5-9)  and  (5-14)  are  linear  functions  in  <f>,  the  stress 
function 


<fr(r,z)  ■    $(m,r, 

in 


z)dm 


(5-48) 


will  solve  Equations  (5-9)  and  (5-14)  when  subjected  to  a  load 


p(m)J  (mr)dm 


(5-49) 


The  inverse  formula  for  the  Hankel  Transform  is 


\  [f (r-0)  +  f (r+0)]  =    m[r  f(r)JQ(mr)dr]  J0(mr)dm     (5-50) 


Therefore  p(m)  can  be  expressed  as 

P(m)  =  m  f  r  f(r)J  (mr)dr 
*  a  J^ma) 


where 


fir)    =/l 


2   2 
r  <  a 


2   2 
r  *  a 


2   2 
r  >  a 


(5-51) 


The  stresses  and  displacements  at  any  depth  in  the  i-th  layer  are 
then  given  by 


i>, 


az  a  a    J^  (ma  )a^dm 


(5-52) 


Vz  =  a 


I  J^majx^dm 


^i 


(5-53) 
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oj.    -    a  [     J^majojdm  (5-54) 

56     "     a  |     J1(ma)ajdm  (5.55) 

u1     -    a  J     J-,(ma)u  dm  (5-56) 


w1     ■    a 


r00 


>-1 


J,(ma)w  dm  (5-57) 

0 


5.4  PROBABILISTIC  N-LAYER  SOLUTION 

Since  the  external  loading,  environmental  conditions  and  internal 
system  parameters  of  a  pavement  structure  are  generally  random  quantities, 
the  underlying  reliability  theories  must  employ  probablistic  concepts. 
However,  most  structural  analyses  are  deterministic  with  uncertainties 
reflected  in  design  safety  factors.  These  safety  factors  are  not 
normally  based  on  mathematical  models  describing  the  interplay  between 
the  various  uncertainties. 

Generally  speaking,  it  is  impossible  to  determine  the  character- 
istics of  an  uncertain  quantity  given  as  a  nonlinear  function  of  other 
uncertain  quantities  soley  by  the  characteristics  of  the  set  of  these 
uncertain  quantities.  Such  a  determination  is  possible  only  if  the 
distribution  function  of  the  set  is  known  and  satisfies  certain 
integrability  requirements  [29]. 

Solving  stochastic  differential  equations  is  one  currently  prev- 
alent method  for  structural  problems  with  simple  geometries  [30].  In 
applying  this  method,  however,  not  only  must  the  problems  be  limited 
to  simple  geometries  and  material  properties,  but  simplifying  assump- 
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tions  must  also  be  made  for  the  random  variables,  such  as  assuming 
normal  distributions,  because  of  the  complexity  of  the  differential 
operators  involved.  Moreover,  the  mathematics  involved  is  often  so 
complicated  that  the  problem  becomes  intractable  [31].  An  appropriate 
probabilistic  model  was  developed  in  this  program  since  the  application 
of  full  probability  theory  has  several  practical  limitations  [4]: 

1)  The  mathematics  becomes  complicated. 

2)  Available  data  and  information  about  uncertainties  are  often 
too  scarce  for  practical  implementation. 

Additional  practical  considerations  which  support  development  of 
approximate  probabilistic  models  include  [4]. 

1)  Near  the  mean,  probabilities  are  insensitive  to  the  distri- 
bution type. 

2)  At  the  tails,  data  and  models  are  too  weak  to  predict  reliable 
absolute  probabilities;  but  many  distributions  give  similar 
relative  changes  in  probability  per  unit  increase  in  the 
standardization  abscissa. 

3)  When  expected  cost  optimization  is  used  to  determine  engineering 
design,  first  and  second  moments  are  sufficient  if  the  cost 
function  is  quadratic. 

4)  Engineering  design  and  cost  are  yery   insensitive  to  small 
changes  in  probabilities  (the  relationship  usually  being 
approximately  logarithmic). 

The  approach  followed  in  the  program  is  based  on  Cornell's  first- 
order,  second  moment  theory  [  4  ].  The  term  "first-order"  indicates 
the  application  of  linearization  to  nonlinear  equations  and  the  term 
"second-moment"  refers  to  the  fact  that  only  the  first  two  moments  of 
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the  random  variables  involved  are  needed. 

This  approach  is  believed  to  be  superior  to  other  approaches,  such 
as  the  use  of  maximum  likelihood  methods  or  moment  generating  functions, 
since  these  approaches  require  knowledge  of  the  distributions  of  the 
random  variables  [32]  whereas  the  first-order,  second-moment  theory 
does  not  require  knowledge  of  the  distributions  of  the  random  variables. 
5.4.1  First  Order,  Second  Moment  Theory 

In  the  application  of  a  first-order,  second-moment  approximation, 
the  term  "first-order:  indicates  linearization  of  nonlinear  equations 
and  the  term  "second-moment"  refers  to  the  fact  that  only  the  first  two 
moments  of  the  random  variables  are  needed.  For  example,  some  random 
variables  X  can  be  characterized  by  their  mean,  X,  and  variance, 
Var  [X]  =  o2  (or  equavalently  their  standard  deviation,  aY,  or  coefficient 

X  A 

of  variation  V»  =  oJX).     The  stochastic  dependence  between  a  random 
variable  X  and  a  random  variable  Y  is  measured  by  their  covariance, 
Cov[X,Y],  or  equivalently  by  their  correlation  coefficient 

=  Cov  [X,Y]  (5_58) 

PXY     oxOy  {0  00) 

Similarly,  a  random  process  X(t)  is  defined  by  its  mean  function 
X"(t)  and  is  covariance  function  Cov[X(t),  X(s)]  (possibly  cross- 
covariance  function). 

Linearization  of  a  nonlinear  equation  is  illustrated  in  the 
following  example.  Let  Y  be  a  continuous  function  of  X, 

Y  =  g(X)  (5-59) 
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where  X  and  Y  are  random  variables.  The  mean  of  X  is  denoted  by  X. 
Expanding  g(X)  about  the  mean  X  in  a  Taylor  series  expansion,  then 

9(X)  ■  g(X)  +  (X-X)  M*l  I  _  +  (x-x)2  Im.   ,    +      (5-6O) 

x=x       dx2   x*x~ 

Retaining  only  the  first  two  terms  of  the  expansion  and  taking  the 
expectation  of  both  sides  then  yields  the  approximation 


E[Y]  =  E[g(X)]  =  g(E(X)) 


(5-61) 


for  E[g(x)]  since  E(X-X)  =  0.  Similarly,  retaining  the  same  terms, 
the  approximation  for  the  variance  is  given  by 


Var[Y]  =  Var[g(X)]  =  Var[X]  (^&  \^_j' 


Since  g(X)  is  a  constant, 


Var[g(X)]  =  0 


and 


Var 


/  dg(x) 
\     dX 


«) 


(X-X) 


=  Var[X]  (MIL  | 


(m*±  I    V 

^dx     'x=x) 


(5-62) 


(5-63) 


(5-64) 


Using  this  same  principle  for  a  variable,  say  Z,  which  is  a 

continuous  function  of  several  variables,  X, ,  X0,  ...,  X  ,  -  whose 

I   2       n 

means,  variances,  and  covariances  are  known  -  the  first  order  approx- 
imation of  the  mean  and  variance  of  Z  can  be  obtained  by  expanding  the 
function  g(X, ,  X«,  ...»  X  )  in  a  multiple-variable  Taylor  series, 
about  the  mean  values  of  the  X's.  Keeping  only  the  linear  terms  in 
the  deviation  from  the  means  and  applying  the  same  linear  transfor- 
mation rules  for  the  mean  and  variance  yields  the  approximations 


E[Z]  =  Z=  g  (Xr  X2,  ....  Xn) 
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(5-65) 


and 

Var[Z]  =    i     (&-  |     _\Z  Var[X.] 
i=l    \  oAi     X=X  1  1 

+  I  j,   (%  'x=x)  fa  W)  C°V[X^  V         (5"66) 

in  which  the  partial  derivatives  are  evaluated  at  the  mean  values 
of  X.  and  X.,  i.e.,  X"  ;  and  X".. 

If  the  random  variables  X.  and  X.  are  independent,  then  Equation 
(5-66)  simplifies  to 

Var[Z]  -  S  / IS-  I   V   Var  [Xj  (5-67) 

1=1  I  8Xi  X=X J  1 

Since  Cov[X.,  X.]  =  0  for  X.  and  X.  uncorrelated. 
l   J         i     J 

Many  details  described  by  full  probability  theory  are  not  described 
by  the  first  order  second-moment  theory;  however,  the  main  idea  is  to 
extend  a  deterministic  description  to  obtain  the  essential  behavior  of 
the  uncertainty  interplay  with  as  few  new  concepts  as  possible.  The 
results  obtained  are  independent  of  any  assumptoins  about  the  form  of 
the  underlying  probability  distributions.  If  appropriate  assumptions 
can  be  adopted,  the  description  of  the  moments  of  the  random  variables 
may  be  extended  to  some  probability  laws. 
5.4.2  Closed  Form  Probabilistic  Solution 

The  development  of  an  approximate  closed  form  probabilistic  solution 
for  the  stresses  and  deflections  given  in  Section  5.3  proceeds  in  the 
following  manner.  Let 

si  ■  H^r   E2 "V  (5"68) 

68 


where  S.  represents  either  the  desired  stress  or  displacement  at  i-th 
layer;  E, ,  E?,  ...,  E  are  the  elastic  moduli,  which  are  assumed  to 
be  the  only  random  variables  in  the  system.  The  E.  are  further 
assumed  to  be  independent.  Then,  from  the  previous  section  the  expected 
value  of  S.  may  be  approximated  by  the  expression 


E[S.] 


•i^i'  V 


i 


a2g« 


•  V  +t. 


j  =  l  8E 


J 


(5-69) 


The  variance  of  S .  is  given  by 


(5-70) 


The  term  g^E,,  E„,  ...,  E  )  represents  the  mean  value  of  S.  ,  i.e., 

the  value  of  S.  obtained  using  mean  moduli  values  E..  The  terms 
l  i 


*2 


3E. 
J 


2 


F,  °f 


and 


are  the  2   and  1   partial  derivatives  of  S.  with  respect  to  E. 

evaluated  at  E-   and  multiplied  by  ol     respectively.  The  E",  and 
J  tj  J 

Op  are  assumed  to  have  been  determined  from  an  experimental 
characterization  program. 
From  Equation  (5-34) 
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3K(v.',    E.) 


-  M(z.,    v.)    Of:) 


"a." 

1 

B. 
i 

C. 
i 

D. 
i_ 

+   KCvk,    E.)   M(z. 


Vi)    D(z) 


3A. 


3E. 
J 

3B. 
l 


3E. 
3C" 


3 


3E. 
J 

3D. 

l 

3E. 
J 


(5-71) 


3E 


3^S, 


3E 


3   K(v.,    E.) 

l        i 

3E/ 


M(z.,   v±)   D(z) 


A. 

l 


+   2   K(v.,    E.)   M(z.,   y/)    D(z) 


+   K(v.,    E.)    M(zif   v.)   D(z) 


3  A, 


r 


9Ej 

3B. 

3E. 
3 

3C. 

i 

3Ej 

3D. 

l 

3  E. 
J  . 

32A. 
1 

3E.2 

•7 

3"B. 
l 

3E/ 

J 

3-C. 
i 

■? 
3E." 
J 

3  -n. 
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(5-72) 


3K(v.,    E.j 

air 


-L- v— v    5. .   m  J.    mr 
Ei."  1J  l 


+(1!V)   6..   m  J      (mr) 
E.2  1J  ° 


(5-73) 


|K(V.,   E.) 

v   l        l 


and 


,.  .   j1 


+2(l+v)    -  T    ,     . 

*-   ,   -  6.  .   m  J,  (mr) 

c     3  !J  X 


E. 


when   i=j 
when  i^j 


m^lfi..   m  J    (mr) 

E.3  1J  ° 

3 


(5-74) 


Differentiate  Equation   (5-46)  and  rearranqinq  the  terms 


n-1  D_1(h.)    X.    D(h   ) 

q   n  — tP- — l — l 

•    ,  4(v.    -    1) 

3=1  3 


0 

0 

3C 


J 


3E. 
3D" 

T 

3E 


-1         *h 

n-1   D      ^  JeT0^    n-1   D_1(h.)    X.    D(h.) 
I     3 n     3 3 3_ 

1=1 


4(v£  -   1) 


n 

3=1 


4  (v.    -    1)   , 
J 


3X. 
where  3E  is  defined  in  Equatioi  (5-76) 
j 


(5-75) 
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From  Equation  (5-40)  we  obtain 
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Since  C   ,  D   .  A.,  B. ,  C,   D.   are  determined  using  the  mean  values, 

_   3C   3D   3A.   3Bi   3C .   3D. 

E.j,  -r=— '  -tf— '  -rp- '   -gr— '  gr— '  gp—  (i=l>  .  ..>  n-1)  may  be  determined  using 

J  iJ  vJ  J  J  «J 

Equations  (5-75)  and  (5-77).  The  second  partial  derivatives  of  S. 
with  respect  to  E-  are  obtained  from  Equations  (5-75)  and  (5-77): 
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Now  that 


3A.   3B.   32C.   32D. 


9E.'  3E.  3En 

J    J    J     J 


t     ^p  2  are  known,  these  nay  be  substituted  into 


Equations  (5-71)  and  (5-72)  and  then  (5-69)  and  (5-70)  to  obtain 
E[S.]  and  Var[S.]. 

In  evaluating  the  computer  code  for  the  means,  pavement  response 
it  was  found  that  the  terns  involving  the  second  deviations  were  not 
significant  and  hence  they  were  deleted  from  the  code. 

To  obtain  the  solutionfor  a  circular  load,  the  stresses  or  dis- 
placement must  be  integrated  from  zero  to  infinity  with  respect  to 
m  and  multiplied  by  a.  Hence 
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0         x       J 
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(5-81) 
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where  J, (ma).,  S . .  denote  J, (ma),  S.  evaluated  at  m  =  m. 
Similarly, 
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5.5  EXAMPLE  SOLUTIONS 

The  probabilistic  responses  for  surface  deflection,  strain  or  the 
bottom  of  the  surface  layer  and  stresses  in  the  surface  layer  of  a 
three-layer  pavement  system  subjected  to  a  circular  wheel  load 
are  illustrated  in  Figures  5-5,  5-6  and  5-7.  The  molulus  of  the  first 
layer  is  assumed  to  have  a  coefficient  of  variation  of  0.25  while  the 
properties  of  the  other  two  layers  are  deterministic.  It  may  be 
observed  that  the  variability  of  the  responses  are  quite  different. 
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CHAPTER  6 
RUT  DEPTH  PREDICITON  MODEL 

6.1  INTRODUCTION 

The  design  framework  proposed  by  Moavenzadeh  [1]  and  continued 
in  VESYS  IIM  is  shown  schematically  in  Figure  6-1.  It  consists  of 
five  components  which,  when  integrated  together,  result  in  a  predictive 
pavement  design  system.  Each  component  consists  of  several  elements 
or  subsystems  which  form  the  rational  basis  for  the  design  pro  edure. 
Among  the  five  subsystems  shown  in  Figure  6-1,  the  structural  sub- 
system is  perhaps  the  most  important  one.  As  shown  in  Figure  6-2, 
the  structural  subsystem  consists  of  three  models  -  the  primary  response 
model,  the  damage  model  and  the  serviciability  model. 

The  primary  response  model  accepts  the  basic  material  properties, 
(i.e.,  creep  compliance,  Poisson's  ratio),  pavement  system  geometry 
(i.e.,  layer  thickness),  traffic  loading  (e.g.,  wheel  load,  tire 
pressure,  frequency  or  distribution  of  load  applications,  etc.)  and 
environmental  variables  (e.g.,  temperature,  moisture  level,  etc.), 
and  computes  the  stresses,  strains  and  deflections  at  every   point  in 
the  pavement  system.  The  primary  response  model  in  VESYS  IIM  accepts 
stochastic  moduli  or  compliances  in  order  to  account  for  the  inherent 
variability  of  the  input  variables  and  the  pavement  system. 
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83 


CO 

o 

D_ 
CO 

LU  — I 
Cd    LU 

>-  o 
or:  :c 


ft 


LU 

CO 

CD 

o 

Q_ 

1 — 1 

CO 

1 — 

LU 

CO 

zz. 

c_> 

Ctl 

CO 

1— 1 

LU 

LU 

<c 

1 

>- 

en 

en 

u_ 

GZ. 

h— 

l — 

LlJ 

<z 

CO 

CO 

Q 

■ — 

en 

Q_ 

CD 


CO 
ST. 


LU  CD 


T 


>- 

LU  h- 
h-'LU 

co  sr: 

>-  CD 

CO  LU 

CD 


CD  >- 

~  CD 
CT  I — 

<n  co 

CD  — 


I 


CO 

>- 

CO 
CO 


CD 

h- 

CO 


84 


The  probabilistic  solution  for  the  primary  responses  i.e.,  the 
stresses,  strains  and  deflections,  obtained  from  the  primary  response 
model  and  the  material  failure  properties  are  inputs  to  the  damage 
model,  which  predicts  pavement  distress.  To  date,  fatigue  cracking, 
rutting  (permanent  deformation),  and  roughness  have  been  considered  to 
be  the  msot  significant  modes  of  pavement  distress.  The  pavement 
distress  modes  are  calculated  separately  in  three  submodels  even  though 
there  is  some  interaction  between  them. 

The  calculated  pavement  distresses  are  input  to  the  serviceability 
model  to  define  the  life  expectancy  of  the  pavement  system. 

Pavement  rutting*  is  the  result  of  a  channelized  traffic  causing 
differential  surface  deformation  under  the  areas  of  intensive  load 
applications  in  the  wheel  paths.  The  amount  of  pavement  rutting 
depends  upon  the  following  parameters: 

1)  The  distribution  of  traffic  loads,  particularly  the  transverse 
distribution  of  traffic  loads,  i.e.,  channelized  traffic; 

2)  The  stresses  induced  in  the  pavement  system,  which  depend 
upon  the  response  characteristics  of  the  layer  materials; 

3)  The  permanent  strain  induced  as  the  result  of  the  stresses 
developed  in  (2).  This  permanent  strain  depend  on  the  permanent 
deformation  characteristics  of  the  layer  materials. 

Although  satisfactory  means  for  rut  depth  prediction  have  not 
yet  been  developed,  the  mechanism  of  rutting  under  repeated  traffic 
loading  is  relatively  well  know  from  various  field  and  laboratory 
investigations.  In  the  AASHO  Road  Test  [33]  an  extensive  study  was 
made  of  the  rutting  phenomenon.  Two  conclusions  were  drawn  from  this 


*  "rutting"  as  used  in  this  report  refers  to  the  depressions  under  a  wheel 
path  to  repeated  wheel  loads  which  cause  consolidation  and/or  lateral 
movement  in  the  pavement  layers.  With  this  definition,  other  mechanisms 
which  may  cause  rutting,  such  as  abrasion,  and  war  due  to  studded  tires 
etc. ,  are  excluded.  85 


study  which  revealed  two  very   improtant  aspects  of  the  mechanism  of 
rutting: 

1.  Rutting  of  the  asphalt  pavement  surface  can  be  attributed 
to  the  change  of  thickness  of  all  layers  of  the  pavement 
system  (see  Figure  6-3). 

2.  Rutting  is  due  primarily  to  lateral  displacement  of  the 
pavement  materials,  rather  than  from  densification. 

These  findings  have  been  confirmed  by  other  field  test  results  [34]. 

Realizing  that  many  varieties  of  materials  are  used  for  paving 
materials,  and  that  the  characteristics  of  the  materials'  response 
properties  and  permanent  deforamtion  characteristics  are  not  well 
known,  a  pavement  system  rut  depth  model  must  be  flexible  so  that 
different  pabing  materials  can  be  in  the  model  with  minimum  inconvenience 
to  the  user.  Thus,  the  rut  depth  model  is  developed  in  modular  form 
so  that  refinements  and  new  developments  related  to  the  permanent 
defomration  characteristics  of  the  various  paving  materials,  as  well 
as  refinements  of  the  primary  response  model  or  other  subsystems 
can  be  readily  adapted  without  haveing  to  re-develop  the  pavement  design 
framework. 

Basically,  a  model  for  predicting  pavement  rutting  consists  of  two 
parts.  The  first  part  is  the  characterization  of  permanent  deformation 
of  pavement  materials  including  materials  to  be  used  for  surface,  base, 
subbase,  and  subgrade,  under  repeated  loading.  The  second  part  is 
the  analytical  pavement  model,  or  primary  response  model,  in  which  the 
stress  distributions  in  the  pavement  system  under  a  vehicle  wheel  load 
can  be  calculated.  These  stresses  together  with  permanent  deformation 
characteristics  of  the  material  are  used  to  estimate  the  permanent 
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Figure  6-3  TRANSVERSE  PROFILES  AASHO  TEST  ROAD 
TRENCH  STUDY  (FROM  REFERENCE  33) 


87 


strains.  Summation  of  the  permanent  strain  with  depth  under  the  wheel 
paths  yields  the  permanent  deformation  of  the  surface.  The  total  rut 
depth  under  a  sequence  of  traffic  loading  (varying  in  magnitude  as  well 
as  in  duration)  is  calculated  by  properly  summing  the  rut-depth  obtained 
from  each  single  load  application. 

The  general  approach  to  the  characterization  of  the  permanent 
deformation  characteristics  of  pavement  materials  under  repeated 
loading  has  been  to  test  the  material  under  conditions  as  nearly 
representative  of  field  conditions  as  possible.  The  test  procedures 
include  the  use  of  the  stresses  (stress  levels  for  axial  and  confining 
stresses)  simulating  those  induced  by  the  wheel  loads  (intensity,  and 
duration),  temperature,  and  moisture  in  the  field.  Using  a  well 
planned  statistical  factorial  design,  the  effects  of  each  parameter 
on  the  permanent  strains  can  be  obtained.  Inasmuch  as  the  relation- 
ships are  derived  for  the  permanent  strains  of  the  paving  materials 
under  simulated  field  conditions,  the  predictions  using  the  equations 
so  derived,  should  be  \/ery   reliable.  A  main  disadvantage  of  this  approach, 
however,  is  that  when  the  number  of  parameters  involved  and  the  levels 
for  each  parameter  are  too  large,  the  experimental  program  becomes 
much  too  cumbersone  and  too  expensive  to  be  of  practical  use. 

Thus,  an  alternate  approach  for  predicting  the  permanent  defor- 
mation characteristics  of  asphalt  concrete,  in  particular,  under 
repetitive  loading  has  been  developed  involving  a  simplified  testing 
procedure  in  which  only  a  relativly  small  number  of  creep  tests  are 
needed  for  the  characterization  of  the  constitutive  equation  for  predicting 
permanent  strains. 
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6.2  GENERAL  FRAMEWORK  OF  THE  RUT  DEPTH  MODEL 

The  rut  depth  model  developed  herein  is  compatible  with  the 
structural  subsystem  in  VESYS  II  or  VESYS  IIM.  As  noted  previously, 
the  rut  depth  model  is  modular  to  provide  for  maximum  flexibility  for 
substituting  various  primary  response  models  and  permanent  deformation 
realtionships  of  the  various  paving  materials.  The  analysis  provedure 
for  calcuating  rut  depth  is  illustrated  schematically  in  Figure  6-4, 
and  a  simplified  flow  diagram  for  the  rut  depth  model  is  shown  in 
Figure  6-5. 

The  geometry  of  the  pavement  system  (e.g.,  layer  thickness),  the 
loading  (considering  a  unit  step  laod  intensity)  and  material  response 
properties  (moduli,  or  compliances)  are  input  to  the  primary  response 
model  which  calculates  stress,  strain  and  deflection  throughout  the 
pavement  structure  for  a  step  load  intensity.  The  stresses  under  the 
wheel  path  thus  calculated  are  the  only  input  from  the  primary  respons 
model  to  the  rut  depth  model.  These  may  be  input  directly  from  the 
primary  response  model  or  read  from  a  file  storage  unit  if  the  two 
subsystems  are  ran  at  different  times.  As  shown  in  Figure  6-5, 
the  vertical  stress  a  and  radial  stress  (confining  stress),  a   ,  at 
several  depths  directly  under  the  wheel  load  are  calculated.  The 
stresses  for  different  load  graoups  (see  Chapter  4)  are  also  calculated. 

The  rut  depth  model  consists  of  several  submodels  for  calculating 
the  permanent  deforamtion  of  each  layer  of  the  pavement  system.  The  rut 
depth  of  the  pavement  system  under  the  given  one  step  loading  (for  a 
given  load  intensity  and  given  load  duration)  is  the  sum  of  the  per- 
manent deformation  accumulated  in  each  layer.  The  total  rut  depth 
under  a  sequence  of  traffic  loading  (varying  in  magnitude  as  well  as  in 
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duration)  is  calculated  by  properly  summing  the  rut  depth  contribution 
from  each  single  load  application.  The  calculation  may  prove  to  be 
quite  lengthy.  A  similar  approach  has  also  been  used  by  Romain  [35], 
McLean  and  Moni smith  [36],  Morris  [37],  and  Barksdale  [38]  for  estimating 
the  rut  depth  of  pabement  systems  due  to  repeated  loading. 

Although  the  logic  of  rut-depth  prediction,  as  mentioned  herein, 
is  rather  straightforward,  the  detailed  calculations  are  not  quite  so 
straight-forward.  First,  the  relationship  of  the  permanent  strains 
to  the  applied  stresses,  temperature,  moisture  levels,  etc.,  for 
most  paving  materials  are  still  relatively  unknown.  Second,  the  problem 
is  complicated  by  the  complex  nature  of  traffic  load  characteristics. 
The  stresses  induced  in  the  pavement  system  due  to  traffic  depend  on 
the  total  weight,  contact  area,  contact  pressure  and  wheel  configuration. 
The  transverse  distribution  of  vehicles  along  the  traffic  lane  will 
affect  the  actual  stresse  induced  and  hence  the  permanent  deformation 
under  the  "wheel  path"*.  Load  duration  also  affects  the  accumulation 
of  permanent  deformation.  Although  under  normal  conditions  only  the 
stresses  induced  from  normal  wheel  loads  are  important  to  the  permanent 
deformation,  the  shear  force  associated  with  acceleration,  braking  and 
cornering  (e.g.,  highway  curves)  could  conceivably  contribute  a  significant 
amount  to  the  rutting.  Third,  the  environmental  effects  of  temperature 
and  moisture  must  be  considered.  Temperature  and  moisture  not  only  affect 
creep  compliances  which  in  turn  affect  the  stress/strain  response  in 
the  primary  response  model,  temperature  and  moisture  also  effect  the 
pavement  deforamtion  characteristics  of  the  layer  materials.  Pavement 


*  The  "wheel  path"  as  used  in  this  report,  is  defined  loosely  as  the 
positions  in  the  transverse  direction  of  the  traffic  lane  in  which 
the  highest  distribution  of  wheel  loads  pass. 
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rutting  due  to  uneven  subgrade  heaving  as  a  result  of  freeze  thaw  action 

is  not  considered. 

6.3  TRAFFIC  LOADING  AND  ENVIRONMENTAL  INPUTS 

The  inputs  for  traffic  loading  and  environment  ( temperature, 
moisture)  for  the  rut  depth  medel  are  the  same  as  for  the  primary 
response  model  and  are  described  in  Chapter  4. 

To  briefly  review,  instead  of  using  a  certain  fixed  intensity 
(such  as  the  equivalent  18  kip  axle  load),  the  traffic  load  input  consists 
of:  (1)  the  actual  wheel  load  distribution,  such  as  the  distribution 
from  the  standard  W-4  form  which  provides  the  number  of  axles  per 
thousand  trucks  (or  vehicles)  in  each  axle  load  group  (10  groups  for 
single  axles,  converting  tandem  axles  into  equivalent  single  axles*); 
(2)  average  daily  traffic,  ADT  (including  mean  and  variance);  (3)  monthly 
variation  of  traffic;  (4)  average  growth  of  traffic  per  year;  (5)  speed 
of  trucks  (mean  and  variation);  (6)  average  radius  of  the  contact  area; 
(7)  channelization  distribution  factor;  (8)  percentage  of  trucks  in 
traffic  stream;  and,  (9)  percentage  of  trucks  in  the  design  lane.  A 
table  which  consists  of  frequency  versus  stresses,  strains  and  de- 
flections per  thousand  vehicles  can  be  obtained  from  input  (1). 
Input  (20),  (3)  and  (4)  are  then  used  to  calculate  the  total  stress, 
strain  and  deflection  in  each  axle  load  group  in  each  month,  each 
year  and  the  total  accumulated  values  for  any  period.  This  infor- 
mation is  then  input  to  the  damage  models  to  calculate  the  pavement 
distresses. 


*  At  the  present  time,  the  tandem-wheel  is  considered  as  two  single 
wneels;  that  is,  the  tandem  axle  load  is  divided  by  two  and  each  is 
considered  as  a  single  axle  wheel  load. 
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The  effects  of  wheel  load  congifurations,  and  shear  force  on  the 
stresses  induced  in  the  pavement  system  are  not  included  in  the  primary 
response  model  of  the  VESYS  II  system.  This  is  due  to  the  assumption 
of  axial  symmetry  for  loading  as  well  as  geometry  of  the  layered  system 
and  automoatically  excludes  multiple  wheel  loads  and  shear  forces. 

The  temperature  is  given  by  the  mean  and  variance  of  air  temper- 
ature for  each  basic  period,  and  moisture  level  is  also  given  by  mean 
and  variance  in  the  subgrade  for  the  same  basic  period.  The  basic 
period  is  usually  taken  to  be  a  month. 

A  heat  transfer  analysis  may  be  used  to  obtain  more  accurate 
estimates  of  the  temperature  in  the  pavement  structure  (e.g.,  Dempsey 
and  Thompson  [39]),;  however,  at  the  present  time  this  type  of  analysis 
is  not  included  for  the  following  reasons:  (1)  inputing  only  mean 
and  variance  of  monthly  temperature  and  moisture  is  not  accurate 
enough  to  warrant  the  use  of  more  accurate  heat  transfer  analyses; 
(2)  only  the  asphalt  concrete  surface  and  asphalt  treated  base  course  are 
strongly  dependent  on  temperature.  The  subgrade  material  is  not  too 
sensitive  to  temperature  level  so  long  as  the  temperature  is  above 
freezing;  hence,  only  a  relatively  rough  estimation  of  subgrade  tem- 
perature with  air  temperature  variation  is  reguired.  Since  temperature 
variations  in  the  surface  and  base  courses  are  relatively  close  to 
the  surface  temperature  variations,  temperatures  in  these  layers  can 
be  estimated  rather  accurately  using  as  approximate  method. 
6.4  MATERIAL  CHARACTERIZATION 

The  material  properties  which  are  pertinent  to  the  analysis  of 
rutting  can  be  catagorized  into  two  groups:  the  material  response 
properties  and  permanent  deformation  properties. 
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Material  response  properties  consist  of  creep  compliances 
(for  viscoelastic  materials)  or  elastic  moduli  (for  elastic  materials) 
and  Poisson's  ratio.  These  properties  of  each  layer  material  are  needed 
in  the  primary  response  modul  for  calculating  the  stresses,  strain 
and  deflections  in  a  pavement  system  under  the  application  of  external 
loadings.  Characterization  of  material  response  properties  of  pavement 
materials  has  been  the  subject  of  extensive  investigation  in  the  past 
ten  years,  (see  for  example,  [24,  40]),  and  therefore  is  not  discussed 
herein.  This  section  concentrates  on  characterization  of  permanent 
deformation  characteristics. 

The  permanent  strain  induced  in  various  paving  materials  under 
an  applied  stress  depends,  in  a  complex  manner,  upon  the  follov/ing 
parameters: 

stress  state,  (axial  and  confining  pressure) 

stress  level 

stress  duration 

number  of  load  repetitions 

environmental  factors  (temperature,  moisture,  etc.) 
Therefore,  a  constitutive  equation,  relating  the  permanent  strains  (e  ) 
to  the  stress  state  and  stress  levels  a..,  load  duration  At,  number  of 
load  repetitions  N,  temperature  T,  and  moisture  content  M,  is 

D 

required  for  each  kind  of  paving  material  to  facilitate  rut-depth 
prediction.  The  general  form  of  this  relationship  is 

£p  =  f\.o.y   At,  N,  T,  M]  (6_1} 

Not  all  materials  necessarily  depend  upon  the  totality  of  these 
parameters,  however, 
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Asphalt  concrete  is  strongly  dependent  upon  temperature 
variation,  but  may  not  be  too  sensitive  to  moisture  level;  whereas, 
the  effect  of  these  factors  upon  clay  soil  could  be  exactly  opposite. 
In  the  same  vein,  if  the  permanent  strains  induced  from  repeated 
loading  are  due  solely  to  "viscous  flow"  such  as  may  occur  in  a  well  compacted 
asphalt  concrete  with  relatively  high  asphalt  content,  the  effect  of 
duration  (At)  under  each  cycle  and  the  total  number  of  load  application 
(N)  can  be  combined  into  a  single  parameter  (At  N).  For  A  material  such  as 
agregates,  which  shows  an  exactly  opposite  behavior,  the  permanent  strain 
due  to  repeated  loading  arises  entirely  from  "time  independent  plastic 
deformation  or  friction",  and  the  permanent  strain  depends  only  on  the 
total  number  of  load  applications  (N)  and  is  independent  of  the  duration 
(At). 

It  is  also  important  to  realize  that,  strictly  speaking,  an 
equation  determined  for  a  given  paving  material  under  a  given  set 
of  conditions  cannot  be  used  to  describe  the  behavior  of  the  same 
material  under  different  conditions.  For  example,  changing  the  asphalt 
content,  compaction,  or  amount  of  air  voids  in  asphalt  concrete  results 
in  a  different  response. 

A  limited  number  of  research  efforts  toward  development  of 
constitutive  equations  for  permanent  strain  for  paving  materials  have 
been  initiated  recently.  At  the  present  time,  the  development  of  the 
constitutive  equation  for  permanent  strains  can  be  categorized  by 
two  approaches: 

(1)  Simulative  statistical  approach, 

(2)  Viscoelastic  -  plastic  approach. 
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In  the  first  approach,  laboratory  samples  of  the  material  are 
tested  under  conditions  as  nearly  representative  of  field  service 
conditions  as  possible.  The  test  procedures  include  use  of  the 
actual  stresses  (actual  stress  levels  for  axial  stress,  confining 
stress)  induced  by  the  wheel  loads,  duration  of  loading,  temperatures, 
and  moisture  levels,  etc.  Based  on  this  preliminary  information,  a 
statistically  designed  test  program  is  undertaken  to  determine  the 
effects  of  the  main  design  parameters  and  their  interactions  on  the 
permanent  strains  using  factorial  analysis  [41].  Using 
multiple  regression  analysis  (linear  or  nonlinear),  an  equation  relating 
the  permanent  strains  to  the  design  factors  can  be  obtained. 

Inasmuch  as  the  relationships  are  derived  for  the  permanent  strains 
of  the  paving  materials  as  a  function  of  the  number  of  load  applications, 
stress  state  and  environmental  conditions  under  simulated  service  conditions, 
the  predictions  of  pavement  rutting  using  these  equations  should  be  \rery 
reliable;  depending  only  upon  the  reproducibility  of  field  service  con- 
ditions in  the  laboratory.  In  principle,  by  using  the  factorial  experiment 
design,  the  effect  of  any  parameter  upon  permanent  strain  can  be  obtained 
as  accurately  as  desired  by  using  small  increments  or  selecting  many  levels 
for  each  factor  in  the  experimental  design. 

The  effort  necessary  for  acceptable  reliability  unfortunately  turns 
out  to  be  the  main  disadvantage  of  this  approach.  If  the  number  of  factors 
involved  and  the  levels  for  each  factor  are   too  large  the  experimental 
program  quickly  becomes  much  too  cumbersome  and  expensive  to  be  of  practical 
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use.  For  example,  for  a  three  factor  design,  with  each  at  three  levels, 
a  factor  design  would  require  3  =  27  tests.  Allowing  for  replications 
to  determine  inherent  material  variabilities,  the  total  number  of  tests 
could  easily  exceed  100.  In  order  to  investigate  permanent  strain  as 
a  function  of  the  number  of  load  applications,  a  minimum  of  10,000 
load  repetitions  per  test  is  required.  Thus,  the  test  proqram  can 
easily  exceed  200-300  hours.  In  addition,  the  apparatus  and  instrumentation 
needed  to  obtain  reliable  results  are  rather  extensive  and  expensive. 
Thus,  the  most  urgent  research  need  in  this  area  is  the  development  of 
simplified  tests  which  decrease  the  total  number  of  tests,  shorten  the 
amount  of  time  required  for  each  test,  and  simplify  the  test  mehtods  and 
instrumentation  requirements. 

The  second  approach  for  development  of  constitutive  equations  for 
permanent  strain  is  based  upon  viscoelastic-plastic  theories.  There  is 
no  universal  theory  that  can  be  used  indiscriminantly  for  all  materials, 
rather,  a  specific  mechanics  theory  can  only  be  used  for  certain 
paving  materials,  which  meet  certain  conditions.  Therefore,  it  is 
anticipated  that  it  may  be  necessary  to  treat  each  material  differently 
with  a  different  theory.  This  fact  makes  this  approach  less  attractive 
and  not  as  easily  understood  as  the  first  approach;  however,  once  a 
particular  theory  is  shown  to  be  applicable  for  a  particular  material, 
the  functional  relationship  between  permanent  strain  and  the  influencing 
factors  can  be  explicitly  developed  from  the  theory.  In  this  approach 
the  size  of  the  test  program  is  reduced  by  an  order  of  magnitude  over 
the  laboratory  based  methods,  and  only  a  relatively  small  number  of 
tests  are  needed  to  determine  the  coefficients  in  the  equation 
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describing  permanemt  deformation.  This  result  makes  this  approach 
very  appealing  from  a  practical  point  of  view.  More  importantly,  the 
characteristics  of  permanent  strains  obtained  in  this  approach  are, 
in  many  cases,  explicitly  obtained  in  terms  of  certain  fundamental 
mechanical  properties.  The  effects  of  physical  characteristics  and 
environmental  factors  on  those  properties  are  relatively  well  known. 
In  contrast,  for  example,  a  change  in  mix  design  of  asphalt  concrete 
requires  a  whole  new  test  program  to  characterize  permanent  strains 
in  the  first  approach. 

Despite  its  obvious  advantages,  the  second  approach  is  not  without 
shortcomings.  For  example,  the  applicability  and  the  accuracy  of  the 
theory  for  predicting  the  permanent  strain  of  a  given  material  is 
undoubtedly  the  central  question  to  be  answered.  As  will  be  seen, 
the  present  state  of  knowledge  in  some  cases,  is  still  not  quite 
sufficient  to  answer  this  question  convincingly.  These  obstacles  will 
certainly  be  overcome  in  the  future  as  more  effort  is  expended  on  this 
approach. 

In  the  following  subsections  the  constitutive  equations  for 
permanent  strain  developed  by  other  researchers  is  briefly  reviewed 
before  presenting  the  specific  viscoelastic-plastic  model  recommended 
for  asphalt  concrete. 
6.4.1  Asphalt  Concrete 

Morris  [37]  used  the  axial  stress  a,,  confining  pressure  o_,  and 

temperature  T  as  the  three  main  design  factors.  Each  has  two  levels  in 

3 
a  2  factorial  design  experiment  for  investigating  the  relationship 

between  the  permanent  strain  of  asphalt  concrete  and  the  number  of  load 
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repetitions.  For  each  niven  set  of  loading  and  temperature  conditions, 
an  ashphalt  concrete  specimen  (4  inches  in  diameter  and  8  inches  long) 
was  subjected  to  a  minimum  of  106  load  repetitions  and  the  permanent 
strains  vs.  number  of  repetitions  were  measured. 


For  the  compression  model  he  obtained 

c 

l09  (tt£n)  =  eo  +  e,  a,  +  e?  a,  +  e,  T  +  e-  o,  t  (6-2) 


log  My    °   H  wl   "2  "3    3  '   "3  w3 
For  the  tension  model, 

where  l     is  the  permanent  strain,  N  is  the  number  of  load  repetitions,  and 
6o  .  .  .  Q  •  B0  .  .  .  B3  are  constants. 

McLean  and  Monismith  [36] used  a  similar  approach  in  investigating  the 
permanent  strain  of  asphalt  concrete  under  repeated  loading,  with  the 
following  results: 

log  ep  *  K[  (a^o3)ce]  ♦  C^log.  M+  C2(log  N)2  +  C3  (log  N)3     (6-3) 

where  e     is  the  axial  permanent  strain,  o,  and  a,  are  the  axial  and  lateral 
stresses,  e  is  the  elastic  strain  and  K  and  a  are  temperature  dependent 
coefficients. 

Soussou,  Moavenzadeh  and  Findakly  [1]  have  suggested  the  following 
equation  to  account  for  the  permanent  strain  of  asphalt  concrete: 

-P  (t)  =  [J  (t)  -  J  (t-At)]  +  (l-vj)J(t-At)e"6e  (6-4) 
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In  this  equation  e  represents  the  permanent  deformation  under  a  qiven 
P  3 

loading  with  a  load  duration  equal  to  At,  J  (t)  is  the  creep  compliance 

of  the  material,  u  is  the  ratio  of  the  creep  compliance  under  unloading 

and  the  creep  compliance  under  loading, and  e    is  introduced  to  account 

for  the  decrease  in  the  accumulation  rate  of  permanent  strain  as  the 

degree  of  compaction  increases.  When  u=l ,  i.e.,  when  the  loading  and 
unloading  behavior  are  the  samo (6-4)becomes 

cp  (t)  =  [J(t)  -  J(t-At)]oo.  (6-5) 

Under  this  situation  the  permanent  strain  is  accounted  for  only  through 
the  rate  sensitivity  of  the  linear  viscoelastic  material.  Although 
linear  viscoelasticity  theory  has  been  used  by  many  researchers  to 
characterize  the  time-dependent  behavior  of  asphalt  concrete,  asphalt 
concretes  are  not,  strictly  speaking,  linear  viscoelastic  materials. 
In  particular,  the  unloading  and  loading  behavior  of  these 
materials  are  quite  different  as  demonstrated  by  Lai  [42,43]. 
Because  of  this,  the  use  of  (5-5)to  describe  the  permanent  strain  in  the 
original  version  of  the  VESYS  II  was  not  ^ery   successful  as  pointed  out 
in  [2,36,37].   Equation  (6-4)should  improve  the  predictive  capability.  In 
(5-4)when  the  stress  pulse  shape  (constant  shape  throughout  or  varying 
shape)  is  qiven,  the  accumulated  permanent  strain  for  any  period  (i.e., 
at  any  number  of  load  repetitions)  can  be  calculated  so  long  as 
J(t),  u,  and  6  for  the  material  are  known.  In  principle,  J(t)  is  the  creep 
compliance  determined  from  the  loading  portion  of  the  creep  test,  and  u 
is  determined  from  both  the  loading  and  unloading  portion  of  the  creep 
curve,  the  constant  6  must  be  determined  from  a  repetitive  test.  It 


01 


is  important  to  recognize  the  fact  that  u,  the  ratio  of  creep  compliances 

under  unloading  J  (t  -  At)5  loading  J  (t  -  At),  as  illustrated  in  Figure  6-6, 

Ju  (t-  At) 
v.  =  J  (t-At)  (6-6) 

depends  upon  the  duration  of  loading  time  (At),  and  even  at  a  given 
duration,  the  ratio  of  these  two  time  dependent  curves  depends  on  time. 

6.4.2  Granular  Base  Material 

The  available  results  on  the  permanent  strain  characteristics  of  granular 

materials  are  relatively  meager.  Barksdale[38]  conducted  rather  extensive 

repeated  load  laboratory  triaxial  tests  on  untreated  base  course  materials. 

The  results  indicated  that  the  effect  of  axial  and'  confining  pressure  on  the 

permanent  strain  can  be  fitted  by  the  following  equation,  which  was  originally 

developed  for  the  static  triaxial  test  for  sand  and  clay  [44]. 

(cyo^)  I   (Ka3n) 

((ara3)Rf        ) 

j  2(c  cos  (J)  +  o3  sin<|>)|  (6~7> 

'    (1  -  sin*)      ) 

where  c  =  axial  plastic  strain  at  100,000  repetitions 
c  =  cohesion 

<J>  =  angle  of  internal  friction 
k,  n,  R,  are  constants 
o-j ,  o0  =  axial  and  confining  pressure  respectively. 

It  was  also  indicated  by  Barksdale  that  the  total  accumulated  plastic 
strain  is  related  linearly  to  the  logarithm  of  the  number  of  load  applications. 
Hence,  the  right  side  of  (6-7)  should  be  multiplied  by  f-  1og,nN.  Thus, 


£P  = 


1 
p    '  5  ,uy10' 


=  F  f  loglnN  (6-8) 
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where  F  is  given  in  (6-7)  and  is  the  actual  accumulated  permanent  strain  at 
100,000  load  applications.  Although  the  factor  F  can  be  estimated  from 
static  tests,  as  noted  in  [44],  a  more  practical  way  of  using  (6-8)  to 
estimate  the  permanent  strain  is  to  determine  F  directly  by  performing 
repetitive  triaxial  tests  (up  to  100,000  cycles  per  test)  for  the  proposed 
base  course  materials  with  the  applied  axial  and  lateral  stresses  chosen 
to  simulate  the  anticipated  actual  stresses  occurring  in  the  pavement  struc- 
ture. With  F  thus  determined,  (6-8)  can  then  be  used  to  estimate  the  permanent 
strain  for  different  load  applications.  If,  however,  the  results  of  triaxial 
confining  repetitive  tests  are  not  available  ,  due  to  lack  of  a  testing 
apparatus  or  other  reasons,  a  static  triaxial  test  may  be  used  to  estimate 
the  F  factor  for  a  given  material. 

So  far,  there  is  no  theoretical  method  available  which  is  capable  of 
predicting  the  accumulative  permanent  strain  of  aggregates  under  repetitive 
loading.  The  work  by  Mandl  and  Lague  [45],  in  which  plasticity  theory  is 
used  to  describe  the  deformation  characteristics  of  aggregates,  may  offer 
promise  for  development  of  a  theoretical  framework  for  prediction  of  the 
permanent  strain  due  to  repetitive  loadings. 

6.4.3  Clay  Soil 

The  results  of  triaxial  load  tests  on  clay  soil  [46,  47]  and  some 
test  results  on  Penn  State  Test  Track  Clay  conducted  by  FHWA,  indicate 
a  linear  permanent  strain  vs.  the  logarithum  of  the  number  of  load 
applications 

z     =  A  log  N  (6-9) 

P 
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The  coefficient  A  depends  on  axial  and  lateral  stresses,  moisture, 
compaction,  and  temperature  for  a  given  type  of  clay.  At  the  present 
time,  it  seems  that  the  most  practical  way  of  estimating  the  permanent 
strain  under  repeated  loading  is  to  use  (6-9)  with  the  coefficient 
A  determined  in  the  laboratory  under  tesing  conditions,  such  as  axial 
and  lateral  stresses,  moisture,  simulating  the  conditions  anticipated 
in  the  field. 

It  has  been  pointed  out  [48,  49]  that  clay  soils  exhibit  creep 
phenomenon  when  subjected  to  a  constant  stress  and  exhibit  partial 
recovery  and  partial  permanent  strain  upon  stress  removal.  This 
phenomenon  seems  to  indicate  tnat  the  theoretical  approaches  discussed 
subsequently  can  be  applied  to  develop  a  simplified  method  for  estimating 
the  accumulated  permanent  strain  of  clay  soils  from  creep  tests.  Figure 
6-7  shows  the  results  of  permanent  strain  induced  from  different  load 
durations  versus  the  laoding  duration  of  a  series  of  creep  tests  of 
Penn  State  Test  Track  Clay.  A  power  law  relation  between  the  permanent 
strains  and  the  duration  of  loading  time  for  this  material  is  clearly 
indicated. 
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For  cement-treated  base,  it  is  assumed  in  VESYS  I IM  that  this 
material  does  not  exhibit  permanent  deformation. 
6.4.4  Spatial  Variations  of  Materials  Properties 

This  coefficient  is  a  measure  of  tne  statistical  correlation  of 
two  points  at  a  given  distance  of  each  other.  This  coefficient  p  may 
be  represented  by  the  following  equation 

n 
c 

p  =  A  +  B  exp  (^-) 
c 

A  +  B  =  1 

where  x  represents  the  distance  between  two  points  to  be  correlated. 

This  equation  implies  that  two  neighboring  locations  (x  smaller)  are 

moer  likely  to  have  similar  properties  (larger  p)  than  two  locations 

which  are  farther  apart. 

This  coefficient  of  spartial  variations  is  used  to  facilitate  the 

roughness  prediction. 

6.5  DEVELOPMENT  OF  A  MODEL  FOR  PREDICTING  PERMANENT  DEFORMATION 
OF  ASPHALT  CONCRETE 

6.5.1  Model  Development 

Asphalt  concrete,  as  noted  previously  and  shown  schematically 
in  Figure  6-6,  and  also  in  Figure  6-3,  exhibits  creep  behavior  under 
a  sustained  load.  When  the  applied  stress  is  small  (about  30  psi 
unconfined,  at  room  temperature),  the  creep  strain  is  approximately 
proportional  to  the  applied  stres.  At  higher  stress,  the  creep  strain 
become  nonlinearly  dependent  on  the  applies  stress  [42]. 

Thus,  the  creep  compliance  J(t,a)  of  that  material,  which  is  defined 
as  the  time  dependent  creep  strain  divided  by  the  applied  constant  stress, 
can  be  characterized  from  Figure  6-8. 
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J(t,a)  =  ^  (6-10) 

Creep  compliance  could  be  stress  dependent  if  the  response  is  non- 
linear. For  example,  from  Figure  6-9,  the  following  results  can  be 
obtained 

J(t,a)  =  A(a)  tn  (6_n) 

A(a)  =  Aq  +  A^  +  A2a2  (6-12) 

The  creep  behavior  of  asphalt  concrete  becomes  more  complicated 
when  the  material  is  subjected  to  a  more  complex  loading  history. 
For  example,  the  recovery  behavior  of  the  material  as  shown  in  Figure 
6-9  indicates  that  the  material  exhibits  permanent  strain  which  depends 
upon  the  duration  as  well  as  the  intensity  of  the  loading.  The 
recovery  curves  as  shown  in  Figure  6-9  cannot  be  predicted  with 
reasonable  accuracy  by  using  existing  viscoelastic  theory.  For  example, 
using  the  superposition  principle,  or  modified  superposition  principle 
[50]  the  recovery  strains,  e  (t),  can  be  predicted  from  creep  compliance 
as  follows 

ef(t)  «  o0[J(a0,t)  -  J(o0,  t-tj)]  (6-13) 

t  >  t1 

where  a  1s  the  applied  constant  stress  and  t,  1s  the  duration  of  the 
o  i 

loading.  Results  of  the  predictions  using  (6-13)  are  shown  in 
Figure  6-9.  Similarly,  using  the  modified  superposition  principle 
which  has  been  used  quite  successfully  on  many  other  materials  [51, 
52]  yeilds  unaccpetable  prediction  of  the  creep  behavior  of  asphalt 
concrete  under  multiple  steps  loading  as  shown  in  Figure  6-10. 
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The  large  portion  of  unrecoverable  strain  of  each  creep  test 
as  shown  in  Fiture  6-9  and  the  observation  that  the  amount  of  unrecov- 
erable strain  is  dependent  on  the  length  of  the  loading  period, 
suggests  that  the  creep  behavior  may  be  represented  by  a  generalized 
nonlinear  Kelvin  model  as  shown  in  Figure  6-11.  The  nonlinear  dash- 
pot  contributes  to  the  irrecoverable  strain,  while  the  series  of  Kelvin 
models  (Kelvin  chain)  contribute  to  the  power-law  creep  behavior  as 
discussed  in  [51].  Therefore,  the  total  creep  strains  (e-r)  are  separated 
into  two  parts,  the  irrecoverable  strains  (e  )  due  to  the  nonlinear 
dashpot  and  the  completely  recoverable  strains  (e  )  due  to  the  non- 
linear Kelvin  chain  (see  Figure  6-12). 

These  two  strain  components  can  be  characterized  from  a  series 
of  creep  tests.  Each  creep  test  has  different  stress  level  and  duration 
of  loading.  During  each  test,  both  the  loading  creep  curve  and  re- 
covery curve  including  the  permanent  strain  have  to  be  measured.  The 
permanent  strain  from  each  test  is  plotted  against  the  duration  of 
loaidng  time  can  be  obtained.  For  example,  (6-14)  was  obtained  from 
the  test  results  discussed  in  a  later  part  of  this  section. 

ep(t)  =  A(a)  ta  (6-14) 

This  equation  can  also  be  interpreted  as  the  amount  of  irrecoverable 
strain  induced  at  any  instant  of  time  t  during  a  creep  test. 

Once  this  is  done,  the  recoverable  strain  in  each  creen  test  can 
be  estimated  by  subtracting  the  e  calculated  from  (6-14)  from  the  total 


me  a 


sured  creep  strain  e-r(t).  Thus 

ev(t)  =  eT(t)  -  e  (t)  (6-15) 
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Again,  using  the  same  technique  of  cross-plotting  e  (t)  versus  t, 
an  analytical  form  for  ev(t)  can  be  obtained.  Equation  (6-16)  was 
thus  obtained 

ey(t)  =  B(a)  t6  (6-16) 

Thus  under  a  constant  stress  creep  and  recovery  test,  the  creep  strain 
eT  and  recovery  strain  c  can  be  predicted  as  follows: 

eT(t)  =  B(a)  t6  +  A(a)  ta  ,  t  <  t  (6-17) 

er(t)  =  B(a)  [t6  -  (t  -  t,)3]  +  A(o)  t^  ,  t  >  t  (6-18) 

Although  the  coefficients  A(a),  B(a),  a,  and  6  in  (6-14)  and 
(6-16),  which  may  depend  upon  test  temperature, were  determined  separately 
as  mentioned  above,  it  seems  that  the  nonlinear  optimization  techniques 
as  proposed  by  Distefano  [53,  54]  and  Hufferd  [55]  can  be  used  to 
evaluate  the  coefficients  from  the  complete  creep  and  recovery  test 
results. 

It  has  been  shown  in  [42]  that  under  time-dependent  stress  input, 
the  recoverable  strain  can  be  predicted  by  using  the  superposition 
principle  (for  linear  behavior)  or  the  modified  superpostion  principle 
(for  nonlinear  behavior)  and  the  irrecoverable  strains  can  be  predicted 


more  accurately  by  using  a  strain-hardening  theory.     Thus 

t 


CT 


j    (t_v  uMo  «)«+  /  yoc* 


(6-19) 


o 

where  the  premanent  strain  rate  can  be  determined  from  (6-14) 

e  -  a  A  t 
P 
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By  eliminating  t  from  the  previous  equation,  e  kpcomes 


ct-1 

=  a  (  f.  ) 
P 


a      a 
i     =  a(e)     (A)  (6-20) 


Equation  (6-19)  and  (6-20)  were  used  in  [42]  to  predict  the  creep  of 
asphalt  concrete  under  the  multi-step  laoding  as  shown  here  in  Figure 
6-10.  The  results  compare  much  more  closer  with  the  test  results  than 
the  predictions  using  strictly  viscoelasticity  theory. 
6.5.2  Behavior  of  Asphalt  Concrete  Under  Repeated  Loading 

Equation  (6-19)  and  (6-20)  can  be  used  to  predict  the  creep 
behavior  of  asphalt  concrete  under  repeated  loading.  In  general,  the 
multiple  steps  stress  inputs  as  illustrated  in  Figure  6-1 3A  can  be 
expressed  in  a  single  algebraic  equation  as  follows: 


N 


(t)  -  ^T  (oi  -   0l-1)  H(t  -  t±)  (6-21) 


i-0 


where  H(t  -  t.)  is  the  Heaviside's  unit  function  which  has  the  value  1 

when  t  >  t.,  o  when  t  <  t  ,  a-1  =  0  and  t  =  0.  For  repeated  loadinq 
1  o  o 

with  variable  stress  amplitude  but  constant  duration  of  loading  (aT) 

and  constant  recovery  period  (l-a)T  in  each  load  application  as  illustrated 

in  Figure  6-13B  (6-21)  becomes 

N-l 
a(t)  =  a2N  H  (t-NT)  +  ^T   o^  [  H(t-iT)  -  H(t-it-aT)]      (6-22) 

i»0 
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for  NT  < t < NT  +  aT  (loading  portion  at  n-th  cycle),  and 
N 


(t)  -  ^2  a2±  [H(t-iT)  *  h<*-it-»t>j 


(6-23) 


i»0 


for  NT+aT  <  t  <  (N+1)T  (recovery  portion  at  n-th  cycle). 

Inserting  (6-22)  and  (6-23)  into  (6-19)  and  (6-20),  the  recoverable 

and  irrecoverable  strains  at  any  instance  of  time  can  be  calculated. 

For  example,  the  strains  at  the  end  of  the  loading  portion  in  the  n-th 

cycle  (t= NT+aT)  are 

N-l 

ev  =  B(o2,T)(aT)6  +  (T)8  ^Z  B(°2i)  |>+a-i)B-<N-i)B  ]  ,    (6-24) 


(aT) 


i«0 


N 


22  A(°2i) 


1/a 


i-0 


(6-25) 


The  strains  at  the  end  of  the  recovery  period  in  the  n-th  cycle, 
t  =  (n+l)T,  are 


N 


s  - (T) 


22  B(o2i)  [ 


(N+l-i)3  -  (N+l-i-a)6J     (f>-26) 


and  c  si  the  same  as  given  in  (6-25).  Some  useful  results  can  be 
p 

obtained  from  (6-24),  (6-25)  and  (6-26)  which  are  discussed  briefly 
in  the  following  sections. 
6.5.3  Effect  of  Load  Duration 

From  (6-25),  it  can  be  seen  that  the  permanent  strain,  z   ,  under 
repeated  loading  is  dependent  on  the  duration  of  loading  time  (aT). 
If,  the  stress  amplitude  is  the  same  throughout  the  test  peiod,  (6-25) 
reduces  to 
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Ep  =  A(a)(aTN)a  =  A(a)(aT)a  (N)a  (6-27) 

This  power  law  relationship  has  been  found  to  be  suitable  in  represent- 

ing  the  accumulative  permanent  strain  under  repeated  stress  cycles 

from  the  test  results  shown  in  the  next  section,  and  has  also  been 

reported  by  Brown  and  Snaith  [56].  In  [56]  it  was  also  confirmed  that 

the  total  permanent  strain  was  dependent  on  the  total  effective  duration 

of  loading  (aTN)  instead  of  the  total  number  of  repetitions  (N). 

It  can  also  be  shown  that  for  a  less  than  one  (true  in  all  cases), 

the  actual  permanent  strain  during  each  loading  period  decreases. 

That  is,  the  accumulation  rate  decreases  as  the  number  of  load  applications 

increases. 

The  cumulative  strain  contributed  from  the  recovery  strains  e  ,  as 

given  in  (6-24)  and  (6-25)  are  relatively  small.  For  example,  under 

constant  stress  amplitude,  (6-26)  becomes 

N 

%-*(o)(T)3    21  <N+1-i)6 "-  (NH-I-a)6  (6"28) 

i=0 

For  0  <  3  £  0.2  which  is  valid  for  most  asphalt  materials  and  for 
0  <  a  <  <  0.5,  the  contribution  of  (6-28)  compared  with  (6-27)  fov 
a  large  number  of  load  repetitions  becomes  very   small  (about  10%).  The 
calculations  of  (6-28)  for  large  N  (such  as  10  )  becomes  very   involved. 
However,  by  assuming  that  a  and  3  in  (6-27)  and  (6-28)  are  the  same, 
the  ratio  of  e  and  e     is  approximately  equal  to 


ev    B(a 


ep  "  A(a 


(crB)  (6-29) 
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6.5.4  Effect  of  Rest  Period 

Equations  (6-27)  or  (6-25)  indicate  that  e  is  independent  on 
the  rest  period,  while  e  as  given  in  (6-26)  or  (6-28)  shows  that,  for 
a  constant  T,  the  longer  the  rest  period  (smaller  a),  the  smaller  the 
recoverable  strain.  For  smaller  a,  or  equivalently,  longer  resting 
period,  (1-a),  the  contribution  of  e  to  the  total  cumulative  strain 
is  small,  particularly  when  6  is  small.  The  fact  that  the  rest  period 
has  little  effect  on  the  accumulation  of  the  permanent  strain  under 
cyclic  loading  has  also  been  reported  in  [56]. 

6.5.5  Effect  of  Mixed  Stress  Amplitude 

Since  the  contribution  of  e  to  the  total  permanent  strain  is 
small,  the  discussion  on  the  effect  of  mixed  stress  amplitude  is  concen- 
trated on  the  e  part.  Equation  (6-25)  indicates  that  under  constant 
duration  of  loading  in  each  cycle,  the  total  accumulated  strain  is 
independent  of  the  dequence  in  which  the  different  stress  amplitudes 
are  applied,  so  long  as  the  combination  of  the  various  stress  amplitudes 
in  the  entire  period  is  the  same. 

Equation  (6-25)  and  (6-27)  can  also  be  used  to  determine  an 
equivalent  stress  amplitude  which  will  yield  an  equivalent  total  permanent 
strain  at  the  same  total  number  of  repetitive  loadings.  These  equations 
can  also  be  used  to  determine  different  combinations  of  stress  levels 
versus  number  of  repetitions  which  will  all  yield  the  same  total 
permanent  strain  under  a  sequence  of  mixed  stress  levels.  This  can  be 
done  so  long  as  the  constant  a  and  the  coefficient  A  as  a  function  of 
stress  have  been  determined  from  a  series  of  creep  tests. 
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6.6  COMPARISON  OF  PERMANENT  DEFORMATION  MODEL  WITH  EXPERIMENTAL  RESULTS 

In  order  to  test  the  accuracy  of  (6-24),  (6-25)  and  (6-26)  in 
predicting  the  accumulative  permanent  strain  under  repetitive  loading, 
a  series  of  creep  tests  and  repetitive  loading  tests  v/ere  conducted. 

6.6.1  Materials  and  Specimens* 

A  single  gradation  of  aggregates  as  shown  in  Table  6-1  were  used. 

Table  6-1  Gradation  of  Aggregates 
Sieve  Size    3/8"    #4    #8    #16    #50    #200 
%   Passing     100    95.0   60,0    33.0    15.0    6.0 

The  asphalt  cement  used  in  preparing  test  specimens  was  from  Phillips 
Petroleum  Company  with  85-100  penetration  grade.  The  asphalt  content 
was  8%  by  weight. 

The  test  specimens  used  in  this  study  were  compressed  asphalt  con- 
crete cylinders  2"  in  diameter  and  3"  in  length.  The  asphalt  mixture  was 
pressed  in  a  2"  diameter  mold  at  300  psi  for  5  minutes  at  a  temperature 
of  200°F.  The  specimens  were  cured  in  an  oven  at  120°F  for  72  hours 
prior  to  the  testing.  Also,  in  order  to  minimize  the  randomness  of  the 
test  results,  each  specimen  was  subjected  to  a  50  psi  prestress  for 
12  minutes. 

6.6.2  Creep  Tests 

A  series  of  creep  tests  were  conducted  with  30  psi  unconfined 
compressive  stress.  Each  specimen  was  tested  at  different  duration  of 
loading  time  (2,  10,  100,  1000  seconds).  At  the  end  of  the  loading 


*  The  materials  and  specimens  used  are  the  same  as  that  used  in  [42] 
except  that  9%  asphalt  content  was  used  in  [42]. 
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period,  the  specimens  were  allowed  to  recover  until  the  recovery  curve 
leveled  off.  Figure  6-14  shows  the  results  of  permanent  strain  induced 
for  different  loading  durations  versus  the  loading  duration.  These 
data  were  fitted  by  a  power  law  of  the  form  for  both  the  irrecoverable 
and  recoverable  strain  components 


e  =   A  (t  )a  ,   =0.4 
p    pv  o' 


e  =  B(t)p   ,   =  0.115 


6.6.3  Repetitive  Tests 

A  series  of  repetitive  tests  were  conducted  at  frequencies  between 
0.2  and  5H  .  The  first  set  of  the  repetitive  tests  had  equal  loading 
and  unloading  time,  though  each  test  had  different  frequency.  The 
accumulative  permanent  strain  versus  number  of  load  repetitions  is  shown 
in  Figure  6-15.  It  is  shown  clearly  in  this  figure  that  the  accumulative 
permanent  strain  depends  on  the  frequency.  However,  replotting  the 
results  of  accumulative  permanent  strains  versus  effective  loading  time 
as  shown  in  Figure  6-16  indicates  that  the  effective  load  duration 
time  is  perhaps  a  better  measure  of  the  accumulative  permanent  strains 
in  repetitive  loading  tests.  Superimposing  the  equation  determined 
from  creep  tests  as  shown  in  Figure  6-14  on  Figure  6-16  indicates 
that  (6-27)  which  can  be  determined  from  a  series  of  creep  and  recovery 
tests,  can  also  be  used  to  predict  the  accumulative  permanent  strain 
under  repetitive  loading  with  resonable  accuracy. 

The  second  series  of  repetitive  tests  were  conducted  in  which  the 
loading  times  were  kept  constant  and  the  recovery  time  were  varied. 
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The  results  shown  in  Figure  6-17  indicate  that  the  longer  the  resting 
time,  the  smaller  the  accumulative  permanent  strain. 

In  the  third  series  of  repetitive  tests,  each  repetitive  test 
consisted  of  mixing  stress  amplitude.  The  results  shown  in  Figure  6-18, 
in  which  three  stress  levels  of  20  psi ,  30  psi  and  40  psi  were  applied 
in  different  sequence,  indicate  that  at  the  end  of  54  cycles  the  total 
acumulation  of  permanent  strain  is  the  same.  Additional  tests  at 
higher  frequencies  and  for  a  larger  number  of  cycles  are  needed, 
however,  to  further  verify  this  result. 
6.7  PROBABILISTIC  PERMANENT  DEF0RAMTI0N  MODEL 

One  of  the  primary  objectives  in  the  development  of  the  rationally 
based  pavement  design  framework  is  to  take  into  account  the  variability 
of  the  material  properties,  traffic  loading  and  environmental  factors. 
In  the  primary  response  model,  as  mentioned  before,  the  input  material 
properties  (creep  compliance,  resilience  modulus,  etc.)  are  considered 
probabilistic  in  nature  and  the  output  stresses,  strains,  and  deflections 
are  also  probability  functions,  even  when  the  input  traffic  loading 
adn  environmental  factors  are  deterministic.  The  variability  of  the 
input  from  the  primary  response  model  is  further  compounded  by  the 
variability  of  temperature  and  moisture. 

Thus,  the  output  from  the  rut-depth  model  is  expected  to  be 
stchastic  inasmuch  as  the  output  from  the  primary  response  model  is 
input  to  the  rut-depth  model  plus  the  fact  that  the  permanent  defor- 
mation characteristics  of  paving  materials  under  repetitive  loadings 
are  expected  to  be  variable. 

In  (6-2),  (6-3),  (6-8)  and  (6-9),  the  variability  of  permanent 
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strain  under  a  given  set  of  conditions  (given  stresses,  temperature, 

etc.)  can  be  obtained  directly  from  a  number  of  repeated  tests  under 

the  given  conditions.  This  is  shown  schematically  in  Figure  6-19. 

In  principle,  if  a  sufficient  number  of  repetitive  tests  can  be  conducted, 

a  constitutive  equation  for  permanent  strains  whcih  take  into  account 

the  variability  as  a  function  of  stresses,  etc.,  can  be  developed,  for 

example 

e  =  A^,  Oy    ...)  log  N  +  DA  (a] ,  Oy   ...)  log  N       (6-30) 

In  (6-30)  Dn  is  the  standard  deviation  of  A.  Both  A  and  D,  are  functions 
of  the  testing  conditions  (a, ,  Oy   etc.).  Since  (6-30)  is  determined  in 
the  laboratory,  a,,  cu,  etc.,  in  (6-30)  are  treated  as  deterministic 
variables. 

One  of  the  disadvantages  of  (6-30)  is  that  even  though  the 
probabilistic  nature  of  the  stresses  obtianed  from  the  primary  response 
model  are  directly  influenced  by  material  variability,  it  is  not  possible 
to  separate  the  influence  of  material  variability  upon  permanent  strain 
from  that  of  the  stresses. 

In  the  primary  response  model  in  VESYS  IIM  the  variability  of 
stress  is  expressed  in  terms  of  its  mean  and  deviation  (a,  D  ).  The 
variability  of  permanent  strain  can  be  estimated  in  terms  of  mean  stress 
and  deviation  as  follows: 


e  =  A(a, ,  ay    ...)  log  N  +  D.(a, ,  a~,  ...)  log  N        (6-31) 
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In  (6-4)  and  (6-27)  however,  the  variabil ity  of  permanent  strain 
can  be  obtained  indirectly  from  creep  tests.  The  permanent  strain 
including  the  variability  of  (6-4)  can  be  rewritten  as  follows: 

e  =  F(t)  a  +  DF  a  (6-33) 

where 

F(t)  =  [J(t)  -  J(t-At]  +  (1-u)  J(t-At)  e"6e             (6-34) 

DF(t)  =  (1-u)  Dj(t-At)  e"3£  (6-35) 

In  (6-33)  and  (6-35),  DF  is  the  deviation  of  the  permanent  strain,  and 
D,  in  (6-34)  is  the  deviation  of  the  creep  compliance.  The  mean  and 
deviation  of  permanent  strain  for  a  variable  stress  is  given  by 


ep  =  F(t)  a  +  DF(t)  a  (6-36) 

and 

/   \2  2    2 

[D£  J  =  [F(t)  =  DF(t)]   (D0)  (6-37) 

In  (6-27)  the  variability  of  the  permanent  strain  is  obtained 
from  the  variability  of  the  irrecoverable  portion  of  the  creep  tests. 

e     =  A  (a)  [N(aT)]a  +  Dft  (a)  [N(aT)f  (6-38) 


In  (6-38)  D„(a)  is  the  standard  deviation  of  the  irrecoverable  part  of 
the  creep  tests.  For  example,  the  mean  values  and  the  standard  deviations 
of  the  irrecoverable  part  of  the  creep  tests  of  the  asphalt  concrete 
cores  taken  from  the  Green  River  Test  Road  [57]  indicate  that  both  of 
these  properties  depend  on  the  mix  design  and  the  age  of  the  materials. 
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Higher  coefficient  of  variations  were  observed  for  low  asphalt  content 
mixtures  and  for  older  asphalt  concrete  cores.  When  the  stress  varies, 
the  permanent  strains  of  (6-38)  become 


ep  =  Ap(a)  [N(aT)]a  +  Dz(5)  [N(aT)f 
(\YM    [(^)2  +  (^)l  (Da>2  WCT)]"      (6-39) 


The  variability  of  the  loading  duration  time  aT  under  each  application 
can  also  be  incorporated  in  (6-39). 

It  seems  that  (6-38)  and  (6-30)  offer  a  more  practical  applicatio 
for  those  materials  which  exhibit  flow  (plastic  deformation)  under 
constant  loading.  Not  only  are  no  additional  tests  needed,  aside  from 
those  required  for  the  primary  response  model,  but  more  importantly, 
the  variability  of  the  permanent  deformation  is  tied  in  directly  with 
the  variability  of  the  materils  properties  used  in  the  primary  response 
model  in  which  the  stresses  and  strains  in  the  pavement  system  are  to 
be  evaluated. 

At  the  present  time,  the  major  obstacle  in  using  these  approaches 
is  the  lack  of  sufficient  experimental  results  to  verify  these  theories. 
For  that  matter,  more  experimental  results  are  needed  ont  only  to  verify 
(6-38)  and  (6-39),  but  also  to  simplify  them  and  to  simplify  (6-31)  and 
(6-32)  as  well.  Eventually,  when  sufficient  experimental  data  has  been 
obtained,  it  may  be  possible  to  decide  on  which  approach  should  be 
taken  and  develop  a  concise  material  characterization  program.  At  the 
present  time,  it  is  best  to  allow  users  the  flexibility  of  selecting 
the  type  of  theory  to  be  used  for  different  materials.  This  selection 
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be  based  upon  user's  perference  and  the  availability  of  resources  in 

conducting  the  necessary  testing. 

6.8  MODULAR  FRAMEWORK  FOR  RUT-DEPTH  PREDICTIONS 

As  mentioned  earlier,  the  VESYS  II  Computer  program  and  design 
framework  have  not  yet  reached  the  final  stage  of  development.  There 
is  still  room  for  modifications  and  improvements  to  the  program.  In 
the  primary  response  model,  different  versions  of  layer  systems  and 
finite  element  methods  can  be  used  to  improve  the  predictive  capability 
of  the  model.  On  the  other  had,  in  the  rut-depth  submodel,  there  is 
still  need  for  a  rather  large  effort  directed  at  development  of  constitutive 
equations  for  permanent  strain  for  various  materials.  Similarly,  other 
damage  submodels  require  further  modifications  and  improvement  as  well. 
Because  of  this,  the  VESYS  IIM  is  developed  in  modular  form  so  that 
any  improvement  or  modifications  on  any  part  of  the  system  can  be 
readily  incorporated  into  the  system  without  having  to  change  the  rest 
of  the  system. 

The  modular  concept  for  the  rut-depth  submodel  is  shown  schematically 
in  Fiture  6-20.  In  this  submodel,  the  input  data  is  the  stresses  (means 
and  deviations)  at  various  prsitions  in  the  pavement  sustem.  The  stresses 
are  the  output  from  the  VESYS  IIM  primary  response  model  or  any  other 
source.  The  calculations  of  the  total  permanent  deformation  within  each 
layer  is  also  a  modular  form;  that  is,  the  modules  for  each  can  be  readily 
replaced  or  omitted.  For  example,  for  a  full  depth  asphalt  pavement, 
the  modules  for  calculating  the  permanent  deforamtion,  other  than  the 
ones  for  asphalt  concrete  and  subgrade,  can  be  omitted.  Within  each 
module  for  calculation  of  the  permanent  deforamtion  of  each  layer,  the 
constitutive  equation  for  permanent  strain  is  also  modular  in  that 
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Figure  6-20    Modular  Form  of  the  Rut-Depth  Submodel 

(Note:  Subroutines  enclosed  by  dashed  lines  are  independent  modules  which  can 
be  replaced  or  in  some  cases  omitted  without  affecting  the  overall 
calculations. ) 


*See  Table  6-2 
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different  constitutive  equations  can  be  substituted  depending  upon  the 
user's  preference  with  minimal  incovenience  to  the  user.  For  example, 
in  the  module  for  calculatio  of  permanent  deformation  of  an  asphalt 
concrete  surface  course,  different  constitutive  equations  such  as  AC-1, 
AC-2,  AC-3,  or  AC-4  (see  Table  6-2)  may  be  used. 
6.9  CONCLUDING  REMARKS 

The  philosophy  and  methodology  for  the  development  of  a  realistic 
model  for  rut-depth  predictions  has  been  presented. 

The  primary  motivation  for  the  procedure  adopted  is  the  need  to 
provide  for  maximum  flexibility  in  the  particular  permanent  deformation 
characteristics  observed  of  different  paving  materials  and  the  different 
layers  of  the  pavement  system,  and  the  establishment  and  verification 
of  simplified  experimental  methods  for  rut-depth  predictions  which  offer 
substantial  appeal  from  a  practical  highway  designer's  point  of  view. 

A  simplified  method  for  predicting  permanent  deformation  of  asphalt 
concrete  under  repeated  loading  using  the  results  of  creep  and  recovery 
tests  is  proposed.  The  prediction  of  permanent  deformation  using  this 
method  compared  with  the  experimental  results  of  asphalt  concrete 
under  repeated  loading  is  excellent.  Further  verification  and  implemen- 
tation of  this  simplified  method  and  the  overall  rut=depth  model 
proposed  herein  is  required,  however. 

The  approach  is  based  on  the  premise  that  the  permanent 
strain  under  repeated  loadinq  is  due  solely  to  "viscous  flow," 
which  for  a  well  compacted  rich  mix  asphalt  concrete  should  be  quite 
accurate.  The  permanent  strain  of  asphalt  concrete  inducted 
during  creep  tests  at  different  loading  periods  is  used  to  predict 
the  permanent  strains  induced  under  repetitive  loadings.  This  is 
accomplished  by  conducting  several  creep  tests  at  several  stress  levels 
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for  different  loading  durations.  During  each  test,  both  the  loading  and 
recovery  portions  of  the  creep  test  are  recorded.  A  creep  compliance 
is  determined  from  the  loaidng  portion  of  the  test.  From  the  un- 
loading portion  of  the  results  for  different  durations  of  loading  time 
an  equation  relating  the  amount  of  total  permanent  strain  to  the  loading 
time  duration  can  be  established.  The  effect  of  load  duration  per 
cycle  and  the  number  of  load  applications  on  the  permanent  deformation 
can  be  combined  into  a  single  "effective"  load  duration  parameter. 
The  approach  discussed  is  more  applicable  to  a  rich  mix  asphalt  concrete, 
for  the  reasons  meantioned  before.  For  asphalt  concrete  with  low 
asphalt  content,  the  total  permanent  strain  due  to  repeated  loading 
can  no  longer  be  considered  as  solely  due  to  "viscous  flow"  and  hence 
(6-24)  -  (6-26)  alone  are  no  longer  adequate  to  describe  the  permanent 
deformation.  It  is  more  likely  that  the  following  equation  should  yield 
a  better  prediction 

e  =  A(AtN)a  +  (N)?  (6-40) 

The  first  term  in  this  equation  is  the  same  at  (6-27)  and  the 
second  term  is  considered  due  to  "time  indipendent  plastic  deformation" 
or  "friction"  as  against  "viscous  flow".  In  principle,  from  the  creep 
and  recovery  test  results,  and  plotting  irrecoverable  strain  versus  the 
load  duration,  similar  to  Figure  6-14,  the  following  representation  can 
be  obtained 


ep  =  A<V  +  y  (6_41) 


The  exponent  for  the  second  term  of  (6-41)  can  be  determined  from  a 
single  additional  repeated  loading  test. 
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CHAPTER  7 

CONCEPTS  FOR  DESIGN  OPTIMIZAITON  OF 
FLEXIBLE  PAVEMENT  SYSTEMS* 

7.1   INTRODUCTION 

The  primary  objective  of  this  report  is  the  development  of  concepts 
for  optimization  of  design  decisions  for  flexible  pavement  systems  based 
on  viscoelastic  structural  models.   The  research  approach  adopted  was 
to  review  previous  work;  study  conceptual  interrelations  of  system  compo- 
nents, specific  decision  criteria  and  optimizing  techniques;  define  and 
recommend  phasing  for  the  future  developments  needed  and  to  propose  a 
conceptual  program  for  the  first  phase.   The  results  of  these  studies 
are  discussed  in  detail  in  subsequent  sections. 

There  has  been  a  very  considerable  amount  of  work  expended  toward 
development  of  a  rational  flexible  pavement  design  and  management  system 
that  would  allow  optimal  design  decisions  based  on  the  predicted  lifetime 
performance  and  economics  of  flexible  pavement  systems.   Such  systems 
exist  and  are  in  use  [  58  -  60  ] .   However,  knowledgeable  persons  under- 
stand that  this  is  a  beginning  and  that  years  of  continued  use,  feedback, 
and  refinement  are  required  to  complete  this  effort.   The  potential 
interim  and  long-term  payback  on  this  investment  is  quite  large  when  one 
considers  the  billions  of  dollars  spent  for  pavement  maintenance  and 
the  huge  savings  in  time  and  money  that  may  be  realized  through  very  small 
improvements  in  our  ability  to  optimize  design  decisions. 

The  successful  development  of  the  VESYS  System  into  a  complete 
Pavement  Analysis  and  Design  System  (PADS)  is  considered  to  be  heavily 
dependent  on  expediting  realization  of  its  capability  as  a  complete 
system.   Recommendations  have  been  made  in  a  subsequent  section  for  accom- 
plishing this  as  the  next  phase  of  research  and  other  recommendations 
are  included  for  further  development  in  subsequent  phases. 


'•'  "rhis  chapter  was  prepared  under  subcontract  to  Austin  Research  Engineers  (ARE) 
by  J.B.  Ranhut  and  W.R.  Hudson.   The  ideas  presented  represent  those  of  ARE.  Important 
contributions  to  the  ideas  presented  were  made  by  Drs.  B.F.  McCullough  and  M.  Darter 
and  Mr.  Harvey  J.  Treybig,  to  whom  the  authors  express  their  appreciation. 
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7.2  REVIEW  AMD  SYNTHESIS  OF  PREVIOUS  WORK 

Conceptual  work  on  the  development  of  a  systematic  approach  to 
pavement  design  was  carried  out  by  Hudson,  Finn  and  McCullough  in  1967 
to  1970  [  60  ]  with  the  support  of  the  National  Cooperative  Highway 
Research  Program  (NCHRP),  Highway  Research  Board.   This  work  modified 
and  extended  previous  work  by  Scrivner,  McFarland  and  Carey  [  61  ] 
that  had  in  1968  culminated  in  a  computerized  flexible  pavement  design 
method.   An  updated  and  improved  system  called  Systems  Analysis  Method 
for  Pavements  (SAMP)  with  a  computer  program  designated  in  [60]  as  SAMP5 
was  one  product  of  the  NCHRP  study  by  Hudson,  Finn  and  McCullough. 

Project  123  was  initiated  in  1969  as  a  joint  effort  by  the  Center 
for  Highway  Research,  University  of  Texas  at  Austin  and  the  Texas  Trans- 
portation Institute  (TTI)  under  the  support  and  leadership  of  the  Texas 
Highway  Department.   This  project  has  provided  continued  development  and 
extension  of  the  flexible  pavement  design  and  SAMP  concepts  and  their 
implementation  for  the  design  of  flexible  pavements  in  Texas.   This  has 
resulted  in  a  series  of  Flexible  Pavement  System  (FPS)  Computer  programs 
currently  through  FPS  13  with  increasingly  sophisticated  submodels  that 
have  advanced  a  long  way  toward  rational  idealization  of  the  total  flexible 
pavement  design  problem.   A  Rigid  Pavement  System  (RPS)  was  also  developed 
at  the  University  of  Texas  [  59  ]  and  its  approach  to  optimization 
approximated  that  in  FPS. 

The  NCHRP  continued  its  work  during  the  period  1972  to  1974  through 
a  contract  with  the  Texas  Transportation  Institute  to  conduct  pilot 
implementation  tests  in  Kansas,  Lousiana  and  Florida  with  SAMP5  and  to 
improve  SAMP5  based  on  the  experience  gained  in  implementation.   Consider- 
able experience  had  already  been  gained  during  implementation  in  Texas, 
but  it  was  desired  to  test  implementation  in  other  states.   The  principal 
results  from  this  effort  as  reported  by  Lytton  and  McFarland  in  [  62  1 
was  a  somewhat  improved  computer  program  called  SAMP  6,  and  the  conclusion 
that  this  type  system  was  readily  implementable  in  other  states,  particul- 
arly those  using  the  AASHO  Interim  Design  guides. 

The  Federal  Highway  Administration  has  also  carried  out  a  research 
program  to  develop  a  mechanistic  structural  submodel  and  concepts  for 
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implementation  within  an  overall  management  system.   This  research  was 
conducted  largely  at  MIT  [1,  18,  63],   and  produced  a  computer  program 
called  VESYS  II   that  considers  pavement  materials  as  linearly  visco- 
elastic  and  predicts  their  initial  life  on  a  probabilistic  basis. 

Additional  work  has  also  been  carried  out  at  the  University  of 
Waterloo,  Canada   [64],    but  will  not  be  discussed  here. 

The  ITS  and  SAMP  pavement  systems  are  basically  the  same  since  the 
overall  pavement  system  development  in  Texas  has  evolved  through  cooperative 
efforts  and  important  improvements  in  one  system  are  generally  applied 
to  the  other.   The  choice  of  a  specific  computer  program  for  use  between 
FPS  14  and  SAMP6  for  example  would  be  primarily  a  choice  between  specific 
features  of  the  particular  programs.   From  an  optimization  viewpoint 
there  are  no  important  differences  at  all.   It  will  be  convenient  for  the 
rest  of  this  report  to  consider  FPS  and  SAMP  as  similar  systems  in  terms 
of  optimization  with  the  understanding  that  specific  computer  programs 
contain  some  variations  in  the  capabilities  of  their  submodels.   This 
group  of  concepts  and  computer  programs  will  subsequently  be  called  the 
Texas  FPS  or  FPS. 

The  goals  of  each  of  the  existing  pavement  systems  are  the  same, 
although  the  approach  toward  achievement  of  these  goals  varies  considerably. 
Eacn  of  these  systems  is  aimed  at  a  rational  design  that  will  provide 
the  riding  public  a  pavement  that  meets  their  requirements  for  comfort  and 
safety.   Fig.  7-1  describes  the  design  process  simply  in  systems  analysis 
terms.   Fig.  7-2  shows  the  conceptual  basis  for  VESYS  II  as  reported  by 
Soussou,  Moavenzadeh  and  Findakly  in   [1]»   Only  Phase  1  of  this  effort 
has  currently  been  developed. 

The  most  basic  difference  between  the  Texas  and  the  FHWA  systems 
is  that  VESYS  II  seeks  to  characterize  material  properties  as  visco- 
elastic  and  to  predict  the  stress  damage  and  life  expectancy  on  that 
basis,  while  the  Texas  FPS  characterizes  materials  on  their  resilient 
properties  after  an  appropriate  number  of  cycles  of  loading.   Both 
use  modification  of  the  equation  resulting  from  the  AASHO  Road  Test  to 
predict  the  present  serviceability  index  with  time.   Both  consider  some 
variability  of  design  parameters,  and  FPS  calculates  rehabilitation 
strategies  required  to  restore  the  serviceability  of  the  pavement  to  a 
higher  level  when  necessary  and  optimizes  the  design  choices  on  the 
basis  of  costs.   Each  method  looks  forward  to  a  development  of  sufficient 
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PHASE 


(1)  DEFINE  THE  PROBLEM 

(State  Objectives  and  Constraints) 


(2)  ESTABLISH  SYSTEM  REQUIREMENTS 

(Inputs,  Outputs,  Costs,  Decision  Rules) 


(3)  MODEL  THE  PROBLEM 

(Response  to  Loads  and  Climate, 
Performance  and  Cost  Models,  Etc.) 


(4)  GENERATE  ALTERNATE  SOLUTIONS 

(Thickness,  Materials,  Construction,  Maintenance) 


(5)  SELECT  BEST  AND  IMPLEMENT 


(6)  CHECK  PREDICTED  PERFORMANCE 


Figure  7-1  MAJOR  PHASES  OF  THE  SYSTEMS  ANALYSIS  METHOD  [1] 
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Figure  7-2   GRAPHICAL  REPRESENTATION  OF  A  PAVEMENT 

DESIGN  PROCESS  (After  [1]. 
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data  on  social  and  user  benefits  to  broaden  the  optimization  model  past 
the  cost  basis  currently  possible. 

Another  very  important  difference  between  the  two  systems  is  the  lack 
of  any  consideration  of  the  effects  of  swelling  clays  on  pavement  service- 
ability by  VESYS  II.   It  is  to  be  expected  and  has  been  clearly  proven  that 
loss  of  serviceability  for  pavement  structures  founded  on  clays  subject  to 
seasonal  or  other  volume  change  is  a  strong  function  of  random  volume 
change  as  well  as  traffic  effects.   A  means  of  including  these  effects 
in  serviceability  predictions  exists  in  the  Texas  FPS  and  means  of  quanti- 
fying input  values  from  laboratory  tests  has  been  developed   [65]. 
These  effects  cannot  be  related  to  an  assumed  elasticity  or  viscoelasticity 
in  the  subgrade  as  the  mechanism  involved  is  not  primarily  dependent  on 
loading  and  response  to  loading,  so  inclusion  of  a  subsystem  such  as  the 
one  in  FPS  will  ultimately  be  necessary  if  meaningful  serviceability 
predictions  are  to  accrue  from  VESYS  where  clay  subgrades  are  involved. 

Development  has  been  much  greater  on  the  Texas  FPS  than  on  VESYS  II. 
The  Texas  FPS  in  earlier  versions  had  used  AASHO  equations  altogether 
(no  elastic  layer  theory)  and  has  been  implemented  in  selected  districts 
of  the  Texas  Highway  Department  for  some  time.  A  more  sophisticated 
model  with  elastic  layer  theory  and  probabilistic  considerations  has  been 
developed  but  not  yet  implemented   [66],    VESYS  II  as  currently  configured 
follows  the  life,  cycle  through  a  probabilistic  prediction  of  the  life  of 
only  one  design  strategy  for  the  initially  constructed  pavement  only. 
While  conceptually  stated  in  [1],   the  models  for  generating  maintenance 
policy,  updating  serviceability  and  life  expectance,  the  cost  model, 
and  the  optimization  model  do  not  yet  exist  for  VESYS  II.   The  FPS  model 
on  the  other  hand  includes  these  capabilities  except  that  cost  is  the  only 
decision  criteria  utilized  directly  for  optimization  at  the  present  time. 
This  does,  however,  include  user  costs  for  user  delays  occurring  during 
overlay  periods. 

7.2.1  VESYS  Conceptual  Decision  Structure  for  Pavement  Maintenance 
Developed  at  MIT 

The  VESYS  concepts  for  arriving  at  maintenance  decisions  as  developed 

at  MIT  are  described  in  Appendix  9  (Part  III)  of  [1].    This  subsection 

is  intended  to  provide  a  condensed  version  of  the  most  pertinent  concepts 

described  therein. 
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The  basis  for  the  conceptual  analysis  is  the  Bayes  approach  to  the 
theory  of  decision-making,  utilizing  a  multi-attribute  utility  analysis. 
The  Markovian  property  of  probabilities  of  future  occurrences  being  dependent 
on  current  state  and  independent  of  past  events  is  assumed  to  apply. 
Strategies  yielding  maximum  utilities  or  minimum  cost/benefit  ratios 
constitute  the  optimal  strategies  to  be  selected  for  maintenance  of  the 
particular  design  configurations  under  study. 

The  factors  to  be  considered  are  "cost,  safety,  comfort,  performance, 
life,  etc."  In  recognition  of  the  subjectivity  of  some  of  the  factors, 
expected  utilities  (expressions  of  judgemental  preferences  of  individuals 
for  certain  commodity)  is  the  proposed  basis  for  analysis.   Tradeoff  is 
then  to  be  conducted  between  the  different  attributes  for  the  selection 
of  the  optimal  set  of  strategies. 

The  AASHO  serviceability  index  defines  performance  at  any  particular 
time  and  the  attributes  safety,  comfort,  and  cost  are  defined  as  decision 
criteria  in  addition  to  serviceability. 

Serviceability  may  be  increased  to  a  higher  level  at  any  time  with 
application  of  maintenance,  the  level  of  serviceability  being  dependent 
on  the  level  of  maintenance  selected.   The  three  levels  of  maintenance 
considered  are  high,  medium  and  low. 

Safety  and  comfort  may  also  be  increased  due  to  upgrading  of  the 
system  such  as  an  increase  in  skid  resistance,  overlay,  etc. 

Maintenance  level  is  taken  to  define  the  effort  (in  terms  of  labor 
and  materials)  which  is  expended  to  upgrade  the  state  of  the  system  from 
its  present  level.   A  future  cost  model  is  assumed  to  assign  costs  for  any 
maintenance  effort  discounted  to  present  worth,  providing  a  criterion  for 
alternative  evaluation. 

Fig.  7-3  shows  the  decision  flow  diagram  from  one  state  with  an 
associated  probability  to  the  next.   "M"  means  a  decision  to  maintain  and 
"N/M"  a  decision  not  to.   If  maintaining,  the  different  levels  of  main- 
tenance each  have  an  associated  cost  and  the  resulting  state  has  a  pro- 
bability of  occurring. 

The  utilities  for  cost,  performance,  safety,  etc.  will  be  evaluated 
for  each  time  period  and  these  utilities  accumulated  to  ammount  to  an 
established  utility  for  each  strategy.   These  "multi-attribute  utilities" 
are  ultimately  scaled  down  to  "monetary  equivalent"  for  optimization. 
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A  method  of  approach  is  described  in  Appendix  9  of  I  1  I   using 
the  "mid-value  splitting  technique"  to  scale  the  individual's  preference 
or  utility  for  each  attribute.   This  depends  on  personal  interviews 
between  engineers  and  prospective  users  on  the  assumptions  that  (1)  the 
description  of  the  state  of  the  highway  is  articulate  and  (2)  "the  user 
is  educated  and  is  fully  aware  of  the  consequences  of  his  needs  and 
responses."  Using  a  rating  system  from  0  to  10,  the  engineer  asks  questions 
and  raises  and  lowers  his  values  until  the  user  becomes  indifferent 
to  the  differences.   The  end  result  from  this  would  be  a  rather  complex 
multi-dimensional  "joint  preference  space".   The  evaluation  of  the  strategies 
is  illustrated  in  Fig.  7-4  which  relates  to  the  "chance  nodes"  of  Fig.  7-3 
Each  chance  node  is  represented  by  a  lottery  with  probabilities  and 
consequences  shown  at  the  top  of  each  branch.   All  multi-dimensional 
utilities  are  reduced  to  scalar  utilities  expressed  in  terms  of  perfor- 
mance indices . 

After  some  decision  in  [1]   that  in  effect  recognizes  that  the 
assumptions  may  not  be  sufficiently  valid  and  that  most  utilities  may  not 
as  yet  be  numerically  compared ,  the  authors  decided  to  scale  utilities  into 
monetary  values  and  deal  with  costs  as  the  only  attribute  directly  compared 
in  decision-making  (as  in  the  Texas  FPS  system)  pending  more  research 
on  multi-attribute  utility  analysis. 

7.2.2  Texas  FPS  System  for  Decisions  and  Optimization 

The  method  for  benefit  analysis  applied  in  the  Flexible  Pavement 
System  currently  in  use  in  Texas  is  the  "cost-effectiveness  method", 
utilizing  the  "equal-effectiveness  criterion".  That  is,  a  range  of 
design  strategies,  each  providing  a  serviceability  level  across  the 
design  life  equal  or  greater  than  some  minimum  PSI,  are  compared  on  a 
cost  basis  and  ranked  in  order  of  ascending  costs.   The  costs  considered 
by  the  economics  subsystem  for  each  strategy  include  those  for  initial 
construction,  routine  maintenance,  overlays,  user  delays  during  overlays, 
and  salvage  costs.   All  of  these  costs  are  reduced  to  present  worth  and 
combined  into  a  total  cost  for  the  design  strategy  for  use  in  ranking. 
SAMP6  also  includes  subroutines  to  include  costs  for  initial  construction 
of  shoulders  and  shoulder  overlays  concurrent  with  pavement  overlays  and 
for  overall  cross-section  effort  on  a  volumetric  rather  than  area  of  pave- 
ment surface  basis. 
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The  actual  design  decision  or  final  selection  of  strategy  Is  accom- 
plished external  to  the  computer  program  by  an  engineer.   The  basis  for 
selection  is  engineering  judgement  in  consideration  of  the  cost  compari- 
sons printed  out,  his  working  experience  with  the  various  materials  and 
their  combinations  and  a  variety  of  other  personal  (or  possibly  political) 
decision  criteria.   (See  Fig.  7-5  for  a  partial  listing.) 

FPS13,  a  version  developed  but  not  yet  implemented,  includes  a 
capability  for  considering  the  stochastic  nature  of  the  design  variables 
leading  to  design  strategies  at  differing  levels  of  confidence  or 
reliability.   Standard  deviations  or  variations  are  input  for: 

1.  Average  daily  traffic 

2.  Percent  trucks 

3.  Axles  per  truck 

A.   Temperature  constant 

5.  PSI  of  initial  structure 

6.  Axles  load/equivalency  parameter 

7.  Lack-of-fit  for: 

a.  Performance  model 

b.  Deflection  model  for  new  construction 

c.  Deflection  model  for  overlay 

8.  Swell  of  clay  subgrades 

9.  Stiffness  coefficient  of  the  subgrade 
10.   Stiffness  of  pavement  base  and  subbase. 

SAMP6  also  includes  stochastics  on  a  limited  basis  by  establishing 
a  reliability  confidence  level  requirement  and  a  single  overall  coefficient 
of  variation  on  traffic  variables. 

A  number  of  decision  criteria  have  been  defined  in  conceptual  systems 
from  various  FPS  reference  sources  including  primarily: 

1.  Costs  7.   Availability  of  funds 

2.  Safety  8.   Time 

3.  Maintainability  9.   Existing  practice 

4.  Reliability  10.   Administrative  requirements 

5.  Function  11.   Public  view 

6.  Riding  quality 

The  decision  criteria  are  conceived  as  having  weighting  functions 
assigned  by  the  design  engineer  (Fig.  7-5)  so  that  the  computer  may  make 
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selections  based  on  assessments  of  these  attributes  accruing  to  the 
separate  design  strategies.   Costs  are  specifically  considered  in  the 
Economic  Subsystem  and  reliability  may  be  established  as  a  design  para- 
meter. 

Riding  quality  is  implicitly  established  through  a  minimum  allowable 
PSI,  but  there  can  be  a  very  considerable  variation  overall  in  levels 
of  riding  quality  between  design  strategies,  as  PSI  is  not  solely  a 
function  of  riding  quality. 

Maintainability  is  implicitly  satisfied  when  a  design  strategy  is 
predicted  to  produce  a  satisfactory  level  of  serviceability  over  the 
design  life  at  an  acceptable  cost,  but  other  maintenance  characteristics 
such  as  ease  or  speed  are  not  considered.   Availability  of  funds  is 
satisfied  since  "funds  available"  is  input  as  a  design  constraint. 
"Function"  is  satisfied  by  the  input  description  of  the  traffic  neces- 
sary to  support  the  function  and  by  the  Terminal  PSI  level  selected. 

While  a  safety  subsystem  has  not  yet  been  fully  implemented  into  the 
FPS,  work  is  underway  to  relate  skid  requirements  in  terms  of  skid  numbers 
to  laboratory  tests  for  rate  of  polish  such  that  skid  numbers  may  be 
predicted  in  terms  of  traffic  applications  to  predict  decay  to  a  minimum 
acceptable  skid  number.   Loss  of  skid  number  to  the  minimum  would  trigger 
seal  coats  with  specific  materials  on  the  basis  of  their  rates  of  polish. 
The  design  strategy  would  then  consider  combinations  of  overlays  (can 
improve  skid  resistance  as  well  as  ride  quality)  and  seal  coats  to  maintain 
minimum  levels  of  PSI  and  skid  number.   While  safety  is  also  dependent 
on  roughness,  the  minimum  allowable  PSI  will  generally  exceed  the  minimum 
requirement  for  safety.   Safety  implications  of  pavement  geometries  would 
necessarily  be  considered  external  to  the  computer  program. 

Time  as  a  decision  criterion  is  not  presently  planned  for  considera- 
tion internal  to  the  computer  program,  except  to  the  extent  that  overlay 
periods  increase  user  costs  in  the  economic  model.   Such  considerations  as 
early  completion  of  initial  construction  could  be  programmed,  but  are 
currently  left  as  a  design  engineer  selection  criterion. 

7.2.3   Comprison  of  VESYS  and  Texas  FPS  Concepts  for  Decision  and  Optimization 

The  comparison  of  the  two  concepts  is  more  specifically  the  compar- 
ison of  the  conceptual  development  for  VESYS  as  described  and  of  an 
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existing  operable  optimization  system  planned  for  expansion  in  the  future 
to  consider  attributes  in  addition  to  those  now  considered. 

While  the  FPS  as  implemented  includes  some  decision  criteria  implicitly 
as  previously  described  and  optimizes  on  costs  (including  user  costs 
during  overlays),  VESYS  cannot  presently  generate  design  strategies  and 
the  optimization  features  exist  only  conceptually. 

The  two  approaches  to  decision  making  and  optimization  appear  to 
vary  greatly  in  theory.   The  VESYS  approach  would  generate  a  "tree" 
with  many  branches  with  nodes  occurring  at  succeeding  increments  in  time, 
while  the  FPS  system  in  effect  only  creates  nodes  in  the  tree  when 
serviceability  reaches  the  minimum  allowed.   Both  would  evaluate  the 
resulting  attributes  for  separate  design  strategies  at  the  "tips  of  the 
decision  tree"  or  at  the  end  of  the  design  life  period. 

The  approach  to  considering  stochastic  variations  is  also  some- 
what different  in  the  two  methods  and  could  result  in  different  selections, 
The  probability  of  attaining  the  condition  state  under  consideration  is 
used  as  a  multiplier  at  each  node  of  the  VESYS  decision  tree  to  reduce  the 
anticipated  benefit  in  relation  to  its  uncertainty  of  attainment.   The  Texas 
FPS,  on  the  other  hand,  establishes  a  required  reliability  level  as  a 
constraint  and  eliminates  design  strategies  of  lesser  reliability  from 
consideration.   The  reliability  for  a  particular  design  strategy  is  a 
function  of  coefficients  of  variation  or  standard  deviation  input  for 
various  design  parameters  as  derived  from  field  data,  testing  analysis, 
and  planned  level  of  construction  control.   This  procedure  is  discussed 
in  [66]   by  Darter  and  Hudson. 

Although  a  direct  comparison  would  be  difficult  without  consider- 
able development  work  on  the  VESYS  concepts,  it  seems  possible  that 
similar  design  strategies  might  result  from  the  two  if  initial  condition 
states  were  the  same,  decisions  to  maintain  for  VESYS  were  limited  to  nodes 
when  serviceability  approached  the  same  minimum,  condition  states  after 
overlay  were  the  same  and  if  the  probabilities  assigned  to  VESYS  approx- 
imated the  confidence  levels  and  variations  input  to  FPS. 

Both  systems  consider  the  same  type  of  decision  criteria  (sometimes 
called  benefits,  utilities,  attributes,  etc.,  in  the  literature)  and 
neither  has  advanced  toward  a  serious  total  analysis  of  all  decision 
criteria,  primarily  because  meaningful  bases  for  quantifiying  them  for 
comparison  do  not  yet  exist. 
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7.3   CONCEPTUAL  INTERRELATIONSHIP  OF  SYSTEM  COMPONENTS 

The  development  of  a  pavement  design  system  involves  as  one  task  the 
organization  of  the  system  components  or  subsystems  in  a  rational  and 
efficient  manner.   Each  must  not  only  satisfy  its  specified  functions 
rationally  and  efficiently,  but  It  must  relate  to  and  interact  with  other 
subsystems  so  that  the  overall  system  objectives  are  achieved. 

As  each  subsystem  has  numerous  components  and  is  complex  within 
itself,  the  level  of  complexity  and  abstraction  as  the  entire  system  is 
assembled  exceeds  the  comprehension  of  the  human  brain  unless  some  tool 
is  used  to  allow  it  to  focus  on  a  comprehensive  amount  of  operations  at 
one  time  without  loss  of  the  rest  of  the  detail.   Interaction  and  flow 
diagrams  are  used  to  graphically  show  the  structure  of  a  system  and  to 
retain  it  generally  intact  while  a  subsystem  is  studied  and  possibly 
revised.   These  diagrams  may  be  simple  with  entire  subsystems  shown  in 
one  block  or  in  great  detail,  depending  on  their  intended  use.   Such 
diagrams  to  describe  conceptual  interrelationships  of  the  system  com- 
ponents for  a  pavement  analysis  and  design  system  require  more  detail  than 
in  Fig.  7-1  and  7-2  but  not  so  much  that  they  cannot  be  easily  comprehended. 

A  simplified  diagram  of  a  pavement  design  system  appears  in  Fig  .  7-5. 
The  process  begins  with  inputs  and  continues  along  the  paths  of  the  arrows. 
In  this  very  general  diagram,  the  structural  subsystem  generates  alternate 
design  strategies  including  overlay  strategies  and  predicts  serviceability 
levels  with  time.   The  various  design  strategies  are  checked  for  safety 
in  terms  of  skid  resistance  and  roughness,  user  delays  during  overlay 
operations  are  evaluated  and  costed,  total  costs  are  calculated  for  the 
various  design  strategies  and  the  results  are  arrayed  by  total  cost  and 
printed  out.   The  decision  criteria  block  is  shown  as  a  question  mark 
due  to  current  inability  to  apply  all  the  major  criteria  effectively. 
A  design  in  this  case  is  selected  by  the  engineer  after  reviewing  the  cost 
optimization  by  the  system  and  applying  the  remaining  decision  criteria 
himself.   At  this  point,  the  design  is  complete,  but  it  is  most  important 
that  the  process  be  closed  by  a  "data  feedback  loop"  to  allow  future 
comparisons  of  actual  responses  to  predicted  responses  so  that  the  sub- 
systems may  be  continually  upgraded  to  improve  their  predictive  ability. 
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KcFarland   [67]   has  developed  a  very  approximate  means  of  evaluating 
user  costs  during  overlays  and  reviews  the  pros  and  cons  of  five  methods 
of  benefit  analysis.   The  method  used  in  the  Texas  FPS  seems  to  be  as  good 
as  or  better  than  any  other  available  until  a  capability  for  meaningful 
quantification  of  subjective  benefits  is  developed. 

The  method  of  approach  proposed  in   [1]  for  establishing  utility 
values  is  theoretically  interesting,  but  is  not  believed  to  be  practically 
implementable  as  described.   It  suffers  from  lack  of  the  relationship 
to  actual  physical  roadways  that  the  usual  "rating  sessions"  have  enjoyed. 
Also,  "utilities"  for  ride  quality  have  already  been  established  at  great 
expense  and  efforts  are  well  along  to  characterize  safety  in  terms  of 
skid  number  and  roughness . 

While  both  pavement  systems  in  concept  approach  optimization  on  the 
basis  of  selecting  the  best  overall  design  strategy  from  an  array  of 
design  strategies  that  have  been  developed,  the  practicality  of  this 
approach  is  strongly  related  to  the  computer  time  required  in  the  parti- 
cular system  for  developing  a  strategy.  As  VESYS  II  utilizes  a  series 
of  solutions  at  points  in  time  to  arrive  at  its  life  predictions,  the 
computer  time  required  for  generating  an  array  of  design  strategies  is 
expected  to  be  many  times  greater  than  for  FPS.  While  FPS  as  implemented 
may  consider  several  hundred  or  even  thousands  of  design  strategies  to 
reach  a  selection  of  the  ten  or  twenty  best  for  consideration,  only  a 
small  fraction  of  that  many  would  be  practical  for  VESYS  II  and  very 
possibly  for  more  sophisticated  versions  of  FPS  as  well.   It  appears  then 
that  VESYS  II  will  require  optimization  within  the  submodels  to  reduce 
considerably  the  array  of  design  strategies  produced  without  loss  of 
generality. 

It  must  be  remembered  that  the  models  leading  to  serviceability 
predictions  include  a  mix  of  theory  and  empirical  developments  that  range 
from  fairly  gross  to  reasonably  accurate  approximations.   Research  in 
progress  will  result  in  improvements  that  must  be  considered.   This 
expectation  of  change  emphasizes  the  necessity  for  a  formalized  data 
feedback  system  as  discussed  in  detail  in  [68]. 
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It  should  be  noted  that  the  decision  criteria  may  be  applied  as  input 
constraints  or  on  a  comparative  basis  after  all  design  strategies  have  been 
processed  or  as  a  combination  of  the  two. 

having  studied  the  simplified  conceptual  system  of  Fig  •  7-5  additional 
detail  may  now  be  considered  as  in  Fig.  7-6.  Note  that  some  decision 
criteria  will  be  applied  as  constraints  and  that  stochastic  concepts  are 
now  included  to  account  for  variability  in  the  input  parameters. 
Probability  distributions  are  output  for  all  the  predicted  responses. 
This  concept  does  not  establish  the  engineer's  specific  involvement  with 
optimization.   He  may  apply  the  decision  criteria  himself  or  share  the 
application  with  the  system. 

The  basic  building  block  to  be  developed  and  expanded  is  the  VESYS  IIM 
system,  which  at  present  has  only  been  developed  through  Phase  I  as 
depicted  in  Fig.  7-2.   This  capability  is  represented  by  part  of  the  inputs 
and  the  structural  subsystem  in  Fig.  7-5,  except  that  VESYS  II  at  present 
does  not  generate  alternate  designs  and  considers  only  the  initial  constructed 
pavement  (no  overlays  or  restorative  maintenance) .   The  relationships 
betweem  the  components  of  VESYS  II  are  portrayed  in  Fig.  7-7.  Material 
properties  fitted  to  a  mathematical  relation,  the  geometry  (primarily 
layer  thickness,  keyed  to  material  properties)  and  the  radius  of  the 
applied  lead  are  used  in  a  static  elastic  layer  analysis  to  arrive  at  the 
primary  responses  of  radial  strains  and  vertical  displacements  directly 
beneath  a  unit  load.   These  responses  are  curve  fit  with  relation  to 
means  and  variances  and  become  input  to  a  repetitive  load  subprogram. 

Using  1)  inputs  from  the  static  unit  load  analysis,  2)  temperature 
data  including  a  system  shift  function  related  to  the  material,  3)  the 
initial  serviceability  and  4)  the  minimum  allowable  serviceability,  the 
repetitive  load  subsystem  applies  fatigue,  rut  depth  and  roughness  models 
to  arrive  at  predictions  of  roughness,  rut  depth  and  cracking.   This 
data  is  then  used  in  a  serviceability  prediction  subprogram  to  produce 
predictions  of  service  life  and  reliability. 

As  there  is  nothing  to  optimize  until  multiple  design  strategies 
are  generated,  optimization  of  VESYS  IIM  implies  its  expansion  to  essenti- 
ally embody  those  components  or  subsystems  and  capabilities  shown  in 
Fig.  7-5  and  7-6. 
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7.4   SPECIFIC  DECISION  CRITERIA 

Decision  criteria  for  a  pavement  system  should  include  all  those 
characteristics  or  attributes  that  have  impact  on  the  safety,  cost  of  or 
benefits  available  from  the  pavement  system  and  that  are  subject  to  control. 
Some  decisions  may  be  reached  internal  to  the  design  system  through 
establishment  of  constraints  (usually  as  input  data)  such  as  maximum 
allowable  cost  and  minimum  time  between  overlays.   Others  must  be  made  by 
comparison  of  the  "outcomes"  or  the  characteristics  of  the  design  strategies 
predicted  by  the  pavement  analysis  and  design  system. 

The  following  are  decision  criteria  identified  by  study  of  the 
literature  with  some  discussion  of  each: 

1.  Costs  -  in  a  very  general  sense,  costs  could  include  vehicular 
wearout  and  damage  to  sensitive  freight  as  a  function  of  lateral 
and  longitudinal  roughness,  tire  wear  as  a  function  of  surface 
texture  and  many  other  such  parameters  as  well  as  those  costs 
previously  identified  as  being  considered  in  the  FPS.   The  reason 
such  parameters  are  omitted  is  that  insufficient  data  is  available 
upon  which  to  make  useful  comparisons.   Even  the  user  cost  data 
developed  by  McFarland   [67]   for  the  Texas  FPS  is  based  on 
limited  data  although  part  of  it  is  in  use.   The  cost  parameters 
that  may  reasonably  be  compared  with  present  data  are  initial 
construction  cost,  routine  maintenance  cost,  cost  for  overlays 
and  seal  coats,  user  costs  due  to  delays  during  overlay  construc- 
tion, and  the  negative  cost  of  salvage  value  after  the  design  period. 

McFarland   [67]   developed  costs  for  vehicle  operation, 
rural  accidents  and  passenger  discomfort  as  functions  of  pavement 
serviceability  level,  but  considered  them  too  limited  for  use 
pending  further  research.   One  limitation  was  that  the  user  costs 
were  related  only  to  PSI  rather  than  to  other  parameters  also 
affecting  them  such  as  type  of  road  and  skid  number.   The  other 
important  limitation  is  the  lack  of  sufficient  reliable  data  upon 
which  to  base  the  benefit  model. 

2.  Safety  -  Ignoring  the  safety  aspects  of  geometries  that  are  not 
considered  in  the  pavement  analysis  and  design,  safety  may  be 
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considered  in  terms  of  the  effects  of  the  pavement  surface  on 
ability  of  a  vehicle  operator  to  control  his  vehicle  in  dry  or 
wet  weather.   Specifically,  skid  resistance  may  be  related  to  a 
skid  number,  which  is  primarily  related  to  the  initial  skid 
resistance  of  the  surface  materials  and  the  rate  of  polish  and 
consequent  decay  of  skid  resistance  with  traffic  applications. 
Rutting  and  longitudinal  roughness  affect  the  driver's  control, 
but  only  seriously  when  rutting  and  roughness  are  severe. 

Substantial  research  is  underway  to  learn  about  skid  resis- 
tance and  to  relate  it  to  safety  and  types  of  surfaces.   Other 
research  is  underway  to  relate  decay  of  skid  number  to  traffic 
applications  and  to  laboratory  testing  to  establish  rates  of 
polish  of  aggregates.   These  studies  may  be  expected  to  allow 
development  of  an  ability  to  predict  decay  of  skid  number  with 
time  for  a  specific  surface  and  traffic  estimate  much  as  service- 
ability is  now  predicted  using  the  AASHO  equation.   The  imple- 
mentation of  this  capability  into  the  Pavement  Analysis  and 
Design  System  (PADS)  will  require  optimization  of  combinations 
of  overlays  that  may  restore  serviceability  and  skid  resistance 
and  intermediate  seal  coats  to  restore  skid  resistance  only. 
The  requirement  will  then  be  to  maintain  irinimum  established 
levels  of  both  serviceability  and  skid  resistance  with  trade- 
offs between  performance,  safety  and  other  benefits. 

It  is  believed  that  the  safety  aspects  of  rutting  and 
longitudinal  roughness  will  be  implicitly  satisfied  by  meeting 
ride  quality  standards. 

Maintainability  -  is  defined  in   [  70,.  71  ]    as  a  measure  of  the 
degree  to  which  effort  may  be  required  during  the  service  life 
to  keep  serviceability  at  a  satisfactory  level.   This  includes 
normal  or  routine  maintenance  such  as  repairing  random  potholes  and 
sealing  random  cracks  and  repair  maintenance  required  to  restore 
adequate  serviceability  when  a  serious  loss  has  occurred  or  seems 
imminent . 

Alexander   [70]    proposes  an  "index  of  maintainability", 
which  amounts  to  the  arbitrary  number  1000  divided  by  predicted 
maintenance  cost  to  avoid  fractions.   A  high  index  number  in  this 
case  implies  low  maintenance  costs.   These  indices  could  be  used 
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as  numerical  measures  of  maintainability  for  comparison  as 
decision  criteria. 

Alexander  also  proposes  a  probabilistic  expected  cost  analysis 
basis  as  follows: 


i=n 
i=l 


=  I       (P  (X,)  x  Vt) 


where: 

EC  =  the  expected  cost  of  a  situation  that  may  have 

any  one  of  n  outcomes 

p(X.)  =  probability  that  the  outcome  will  be  X 

V.  =  cost  of  outcome  X. 
1  1 


While  the  definitions  above  imply  a  comparison  of  the 
maintenance  required  for  various  design  strategies,  a  design 
strategy  must  practically  meet  some  established  rules  to  compete 
as  to  maintainability  and  these  rules  will  vary  in  different 
environments  and  traffic  situations.   For  instance,  the  freeze- 
thaw  cycles  in  the  Northeast  cause  a  generally  higher  level  of 
acceptable  maintenance  than  in  the  Southwest,  but  an  expansive 
clay  subgrade  in  the  Southwest  would  also  markedly  increase  the 
maintenance  frequency  expectation.   Frequent  maintenance  on  a 
major  urban  freeway  would  be  less  pallatable  than  the  same 
frequency  of  maintenance  on  a  rural  highway. 

Some  of  the  rules  related  to  maintenance  used  in  FPS  as 
input  constraints  follow: 

a.  Minimum  period  before  the  first  overlay  is  required. 
Design  strategies  that  are  predicted  to  deteriorate  to 
the  minimum  established  serviceability  in  a  lesser  time 
than  that  allowed  are  eliminated  from  further  consid- 
eration. 

b.  Minimum  period  between  overlays. 

c.  Time  to  first  seal  coat. 

d.  Time  between  seal  coats. 

e.  Allowable  thickness  each  layer. 

f.  Minimum  overlay  thickness. 

g.  Maximum  total  thickness  of  all  overlays. 

h.  Maximum  funds  available  for  initial  construction. 
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Another  rule  suggested  by  Alexander   |7()|    is  that  required 
maintenance  should  not  exceed  that  amount  that  can  practically 
be  made  available. 

With  the  growing  shortages  of  some  aggregates  and  other 
materials,  constraints  on  excessive  uses  of  materials  in  short 
supply  for  maintenance  might  be  appropriate  in  some  locales. 

It  appears  that  maintainability  may  best  be  insured  for  the 
present  through  input  constraints  that  insure  that  candidate 
design  strategies  satisfy  the  rules.  There  is  a  possibility  that 
a  pavement  requiring  frequent  maintenance  might  have  other 
attributes  warranting  its  consideration,  but  the  probability  is 
so  low  that  the  inefficiency  of  processing  these  marginal  strategies 
through  the  system  does  not  appear  warranted.   In  general,  excessive 
maintenance  will  generate  so  much  construction  and  user  costs  that 
the  strategy  in  question  cannot  be  a  strong  contender.  However, 
if  desirable,  maintainability  could  be  developed  along  with  other 
decision  criteria  for  utility  optimization  in  the  future. 
4.   Reliability  is  defined  in  [71]   as  a  measure  of  the  proba- 
bility that  serviceability  will  be  at  an  adequate  level  throughout 
the  design  service  life.   Darter  and  Hudson   [66]   extended 
this  definition  to  make  it  more  complete  as  follows: 

"Reliability  is  the  probability  that  the  pavement  system  will  perform 
its  intended  function  over  its  design  life  (time)  and  under  the 
conditions  (or  environment)  encountered  during  operation." 

Reliability  may  also  be  viewed  as  the  confidence  level  that 
the  pavement  system  will  perform  in  accordance  with  its  stipulated 
design  functions. 

Virtually  every  aspect  of  pavement  design,  construction, 
operation  and  maintenance  is  subject  to  uncertainty.   Traffic 
estimates,  environmental  predictions,  anticipated  material  char- 
acteristics, the  design  models  themselves,  costs  for  construction 
and  maintenance  and  many  more  variables  are  approximations  at  best 
and  educated  guesses  in  some  cases.   While  many  are  subject  to 
some  degree  of  control,  even  they  are  sufficiently  variable  that 
probabilities  could  appropriately  be  assigned  to  all  the  design 
parameters.   The  limitation  is  lack  of  data  upon  which  to  assign 
reasonable  values  of  means  and  variances  for  use  with  the  usually 
assumed  normal  probability  distribution. 
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Increases  in  reliability  may  generally  be  obtained  at  addi- 
tional cost  by  use  of  better  quality  materials,  closer  control 
of  material  variations,  providing  more  maintenance  and  increasing 
pavement  layer  thicknesses  in  general   [66].    Fig.  7-8  illu- 
strates conceptually  the  relationship  between  reliability,  per- 
formance and  costs.   As  can  be  seen,  practical  limits  exist 
past  which  little  additional  serviceability  or  reliability  may 
be  obtained  regardless  of  additional  expense. 

There  are  two  approaches  that  may  be  taken  in  the  application 
of  reliability  as  a  decision  criterion: 

a.  It  may  be  established  at  an  appropriate  level  for  the 
system  considered  and  input  as  a  constraint,  thus 
eliminating  design  strategies  of  lesser  reliability 
(this  capability  already  exists  in  FPS)  or 

b.  It  may  be  defined  in  time  and  considered  on  a  compar- 
ative basis  along  with  other  decision  criteria. 

Lemer  and  Moavenzadeh   [71]    offer  a  very  general  mathe- 
amtical  expression  for  failure  in  terms  of  reliability  and  N 
possible  failure  modes: 


DL  (e1,  e2,  ...,  eL)  £  R^c^  c2>  ...,  cM) 


D2  (el'  *  * •'  eL*  *  R2  (cl'  "  * '  CM* 


DN  (V  •••'  eL}  *  \    (C1>  •••'  CM} 


where: 

D.  =  environmental  loads  placed  on  the  pavement  as 
determined  by  a  set  of  environmental  qualities 


(e  ,  e_ ,  ...,  e  )  such  as  wheel  loads,  temperature 
and  vehicle  speeds 


R.  =  capacity  of  the  pavement  to  resist  environmental 

load  D  .  . 

l 

c.  =  system  characteristics  that  determine  pavement  response 
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Figure  7-8  ILLUSTRATION  OF  RELATIONSHIPS  BETWEEN 
RELIABILITY  (R) ,  PERFORMANCE  (P),  AND 
COST  (C)  (After  [66]). 


162 


Failure  is  the  condition  in  which  one  or  more  of  the  inequali- 
ties D.  ^  R.  is  not  satisfied.   The  reliability  is  then  the 
probability  that  all  the  N  inequalties  are  true  and  may  be  expressed 
as: 


R  =  p(no  failure)  =  p(D.  -  R.) 
for  i  =  1,  . . . ,  N 

In  FPS,  a  much  more  specific  mathematical  definition  is  used: 

R  =  p(N  >  n) 

wh  ere : 

N  =  number  of  load  applications  that  a  section  of  pave- 
ment can  withstand  before  minimum  allowable  service- 
ability is  reached  within  a  limited  maintenance 
input . 

n  =  number  of  load  applications  which  are  applied  to  the 
pavement  section. 

This  definition  of  reliability  is  developed  from  the  basic  concept 
that  a  no-failure  probability  exists  when  the  number  of  load 
repetitions  to  minimum  acceptable  serviceability  is  not  exceeded 
by  the  number  of  load  applications  applied  n. 

In  FPS,  R  is  taken  to  represent  the  number  of  0.2  mile  sections 
that  will  survive  the  design  period  and  it  is  assumed  that  an  over- 
lay will  be  applied  when  1-R  sections  have  failed.  After  an  over- 
lay, each  section  is  assumed  to  have  an  equal  R  chance  of  success 
during  the  next  performance  period  and  the  same  for  still  other 
performance  periods. 

A  choice  of  six  reliability  levels  from  50  to  99.99  percent 
are  available  in  FPS.   Selection  of  a  reliability  level  is  based  on 
the  functional  classification  of  the  highway  traffic  volume  and 
18  kip  equivalent  single  axle  load  applications.   Table  7-1  p,ives 
recommendations  by  Darter  and  Hudson   [66],   where  C,  D,  and  E 
represent  confidence  or  reliability  levels  of  95,  99  and  99.9 
percent,  respectively. 
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TAJiLE  7-1  RECOMMENDED  DESIGN  RELIABILITY 

LEVELS  FOR  FPS  PROGRAM  (After  [66]) 
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Referring  again  to  the  possible  approaches  of  implementing 
reliability  as  a  decision  criterion,  it  appears  much  more  practical 
to  utilize  it  as  a  constraint  for  the  present  until  a  broader 
capability  exists  for  assessing  and  applying  utilities  from  a 
full  range  of  benefits.   As  the  state-of-the-art  advances,  a 
basis  for  applying  reliability  as  a  decision  criterion  should  be 
developed  along  with  other  weighted  decision  criteria. 

Reliability  will  increase  in  complexity  as  other  decision 
criteria  are  applied.   For  instance,  FPS  considers  in  effect  one 
failure  mode  (or  N  =  1  in  the  Lemer  and  Moavenzadeh  formulation) . 
As  it  becomes  possible  to  qualify  and  define  safety  requirements 
(in  terms  of  material  characteristics,  traffic  applications  and 
other  parameters) ,  a  second  failure  mode  or  inequality  must  be 
considered.   If  other  rules  for  failure  are  applied,  other 
inequalities  must  be  considered  and  probabilities  defined  for 
each. 

Reliability  is  only  meaningful  to  the  extent  that  the  proba- 
bility estimates  themselves  are  reliable.   Considerable  work 
has  been  done   [66]   to  learn  the  sensitivity  of  design  life 
predictions  to  variations  in  design  parameters  from  actual  data 
collected  in  Texas  and  elsewhere.   Much  more  research  will  be 
required  to  develop  more  data. 
5.   Function  of  a  pavement  system  involves  the  purpose  it  is 

expected  to  serve.   For  instance,  it  may  be  a  high  speed  express- 
way, an  access  road,  a  farm-to-market  road,  an  interstate  highway, 
a  collector  street  in  a  city  street  system  or  others. 

The  varying  functional  requirements  of  the  users  for  different 
types  of  highway  or  street  systems  implies  variations  in  perfor- 
mance requirements  and  benefits.   For  instance,  the  amount  of 
tolerable  variations  in  vertical  profile  (or  roughness)  for  a 
high-speed  expressway  would  be  much  less  than  that  for  a  neighbor- 
hood street.   Maintenance  frequency  may  be  greater  on  a  rural 
road  with  little  traffic  than  an  interstate  highway  with  thousands 
of  vehicles  per  day  to  be  delayed. 

It  appears  that  the  important  effects  of  pavement  system 
functions  are  "imbedded"  in  other  decision  criteria  and  that  it 
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will  serve  as  an  implicit  decision  criterion  through  its  effect 
on  them.   The  other  decision  criteria  for  which  function  is  a 
consideration  are: 

a.  Safety  -  minimum  skid  number  and  allowable  variations 
in  lateral  and  longitudinal  profile  will  be  a  function 
of  traffic  velocities. 

b.  Maintainability  is  dependent  on  the  level  of  service- 
ability required,  which  is  a  variable  with  function. 

c.  Reliability  levels  selected  should  reflect  function  in 
terms  of  the  impact  that  failures  and  consequent  main- 
tenance would  have  on  satisfactory  support  of  functional 
requirements. 

d.  Riding  quality  in  terms  of  allowable  minimum  serviceability 
levels  is  directly  dependent  on  the  function  the  pave- 
ment is  to  serve. 

e.  Availability  of  funds  is  directly  dependent  also  on  the 
function  the  pavement  is  to  serve.   Much  higher 
funding  is  required  and  justifiable  for  an  expressway 
to  serve  masses  of  vehicles  than  for  a  small  collector 
street. 

6.   Riding  quality  primarily  expresses  the  degree  of  comfort  or 

lack  of  discomfort  experienced  by  a  highway  user  when  riding  in 
an  average  automobile.   It  is  represented  in  the  PADS  by  the 
by  the  Present  Serviceability  Index  (PSI)  calculated  at  any  time 
by  modified  AASHO  equations  in  terms  of  rut  depth,  cracking  and 
slope  variance  (and  also  characteristics  of  a  clay  subgrade  in 
FPS) .   While  this  is  a  reasonable  parameter  for  interim  use,  it 
should  be  understood  that  the  data  from  the  AASHO  Read  Test  are 
basically  50  m.p.h.  data  and  that  discomfort  is  dependent  on  a 
amplitudes  of  roughness  and  frequency  of  occurrence,  which  depends 
on  vehicle  velocity.   That  is,  the  same  profile  on  a  30  m.p.h. 
city  street  and  a  55  m.p.h.  expressway  would  not  be  expected  to 
yeild  the  same  PSI.  rating  by  the  same  raters  in  the  same  auto- 
mobiles traveling  at  the  different  velocities. 

Riding  quality  may  be  implicitly  enforced  by  input  of  a 
minimum  allowable  PSI  as  in  FPS,  but  should  also  be  considered 
in  terms  of  performance  across  the  design  period.   This  would 
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allow  selection  of  the  design  strategy  with  higher  overall  perfor- 
mance if  costs  and  other  benefits  were  approximately  equivalent. 
Performance  may  be  measured  as  the  area  under  the  Pavement 
Serviceability  Index  vs.  Time  Curve: 


final 


Performance  =  j        PSI(t)dt 


Potential  performance  past  the  end  of  the  design  period 
should  also  be  considered  as  a  part  of  the  salvage  value  or  in 
some  other  way  to  take  this  benefit  into  account  during  optimi- 
zation.  The  degree  of  importance  to  be  attached  to  this  secondary 
benefit  could  be  controlled  by  an  input  weighting  function.   This 
benefit  could  be  considered  separately  from  ride  quality  if 
desired. 

7.  Availability  of  funds  for  initial  construction  is  presently 
imposed  in  FPS  as  a  constraint  on  the  allowable  initial  con- 
struction cost  of  the  pavement  system.   Design  strategies  are 
eliminated  that  require  more  funding  to  construct  than  that 
available.  While  this  appears  generally  reasonable,  there  is 
a  good  probability  that  the  total  costs  over  the  design  period 
might  be  considerably  less  for  some  design  strategies  eliminated 
because  of  initial  construction  costs  than  for  the  optimum 
strategy  selected.   Knowledge  of  such  an  advantage  might  warrant 
some  realignment  of  funding  to  attain  it. 

Study  might  be  warranted  to  determine  whether  high  initial 
cost  strategies  should  be  considered  or  eliminated  as  in  FPS. 

The  availability  of  funds  for  overlay  and  seal  coat  might 
also  warrant  consideration  in  selection  between  design  strategies 
having  fewer  increments  of  costlier  major  repair  or  those  having 
more  increments  of  cheaper  major  repair. 

8.  Time  is  designated  in  the  literature  as  a  decision  criterion. 
It  is  indeed  a  very  important  parameter  and  appears  as  input 
constraints  (time  between  overlays,  etc.),  as  a  parameter  for 
calculating  other  design  criteria  (such  as  user  delay  costs) 
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and  as  a  calculated  output  (such  as  time  for  PSI  to  deteriorate 
to  minimum  allowable) . 

Although  not  currently  planned,  it  could  be  considered  as 
an  independent  design  criterion  for  possible  future  design  criteria 
such  as  initial  construction  time  and  maintenance  periods  on  a 
utility  basis  rather  than  user  costs. 
9.   Existing  Practice  is  the  empirical  result  of  years  of  good  and 

bad  experience  with  various  materials  and  designs  over  a  range  of 
conditions.   It  exists  in  codes,  design  manuals,  research  reports, 
textbooks,  in  the  minds  if  practicing  engineers  and  in  other 
locations. 

It  would  be  very  difficult  to  devise  a  means  of  applying 
existing  practice  as  a  decision  criterion  internal  to  the  PADS, 
but  it  may  be  applied  externally  by  the  design  engineer  as  he  reviews 
the  design  strategies  and  makes  his  selection.   Kenis  and  McMahon 
[18]    propose  initial  use  of  the  VESYS  II  as  a  design  check 
procedure  for  the  various  states'  present  design  procedures. 
Such  an  interaction  between  practice  and  the  new  PADS  may  result 
in  changes  in  the  PADS  that  would  amount  to  an  imbedded  consider- 
ation of  the  useful  knowledge  reflected  in  the  existing  practice. 
10.   Administrative  requirements  -  Within  any  public  agency  there  has 
built  up  a  set  of  administrative  rules  and  procedures  which 
govern  one  or  more  phases  of  the  operation  of  the  agency.   The 
rules  may  be  derived  from  state  or  federal  laws  as  part  of  the 
system  of  rules  associated  with  a  particular  funding  source  such 
as  the  governing  Eoard  of  the  Agency;  i.e.,  the  State  Highway 
Engineers  or  the  Highway  Commission. 

These  rules  may  deal  with  political  issues  or  issues  of 
special  local  concern  such  as  "buy  American"  policies,  which 
dictate  the  use  of  local  products  in  lieu  of  imported  substitutes. 
Energy  sensitive  issues  may  also  be  treated  in  this  light  until 
they  can  become  quantifiable  in  direct  "cost"  terms  in  the  program. 
In  any  event,  these  administrative  rules  must  be  considered 
implicitly  or  explicitly  in  the  decision  process  of  pavement 
design. 
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11.   Public  view  -  the  condition  of  a  pavement  system  sometimes  becomes 
a  political  issue  with  strong  expression  of  public  dissatisfaction. 
Such  a  situation  might  warrant  replacement  of  a  particular  pavement 
at  a  high  level  of  initial  serviceability  and  perhaps  maintenance 
of  a  relatively  high  level.   Other  kinds  of  public  pressure  might 
also  be  brought  to  bear. 

The  implementation  of  this  decision  criterion  would  be  most 
complex  in  any  manner  other  than  as  consideration  with  other 
decision  criteria.   Therefore,  it  is  proposed  that  the  public 
view  be  taken  into  account  implicitly  rather  than  as  a  separate 
decision  criterion. 

Another  contemporary  criterion  is  availabilty  of  materials. 
As  many  aggregates  and  asphaltic  materials  are  in  short  supply 
around  the  country,  it  is  appropriate  to  establish  a  vehicle  for 
utilizing  the  materials  in  short  supply  only  when  necessary. 
As  an  example,  let  us  say  that  two  alternatives  for  base  material 
are  to  be  considered  and  the  total  cost  and  other  benefits  are 
about  equal.   It  is  obviously  advantageous  to  consume  from  the 
largest  supply  rather  than  the  limited  one.   The  relative  impor- 
tance of  such  conservation  of  resources  in  short  supply  could  be 
established  by  an  input  weight  for  use  with  this  design  criterion. 
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7.5   OPTIMIZING  TECHNIQUES 

A  complete  Pavement  Analysis  and  Design  System  (PADS)  implies  a 
:'black  box"  such  as  that  illustrated  in  Fig.  7-9  that  can: 

1.  Accept  reaJ  life  inputs  on  a)  available  materials  characteristics 
and  costs,  b)  traffic  and  load  predictions,  c)  performance  require- 
ments, d)  environment  characteristics,  e)  construction  quality 
control  expectations  and  f)  means  and  variances  for  appropriate 
input  parameters. 

2.  Analyze  efficiently  a  range  of  material  combinations  and  thick- 
nesses on  a  probabilistic  basis  to  arrive  at  reasonable  predictions 
of  pavement  responses  (roughness,  rutting,  cracking  and  loss  of 
skid  resistance)  with  time  and  traffic  applications. 

3.  Translate  pavement  response  predictions  into  sufficiently  realistic 
predictions  of  performance  and  reliability  for  the  various  design 
strategies  in  terms  of  both  ride  quality  and  safety. 

4.  Generate  appropriate  maintenance  policies  (generally  comprising 
overlays  and  seal  coats  of  various  thicknesses  and  using  various 
materials)  for  consideration  with  the  various  initial  designs. 

5.  Analyze  and  predict  responses  and  performance  on  a  probabilistic 
basis  for  the  upgraded  pavement  system. 

6.  Repeat  Item  5  as  required  to  provide  the  required  performance 
over  the  design  period. 

7.  Calculate  the  costs  reduced  to  present  worth  for  initial  construc- 
tion, routine  maintenance,  overlays  and  seal  coats,  user  costs 
during  maintenance,  and  the  negative  salvage  value  of  the  materials 
at  the  end  of  the  design  period. 

ft.   Assess  on  a  numerical  basis  the  utility  cf  each  overall  design 

strategy  for  each  established  design  criterion. 
9.   Make  comparisons  between  the  design  strategies  on  the  basis  of 

combinations  of  numerically  defined  and  weighted  utilities  for  the 
established  design  criteria. 
10.   Display  (printout)  the  results  of  the  analysis  and  evaluation  of 

the  various  design  strategies  so  that  a  final  selection  and  design 
implementation  may  be  made  by  the  design  engineer  or  administrator. 
VESYS  II  is  capable  of  the  first  three  items,  except  that  it  only 
handles  one  design  strategy  at  a  time  on  an  input  basis  and  does  not  consider 
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Figure  7-9   BLACK  BOX  SYSTEM  (After  [72]). 
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effects  of  subgrade  volume  change  or  skid  resistance  at  all.   A  later 
version  of  FPS  developed  but  not  yet  implemented  is  capable  of  the  first 
eight  items,  except  it  does  not  yet  consider  skid  resistance  and  optimi- 
zation is  only  on  costs.   The  accuracy  of  predictions  for  either  depends 
en  the  validity  of  its  materials  characterizations,  other  input  parameters, 
and  its  models,  and  these  have  yet  to  be  fully  tested. 

The  task  to  be  accomplished  is  not  simply  one  of  optimization,  but 
includes  extension  of  the  analysis  model  to  produce  a  general  set  of 
design  strategies  for  initial  construction,  overlays  and  seal  coats  instead 
of  one  very  specific  analysis  and  set  of  probabilistic  predictions. 

The  general  process  of  system  analysis  and  optimization  as  defined 
in   [73]   and  illustrated  in  Fig.  7-10  includes  essentially  the  following 
steps: 

1.  The  outcomes  for  each  possible  action  or  feasible  solution  are 
predicted  on  the  basis  of  analysis. 

2.  These  outcomes  are  then  evaluated  according  to  some  scale  of 
value  or  desirability. 

3.  A  criterion  of  decision,  based  on  the  objectives  of  the  system, 

is  used  to  determine  the  most  desirable  action  or  optimal  solution. 

Viewing  the  effort  required  to  produce  a  complete  PADS  from  this  angle 
indicates  that  we  are  starting  with  the  capability  for  Block  1  and  the  basis 
for  Mock  2,  but  must  develop  the  capability  for  generating  multiple 
feasible  solutions,  for  scaling  values  or  desirability,  for  implementation 
cf  decision  criteria,  and  for  selecting  the  optimum  solution. 

The  objectives  of  this  section  are  to  discuss  optimization  briefly 
and  some  working  techniques  available  for  use,  recommend  the  next  phase 
of  effort  to  be  undertaken  to  develop  an  optimization  capacity  for 
VESYS  II  and  to  discuss  some  possibilities  of  subsystem  optimization  to 
increase  the  efficiency  of  the  PADS  as  it  develops. 

7.5.1  Mechanics  of  Optimization 

Decision  criteria  may  be  applied  in  the  form  of  input  constraints 
that  limit  the  generation  of  design  strategies  of  no  interest  or  in  an 
optimization  subroutine  after  the  strategies  are  generated.   The  best 
approach  will  undoubtedly  be  application  of  some  factors  as  constraints 
and  others  in  an  optimization  subroutine. 
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7-10   SYSTEM  ANALYSIS  AND  OPTIMIZATION  (After  [73]). 
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Fig  .  7-11  is  an  FPS  flow  diagram  that  illustrates  how  the  constraints 
are  applied.   This  system  hoth  limits  the  range  of  parameters  analyzed  and 
the  number  of  solutions  processed  further  in  the  system.   Such  constraints 
arc  required  to  limit  consideration  to  suitably  practical  design  strategies. 
While  the  number  of  strategies  generated  is  net  directly  limited  to  some 
practical  number,  limitations  may  be  applied  by  selection  of  magnitude 
of  incremental  layer  thickness  variations  or  as  proposed  previously  for 
subsystem  optimization.   Other  constraints  not  known  in  Fig.  7-11  such  as 
specified  reliability  and  maximum  average  annual  maintenance  costs  can 
be  applied  in  the  same  manner. 

The  bottom  block  of  Fig.  7-11  represents  the  optimization  on  total 
cost  of  the  whole  range  of  feasible  designs  generated  within  the  frame- 
work of  the  constraints. 

Fig.  7-12  portrays  an  optimization  system  used  in  a  Rigid  Pavement 
System  (RPS)  ,  which  could  also  be  applied  to  a  flexible  pavement  system. 
This  system  accumulates  only  a  specified  input  number  of  the  least  expensive 
design  strategies.   Design  strategies  are  generated  up  to  the  number 
specified  without  any  optimization.   The  most  expensive  is  then  identified 
and  subsequent  design  strategies  compared  to  it.   More  expensive  design 
strategies  are  not  stored,  but  less  expensive  strategies  replace  the  current 
most  expensive  in  storage.   New  comparisons  are  made  as  each  design  strategy 
is  processed  to  determine  which  strategy  is  the  most  costly. 

While  both  Fig.  7-11  and  7-12  represent  systems  to  be  optimized  en  the 
basis  of  total  costs,  there  is  no   reason  that  the  same  general  procedures 
cculd  not  be  applied  even  if  the  basis  for  optimization  was  a  number 
representing  a  combination  of  weighted  utilities. 

7.5.2   Subsystem  Optimization  to  Improve  Efficiency  and  Reduce  Computation 
Time  Requirements 

While  not  directly  applicable  to  the  primary  objective  of  optimizing 
the  selection  of  design  strategies,  the  relative  optimization  of  some  of 
the  subsystems  in  the  PADS  would  allow  more  flexibility  in  design  strategy 
generation  and  consideration  within  the  constraints  of  practical  computer 
time  requirements. 

Versions  of  the  Texas  FPS  implemented  were  deterministic  rather  than 
stochastic  and  are  not  as  complex  and  detailed  as  those  subsequently 
developed  and  the  VESYS  program.   There  was  no  constraint  to  solving 
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Figure  7-11   SUMMARY  FLOW  CHART  ILLUSTRATING  MECHANICS 
OF  THE  FPS  PROGRAM  (After  [58]). 
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Figure  7-12   OPTIMIZATION  PROCEDURE  RPS1  (After  [59]) 
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hundreds  of  design  strategies  that  would  maintain  a  minimum  serviceability 
level  and  optimizing  on  the  basis  of  costs.   As  the  subsystems  increase 
in  complexity,  there  is  more  need  for  means  of  eliminating  strategies  of 
little  promise  without  processing  through  the  entire  system. 

There  are  a  number  of  promising  possibilities  that  may  be  considered 
and  implemented  in  later  research  efforts  to  optimize  strategy  generation 
such  that  only  strategies  with  a  reasonable  probability  of  successful 
competition  with  other  strategies  would  be  generated  in  detail  and 
analyzed.   Some  possibilities  follow: 

1.  It  appears  possible  to  run  a  large  factorial  of  representative  and  selected 
solutions  to  develop  a  regression  equation  capable  of  a  reasonably 
accurate  prediction  of  serviceability,  thus  temporarily  bypassing 
most  of  the  time-consuming  portions  of  the  structure  model. 

The  predicted  serviceability  in  combination  with  other  parameters 
could  be  used  as  a  basis  for  elimination  of  the  majority  of 
possible  strategies  prior  to  use  of  the  viscoelastic  structural 
model.  A  similar  analysis  could  be  used  to  develop  regression 
curves  for  overlay  strategies  to  reduce  those  considered  for  each 
of  these  surviving  initial  construction  strategies. 

2.  A  similar  but  somewhat  different  approach  might  be  used  to  develop 
probabilities  of  successful  competition  for  different  strategies 
after  several  viscoelastic  performance  predictions  had  been 
completed.   For  instance,  after  three  solutions  at  varying  thick- 
ness of  a  layer  are  completed,  considerable  insight  may  be  gained  as 
to  the  feasibility  cf  using  other  thicknesses.   This  approach  would 
basically  amount  to  creating  trends  with  a  few  solutions  and 
extrapolating  this  data  to  eliminate  other  variations  of  a  particular 
parameter.   This  might  be  accomplished  for  several  parameters 
simultaneously  or  separately. 

3.  Use  the  AASHO  equations  for  predicting  primary  responses  as  in 
FFS9  to  establish  bounds  on  the  array  before  proceeding  with 
viscoelastic  solutions. 

A.   Use  a  few  viscoelastic  solutions  at  large  increments  of  thickness 
to  "home-in"  on  the  layer  thicknesses  for  initial  construction, 
AASHO  equations  over  the  limited  array  including  overlays  and 
then  visco-elastic  solutions  again  over  the  final  small  array  of 
strategies  for  selection. 
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5.   Use  some  version  of  the  Texas  FPS  system  to  optimize  to  a  small 

array  of  design  strategies  and  then  apply  VESYS  to  that  small  array 
of  strategies  to  optimize  on  the  viscoelastic  basis. 

Another  very  attractive  option  would  be.  to  maintain  both  the  FPS 
and  VESYS  structural  models  in  one  computer  program  and  to  develop  the 
flexibility  to  utilize  either  or  both,  preferably  with  the  capability  of 
using  one  to  reduce  the  array  of  possible  design  strategies  and  the  other 
or  both  to  generate  the  array.   Such  a  computer  program  would  be  a  powerful 
tool  in  the  continuing  refinement  of  the  models  toward  increasingly  accurate 
predictions  of  real-life  responses  and  performance  of  pavement  systems. 
A  similar  multiplicity  could  be  added  for  candidate  optimizational  techniques 
to  allow  their  comparison  for  the  same  set  of  design  strategies. 
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7.6   CONCEPTUAL  PROGRAM 

It  is  quite  clear  that  development  of  full  optimization  of  the  VESYS 
methods  must  necessarily  be  phased,  both  due  to  the  complexity  involved  and 
because  the  means  of  defining  all  the  decision  criteria  identified  do  not 
yet  exibt.  in  such  manner  that  they  may  be  compared  and  utilized  directly. 
This  does  not  mean,  however,  that  useful  and  implementable  systems  may 
not  be  produced  at  each  phase  beginning  now. 

The  purpose  of  this  section  is  to  explore  the  additional  work  which 
appears  to  be  appropriate,  to  develop  reasonable  recommendations  for 
phasing  this  work  and  to  define,  in  some  detail  the  proposed  first  phase 
for  continued  development. 

7.6.1   Definition  of  Effort  Required  to  Achieve  a  Rational  Flexible 
Pavement  Analysis  and  Design  System 

The  required  capabilities  of  the  PADS  are  illustrated  broadly  in 
Fig.  7-2.  It  was  pointed  out  previously  that  VESYS  IIM  essentially  represents 
the  completion  of  Phase  I  and  that  no  optimization  could  occur  until 
Phases  II  and  III  had  been  completed.   Also,  VESYS  IIM  has  no  submodel 
for  considering  the  effects  of  differential  volume  change  in  expansive 
clay  subgrades. 

There  appears  to  be  no  reason  not  to  take  advantage  of  the  state-of- 
the-art  efforts  already  expended  for  the  Texas  FPS  to  expand  VESYS  so  that 
the  VESYS  "black  box"  may  with  minimal  effort  encompass  the  Phase  II  and 
III  capabilities  shown  in  Fig.  7-2. Other  than  the  disparity  in  state  of 
development,  the  principal  difference  between  VESYS  and  FPS  is  in  the 
nature  of  materials  characterizations,  the  method  of  considering  traffic  and 
the  model  for  developing  the  primary  responses  of  rutting,  cracking  and 
roughness.   There  should  be  no  major  constraints  to  replacing  the  FPS 
structural  model  with  the  five-layer  version  of  VESYS  developed  under  this 
project  as  the  first  phase  of  continued  development  effort.   The  resulting 
version  of  VESYS  would  have  the  capabilities  presently  available  in  FPS 
encompassing  the  visco-elastic  structural  model  and  would  provide  a  "platform" 
for  the  next  phase  of  optimization.   Care  should  be  taken  to  maintain 
"modularity"  in  the  ensuing  computer  programs  so  that  subsequent  refinements 
may  be  made  as  easily  as  possible. 
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The  required  effort  for  continued  development  of  a  PADS  is  discussed 
above  with  the  exception  of  implementation.   Implementation  itself  is  a 
fairly  major  undertaking  involving  1)  careful  education  of  those  involved, 
2)  salesmanship,  3)  interfacing  with  existing  data  collection  and  design 
systems  and  4)  establishment  of  a  data  feedback  loop  and  means  of  monitoring 
it  and  refining  the  models.   Also  included  with  the  major  tasks  will 
be  a  number  of  related  tasks  such  as  these: 

1.  Further  development  of  information  on  important  design  parameters. 
The  relative  importance  of  the  parameters  will  have  for  the  most 
part  been  established  through  sensitivity  studies  such  as  those 
reported  in  [66]   and  those  for  VESYS  IIM  in  progress  by  Austin 
Research  Engineers,  Inc  under  FHWA  Contract  DOT-FH-11-8258 .   Those 
with  the  greatest  impact  on  design  selections  should  be  emphasized 
in  research. 

2.  Optimization  of  the  subsystem  to  improve  efficiency  and  reduce 
computational  time  requirements  as  previously  discussed. 

3.  Comparative  studies  of  the  results  obtained  from  VESYS  and  FPS 
prediction  models  for  rut  depth,  fatigue  and  roughness  to  identify 
causes  of  any  variations.   This  should  lead  to  possibilities  for 
refinement  of  the  models. 

4.  Continued  development  of  stochastic  data  on  design  parameters 
for  use  in  design  on  a  probablistic  basis.   Primary  emphasis 
should  be  given  those  parameters  that  have  been  found  to  have 
strong  importance  to  design  results. 

7.6.2   Phasing  of  Continued  PADS  Development 

The  first  phase  recommended  is  the  "marriage"  of  the  five-layer 
VESYS  system  with  the  most  advanced  version  of  the  Texas  FPS,  Computer 
Program  FPS-13.   Appropriate  features  from  SAMP6  such  as  the  cost  sub- 
routines for  consideration  of  shoulders  and  cross-section  volumes  may  also 
be  included.   The  resulting  computer  program  would  be  a  quite  advanced 
system  capable  of  developing  full-life  design  strategies  to  particular 
levels  of  desired  reliability.   Such  a  program  could  serve  both  as  an 
in piemen  table  system  and  a  platform  for  continued  research. 

Phase  II  should  include  the  following: 

1.   Improvement  of  the  computer  program  developed  in  Phase  I  to 
increase  its  efficiency  through  application  of  some  or  all  of 
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the  concepts  described  previously  as  subsystem  optimization. 

2.  Comparative  studies  of  results  obtained  from  the  VESYS  and  FPS 
prediction  models  for  rut  depth,  fatigue  and  roughness.   Refine- 
ment of  models  as  indicated  by  this  or  other  studies. 

3.  Studies  aimed  at  rationally  including  seal  coat  or  other  surface 
repair  strategies  for  maintaining  acceptable  levels  of  skid 
resistance. 

4.  Other  studies  or  developments  identified  as  requirements  during 
Phase  I. 

Phase  III  could  involve  several  separate  efforts  to  include: 

1.  Implementation  in  cooperating  state  highway  departments. 

2.  Development  of  a  basis  for  consideration  of  other  decision  criteria 
and  for  assigning  relative  weights  to  reflect  their  importance. 

The  expanded  optimization  capability  should  also  be  added  to  the 
latest  version  of  the  computer  program  and  experiments  run  to 
study  its  effects. 

3.  Study  of  other  distress  mechanisms  and  limiting  responses  to 
continue  improvement  of  the  predictive  model.   Refine  the  predic- 
tive model  to  reflect  the  results  of  these  studies. 

A.   Further  study  of  the  stochastic  characteristics  of  the  important 
design  parameters. 

Future  work  past  Phase  III  toward  realization  of  an  accurate  and 
reliable  PADS  will  generally  be  defined  as  work  progresses,  but  it  should 
include  improvement  of  the  implemented  versions  as  new  refinements  warrant. 

7.6.3  Concept  for  Phase  I  PADS 

A  summary  flow  diagram  is  shown  in  Fig.  7-13  for  the  proposed  combination 
of  computer  program  FPS-13  with  the  five-layer  VESYS  computer  program 
developed  under  this  project.   While  this  flow  diagram  indicates  fairly 
clearly  the  nature  of  the  procedure,  it  may  be  useful  to  follow  through 
the  various  steps  and  discuss  them  separately. 

As  indicated  in  the  upper  right-hand  corner  of   Fig.  7-13,  a  series  of 
initial  designs  may  be  considered  including  variations  in  thicknesses  and 
properties  of  the  asphalt  concrete,  base  and  subbase  (if  any)  and  in 
properties  for  the  subgrade  (i.e.,  treated  and  untreated).   The  specified 
ccmbinations  of  layers  may  be  considered  iteratively. 
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The  model  fot  volume  change  and  resulting  roughness  due  to  clay 
subgrades  could  also  be  included  in  VESYS  at  the  same  time  to  refine  the 
performance  predictions  for  this  fairly  strong  parameter. 

The  objective  of  this  recommendation  is  to  produce  a  PADS  as  quickly 
as  possible  so  that  the  business  of  refinement  and  development  may  proceed. 
The  cost  base  procedure  in  use  in  the  Texas  FPS  and  proposed  above  for 
VESYS  provides  a  reasonable  selection  procedure  pending  development  of 
more  sophisticated  benefit  or  attribute  analysis  and  the  result  of  this 
"marriage"  should  represent  an  implementable  PADS. 

The  structural  model  itself  needs  to  be  subjected  to  continued  develop- 
ment and  scrutiny  to  include  other  distress  mechanisms  and  refinement  of 
those  presently  considered  as  new  knowledge  is  developed  from  other  studies. 
One  facet  of  this  effort  should  be  continuing  long-term  comparisons  between 
predicted  and  actual  performance  within  a  "data  feedback"  system  such  as 
that  described  in   [68].   Study  of  causes  of  disparities  could  lead  to 
identification  of  deficiencies  in  the  models  that  could  then  be  corrected. 
These  efforts  may  probably  best  be  made  cooperatively  with  states  where  the 
PADS  is  implemented. 

The  optimization  for  the  Texas  FPS  system  is  currently  limited  to 
indirect  imposition  of  several  of  the  desired  decision  criteria  as  constraints 
and  arraying  the  best  design  strategies  on  the  basis  of  costs  so  that  a 
designer  may  apply  this  and  other  decision  criteria  of  his  choice  in  making 
a  selection.   As  indicated  previously,  it  is  conceptually  intended  that 
the  rational  PADS  will  include  a  broader  optimization  involving  weighted 
decision  criteria  other  than  the  costs  now  considered.   Such  optimization 
would  involve  some  means  of  assigning  weights  to  the  various  social,  political 
and  cost  criteria  in  terms  of  their  relative  benefits  or  utilities. 

The  constraint  on  implementation  of  such  an  expanded  optimization 
system  is  net  in  the  mechanics,  but  in  the  identification  of  all  the 
criteria  that  should  be  included  and  the  relative  importance  to  attach 
to  each.   The  identification  and  quantification  of  these  criteria  represents 
a  research  phase  of  its  own,  which  could  be  carried  on  separately  and 
concurrently  with  other  development  of  the  PADS.   The  approach  here  should  be 
one  of  continued  refinement  of  a  cost  optimization  system  that  may  be 
reasonably  and  effectively  used  in  the  interim  as  it  represents  a  much- 
improved  extension  of  the  cost-based  design  decision  well-established  in  the 
profession. 
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For  a  particular  initial  design  i,  the  total' initial  thickness  is 
compared  to  a  specified  maximum  value.   If  too  thick,  this  design  is 
discarded  and  the  next  lower  thickness  tried.   If  not  too  thick,  design  i 
enters  the  five-layer  VESYS  structural  and  serviceability  prediction  model. 

The  five-layer  VESYS,  which  is  illustrated  in  more  detail  in  Fig  7 
and  previously  described,  will  predict  periods  of  satisfactory  service- 
ability at  different  levels  of  reliability,  but  it  is  necessary  to  select 
one  period  for  use.   The  most  straightforward  method  of  selection  is  to 
establish  a  required  level  of  reliability  and  interpolate  for  a  period  of  time  to 
reach  the  minimum  serviceability  level  between  pairs  of  time  to  minimum  service- 
ability level  and  reliability  values  calculated  by  the  five-layer  VESYS  subsystem. 

The  initial  life  (period  for  serviceability  to  decay  to  the  minimum 
allowable  level)  would  then  be  compared  to  a  specified  minimum  time  to  first 
overlay.   If  the  initial  life  is  less  than  the  minimum,  design  i  is  discarded 
and  the  next  stronger  design  is  considered.   If  equal  or  greater,  the 
initial  construction  cost  is  calculated  and  compared  to  the  specified 
maximum  initial  funds  available. 

If  the  estimated  initial  cost  exceeds  the  maximum  initial  funds 
available,  design  i  is  discarded  and  a  new  initial  design  is  considered. 
If  the  estimated  initial  cost  is  within  the  initial  funds  available,  a  check 
is  made  to  see  if  the  initial  life  predicted  equals  or  exceeds  the  analysis 
period.   If  so,  no  overlay  is  required  and  the  initial  construction  cost 
plus  annual  maintenance  costs  reduced  to  present  worth  represent  the  total 
overall  cost  for  the  analysis  period. 

If  the  initial  life  is  less  than  the  analysis  period,  various  overlay 
strategies   must  be  considered,  also  on  an  iterative  basis  for  increasing 
overlay  thicknesses  and  material  properties.   The  present  serviceability 
index  after  an  overlay  must  be  estimated  and  specified  as  for  the  initial 
pavement. 

The  "life"  or  period  of  time  before  the  serviceability  of  the  overlayed 
pavement  structure  deteriorates  to  the  established  minimum  level  may  be 
predicted  using  the  five-layer  VESYS  subsystem  as  for  the  initial  design, 
but  it  will  be  necessary  to  decide  how  to  model  the  various  layers.   For 
instance,  an  initial  design  having  a  subgrade,  a  treated  subgrade  layer,  a 
subbase  layer,  a  base  layer  and  the  asphaltic  concrete  layer  does  not  leave 
another  available  layer  in  the  "model"  for  an  overlay.  Also,  a  second  or 
third  overlay  may  be  required.   It  is  possible  to  consider  both  the  original 
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asphaltic  concrete  and  subsequent  overlays  as  one  thicker  layer  if  their 
material  properties  are  similar.   Similarly,  the  subbase  and  base  may  be 
combined  if  appropriate.   Another  decision  is  required  whether  the  initial 
material  properties  will  be  considered  to  apply  or  should  they  be  revised 
to  reflect  the  anticipated  effects  from  repetitive  loading  and  environment 
prior  to  overlay.   These  problems  are  subject  to  solution  with  a  reasonable 
amount  of  effort,  but  must  be  studies  during  the  Phase  I  development. 

After  the  life  for  the  overlay  is  computed,  it  would  be  compared  to 
the  specified  minimum  period  between  overlays.   If  less  than  the  minimum, 
the  overlay  strategy  is  abandoned  and  another  thickness  or  material  tried. 
If  greater  than  the  minimum,  the  combined  lives  for  the  initial  construc- 
tion and  the  overlayed  pavement  system  are  compared  to  the  analysis  period. 
If  the  combined  total  life  is  greater  than  the  analysis  period,  design  i 
and  the  one  overlay  enter  the  subsystem  for  calculation  of  total  overall 
costs.   If  lesser,  another  overlay  is  required  and  the  procedure  is  repeated. 

When  sufficient  overlays  have  been  included  to  extend  the  life  through 
the  analysis  period,  the  total  overall  cost  for  this  design  strategy  is 
computed.   The  total  overall  cost  will  include  the  initial  construction 
cost,  overlay  costs,  annual  maintenance  costs  and  traffic  delays  less  the 
salvage  value.   All  future  costs  are  reduced  to  present  worth  for  comparison, 

The  total  overall  cost  is  compared  to  those  for  other  design  strategies 
as  illustrated  in  Fig.  7-12.  The  new  strategy  will  be  retained  in  storage  if 
fewer  than  the  number  of  design  strategies  specified  for  printing  out  are 
stored  or  if  it  is  cheaper  than  some  strategy  in  storage  that  it  may  replace, 

After  the  comparison  with  other  design  strategies,  a  new  overlay 
iteration  begins  with  a  new  overlay  thickness  or  material.  After  all 
specified  overlay  strategies  have  been  considered  and  compared,  a  new 
iteration  of  the  entire  procedure  begins  with  a  new  initial  design.  After 
the  array  of  initial  designs  and  overlay  strategies  for  each  have  been 
processed,  the  surviving  overall  design  strategies  are  printed  out  and 
the  designer  may  study  them  and  make  the  final  design  selection. 

It  would  be  helpful  for  the  reader  to  read  over  the  data  input  guide 
for  FPS-13   [66]   or  SAHP6   [62]    to  become  familar  with  the  input 
data  required  and  its  organization.   The  majority  of  the  data  could  be 
input  without  revision,  but  some  revisions  will  be  required  including  the 
following: 
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1.  FPS-13  assumes  an  analysis  period  of  twenty  years  and  flexibilty 
for  consideration  of  other  specified  periods  if  desired,  would 
require  revision. 

2.  FPS-13  operates  with  18-k.ip  axle  equivalencies,  while  VESYS  would 
require  specific  tire  pressures  and  radii  of  the  loaded  areas. 

3.  VESYS  required  stochastic  creep  compliance  data  while  FPS-13 
uses  strength  coefficients  for  characterization  of  the  layer 
materials . 

4.  VESYS  requires  a  vector  of  temperature  values  that  are  used  to 
describe  an  annual  cycle,  a  reference  temperature  for  the  master 
creep  compliance  curves,  a  time-temperature  shift  factor  and 
other  inputs  for  considering  the  effects  of  temperature  on 
material  characteristics,  whereas  FPS-13  uses  a  stochastic 
district  or  reginal  temperature  constant  with  internal 
subroutines . 

5.  Traffic  data  is  required  by  VESYS  in  vehicles  per  day  (usually 
only  trucks  considered)  and  a  stochastic  repeated  load  intensity 
and  duration  while  FPS-13  uses  average  daily  traffic  and  percent 
trucks  en  a  stochastic  basis. 

6.  Other  VESYS  input  requirements  would  have  to  be  provided  for. 

The  expansion  to  include  the  shoulder  and  cross-section  volume  costs  in 
optimization  will  involve  integration  of  these  cost  subroutines  from  the 
SAMP6  computer  program  into  the  cost  subroutines  for  FPS-13.   This  should 
not  be  difficult  due  to  the  similarity  of  the  two  computer  programs.   Detail 
on  these  SAKP6  subroutines  and  input  requirements  appears  in  Chapter  9  of 
Appendix  C,  of  [ 62 ] . 

It  would  be  desirable  to  develop  a  research  version  of  the  combined 
Phase  I  computer  program  as  discussed  previously  to  include  the  structural 
and  predictive  models  for  both  the  VESYS  and  FPS  systems  with  as  few 
variations  in  input  as  necessary.   This  would  allow  relatively  easy  comparative 
studies  that  might  lead  to  refinements  of  the  prediction  models  for  rut  depth, 
fatigue  and  roughness  or  other  improvements.   The  effort  involved  in  develop- 
ment of  such  a  research  program  concurrently  with  the  development  of  the 
the  implcmentable  VESYS/FPS-13  combination  should  be  rather  nominal.   If 
this  research  version  of  the  combined  Phase  I  program  was  not  produced, 
comparative  studies  recommended  in  Phase  II  could  still  be  conducted  using 
the  separate  computer  programs. 
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7 .  7   SUTMARY 

The  very  considerable  body  of  previous  work  toward  development  of 
rational  analysis  and  design  systems  for  flexible  pavements  has  been 
discussed  and  synthesized.   The  structures  of  the  systems  and  inter- 
relationships of  the  various  system  components  were  displayed  in  diagram 
form  and  the  flow  of  system  operations  discussed. 

Specific  decision  criteria  for  use  in  optimization  of  design  were 
defined  and  discussed  in  detail.   While  it  is  possible  to  define  the 
majority  of  the  decision  criteria,  no  feasible  basis  for  assigning 
utilities  and  comparison  has  been  developed  for  most  of  them.   This  has 
led  to  the  use  of  total  cost  as  a  basis  for  optimization.   Total  cost 
as  currently  used  includes  initial  construction  costs,  routine  mainten- 
ance costs,  costs  for  overlays  and  seal  costs,  user  costs  due  to  delays 
during  overlay  construction  and  the  negative  cost  of  salvage  value  after 
the  design  period  all  considered  at  their  "present  worth". 

Techniques  for  design  optimization  were  discussed,  but  there  can  be 
no  optimization  for  the  VESYS  system  until  capabilities  for  considering 
multiple  initial  design  strategies  and  generating  maintenance  policies 
has  been  incorporated.   A  "marriage"  of  the  new  five-layer  VESYS  computer 
program  is  proposed  as  the  next  phase  of  research  to  provide  a  working 
VESYS  Pavement  Analysis  and  Design  System.   Other  research  and  develop- 
ment needs  are  defined  and  their  phasing  proposed. 

There  appears  to  be  no  serious  constraint  on  the  development  of  an 
imp lemen table  PADS  during  the  next  research  phase  proposed. 
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CHAPTER  8 
MODULAR  FRAMEWORK  OF  VESYS  IIM 

Modifications  made  to  VEYSY  II  to  develop  VSYSY  IIM  have  been 
made  in  such  a  fashion  as  to  facilitate  future  changes  and  in  corpor- 
ation of  new  subsystems.  The  VESYS  IIM  computer  routines  have  not 
been  completely  coupled  to  VESYS  II,  however.  Also,  a  significant 
protion  of  VESYS  II  is  not  modular.  Coupling  our  portions  of  VESYS 
IIM  with  the  remaining  VESYS  II  distress,  serviceability  and  maintain- 
ability subsystems  in  a  modular  fashion  is  not  difficult,  however, 

The  various  computer  programs  in  VESYS  IIM  have  been  assembled 

in  a  modular  fashion  into  three  program  elements: 

READ  PROGRAM 
PRIMARY  RESPONSE 
DISTRESS  ANALYSIS 

Rather  than  rely  on  machine  dependent  system  commands  to  control 

execution  of  the  various  program  elements,  a  simple  executive  program 

written  if  Fortran  languate  is  used  to  execute  the  various  subsystems. 

This  procedure,  while  not  as  efficient  for  a  specific  computer, 

greatly  simplifies  the  task  of  transferring  the  program  from  one 

machine  to  another.  Figure  8-1  presents  a  shcematic  of  the  organization 

of  VESYS  IIM. 
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The  executive  program  controls  executions  of  all  program  elements. 
The  READ,  PRIMARY  RESPONSE,  AND  DISTRESS  ANALYSIS  programs  are 
sotred  on  separately  accessible  files.  Each  subsystem  of  these 
programs  represents  a  subroutine  which  may  contain  further  subelements 
depending  upon  the  calculations  being  made. 

The  EXECUTIVE  program  calls  the  READ  program  and  tells  it  what 
input  data  to  read.  This  data  is  stored  on  a  file  or  on  a  tape  and/or 
printed.  Control  is  returned  to  the  EXECUTIVE  program  which  may  call 
for  an  exit  or  execution  of  any  of  the  subsystems  of  the  PRIMARY 
RESPONSE  OR  DISTRESS  ANALYSIS  systems.  If  the  PRIMARY  RESPONSE 
system  is  called,  the  reguired  input  data  is  read  from  FILE  1. 
The  output  of  PRIMARY  RESPONSE  is  stored  in  FILE  2   and/or  printed, 
and  control  is  returned  to  EXECUTIVE.  If  the  DISTRESS  ANALYSIS  system 
is  called,  the  input  data  required  is  read  from  FILE  1  and  FILE  2  and 
the  output  is  stored  in  FILE  3  or  printed. 

With  the  organization  outlined  above,  it  is  possible  to  execute 
certain  portions  of  the  program  separately  as  well  as  make  modifi- 
cations to  any  system  in  a  straight  forward  manner  as  long  as  consistent 
input-output  formats  are  maintained,  without  having  to  change  significant 
numbers  of  common  blocks. 
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CHAPTER  9 
SUMMARY  AND  RECOMMENDATIONS 

The  significant  accomplishments  during  this  contract  consist  of: 

1.  Extension  three-layered  viscoelastic  solution 
to  an  N-layered  viscoelastic  solution. 

2.  Development  of  an  approximate  flosed  form  probabilistic 
solution  for  stresses,  strains  and  deflections  for 

an  N-layered  viscoelastic  system. 

3.  Development  of  an  improved  rut  depth  model  which  uses 
simplified  laboratory  characterization  tests. 

Additional  work  is  required,  however,  to  complete  development 
of  a  rational  pavement  design  framework.  First,  VESYS  II  and  VESYS  I IM 
should  be  merged  completely  and  a  series  of  sensitivity  analyses 
conducted  on  the  resulting  design  system. 

An  approximate  probabilistic  model  developed  for  the  primary 
response  node!  should  be  validated  and  the  limits  of  applicability 
established.  If  the  range  of  applicability  is  not  as  large  as  desired, 
the  first  four  moments,  rather  than  the  first  two  used  in  VESYS  IIM, 
may  be  necessary  to  determine  probability  distribution  functions. 

The  primary  response  model  should  also  be  extended  to  include 
multiple  wheel  loads  with  variable  contact  areas  and  variable  spacing. 
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The  extension  is  straight  forward  assuming  linear  behavior;  however, 
it  may  be  worth  while  to  investigate  a  finite  element  formulation 
to  allow  even  greater  flexibility  in  the  loadinq  and  material 
behavior  characteristics.  Particularly,  if  the  non-linear  behavior 
of  the  base  and  subgrade  materials  is  to  be  adequately  accounted  for 
in  evaluating  pavement  performance. 

The  rut  depth  model  developed  under  this  contract  should  be 
validated  with  field  studies.  Additional  validation  of  other 
distress  systems  in  VESYS  II  is  also  required,  including  incorporation 
of  viscoelastic  fracture  mechanics  concepts  rather  that  the  elastic 
approaches  currently  used. 

Also, due  to  the  importance  of  pavement  roughness  in  evaluating 
pavement  performance,  a  mechanistic  roughness  model  is  needed.  A 
simple  approach  is  this  problem  is  to  model  the  interaction  between 
vehicle  and  pavement  as  a  simple  two  degree-of-freedom,  spring-mass 
system  moving  on  a  Maxwell  model  viscoelastic  pavement.  Roughness 
may  then  be  characterized  by  the  root  mean  square  value  of  the 
pavement  profile  assuming  that  the  pavement-wheel  interaction  is  the 
main  contributor  to  the  change  in  pavement  roughness.  The  change  in 
the  random  pavement  profile  can  be  investigated  with  spectral  and 
correlation  analyses. 

Finally,  simplified  material  characterization  procedures  and 
automated  data  reduction  routines  are  needed  to  increase  the  appeal 
and  usefulness  of  a  rational  pavement  design  framework  to  highway 
designers.  A  simplified  testing  procedure  is  incorporated  in  the 
permanent  defermation  model  developed;  however,  data  reduction 
may  be  more  complicated  requiring  more  training  of  laboratory 
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technicians.  Accordingly,  it  is  recomnended  that  a  material  charac- 
terization subsystem  based  on  numerical  optimization  techniques  be 
incorporated  into  VESYS  I IM  to  perform  the  necessary  data  reduction 
for  materials'  response  and  damage  behavior. 
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APPENDIX  A 


PRIMARY  RESPONSE  MODEL 
PROGRAM  DOCUMENTATION 


I.   SUBROUTINES  AND  VARIABLES 
MAIN 

MAIN  is  the  executinve  control  routine.  All  subroutines  are 
called  from  this  program  and  all  input  data  and  output  are  printed. 

Several  variables  pertaining  to  control  of  data  executions 
are  defined  in  this  program: 

ITN  =  number  of  terms  used  in  numerical  integration 
(default  value  =  46) 
ITN4  =  ITN  x  4  (internally  calculated) 
REED 

All  data  are  read-in  in  this  subprogram.  The  input  is  then 
passed  to  MAIN. 

TITLE  (I)  =  job  title 

NS  =  number  of  layers 
N  =  NS-1  (calculated  internally) 

IR  =  number  of  radial  points  where  output  is  desired 
IT  =  number  of  axial  points  where  output  is  desired 
WGT  =  total  wight  of  load 
AR  =  load  radius 
PSI  =  tire  pressure 
V(I)  =  Poisson's  ratio  for  i-th  layer  elastic 
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E(I)  =  modulus  or  compliance  for  i-th  layer 
PR(I)  =  radial  coordinate 
AA(I)  =  axial  coordinate 
COVE(I)  =  coefficient  of  variation  of  elastic  modulus 
or  compliance  of  i-th  layer 
NT ( I )  =  number  of  time  for  which  lab  creep  data  is 
curve  fit  for  i-th  layer 
TYME(I,J)  =  vector  of  times  for  lab  creep  data 
VALVE(J,I)  =  creep,  compliance  at  time  J  for  layer  I 
C0V(J,I)  =  coefficient  of  variation  of  creep  compliance 
or  modulus  at  time  J  for  layer  I 
NTO  =  number  of  times  for  which  output  is  desired 
T1(I)  =  vector  of  times  for  which  output  is  desired 
COF(K)  =  coefficient  in  power  law  curve-fit  of  compliance 
DN(K)  =  exponent  in  power  law  curve-fit  of  compliance 
A  flag  is  read  to  determine  if  an  elastic  or  viscoelastic  analysis 
is  to  be  performed  and  if  laboratory  data  are  to  be  curve-fit. 
Either  modulus  or  compliance  may  be  input.  The  program  compares  the 
magnitude  with  unity  to  determine  which  has  been  input.  If  the  coef- 
ficients of  variation  input  are  zero  for  all  layers  then  a  deterministic 
solution  is  generated. 

If  compliances  are  input,  then  the  mean  and  variance  of  moduli 
used  internally  in  the  program  are  calculated  from* 


E[Ei]  - 


E. 
l 


*  Our  thanks  are  due  Dr.  R.L.  Lytton  of  Texas  A 
expressions. 


&  M  University  for  these 
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2  r 


VarCE^  = 


E.1* 
l 


PART 


th 


This  subroutine  performs  the  numberical  partition  of  {X.}  and 
e  four  points  in  each  interval  (X.  ,,X.)  at  which  the  integration 


of  /  T(m)drn  is  evaluated.  The  set  of  X.'s  are  so  chosen  from  the 

zeros  of  T(m)  such  that  in  the  interior  of  any  interval  (X.  -,,X.) 

there  is  at  most,  one  zero  of  T(m). 

The  variables  in  this  subroutine  are: 

BZ(I)  =  the  roots  of  the  Bessel  functions  J  and  J. 

o     1 

GT,G2  =  the  roots  of  the  Legendre  polynomial 

AZ(K)  =  the  parameter  m  required  in  the  evaluation  of 

/o°°  T(m)dm 

The  value  of  ITN,  which  represents  the  number  of  terms  used  in 

the  series  approximation  to  the  integral  is  set  at  46.  This  value 

runs  on  the  UNIVAC  1108  Computer,  however,  the  size  of  a  double 

precision  word  is  less  on  IBM.  Therefore,  the  maximum  number  of  terms 

may  have  to  be  decreased. 

CQEF 

3A.   8B.   8C. 
The  constants  A.,B.,C.  and  [).,  and  -^-   ,  t4  ,  -^-  and 

~n  Ill        1        at  •    at  •    dt  • 

oU.  J      J      J 

tt  are  computed  in  this  subroutine.  Most  of  CPU  time  spent  is 
due  to  the  evaluation  of  these  constants: 
The  variables  defined  here  include: 


PI 


n-1  D'^hJX.Dfh.) 
5     4(JVi) 


202 


DPI 


-1  3XJ 


DEN 


the  determinant  of  the  2x2  matrix  for 


solving  for  C     and  D 


A(M),   B(M),   C(M)  and  D(M) 

=  constants  A.,  B . ,  C.  and  D . ,  respectively, 


CN 
DN 


where  M=l,2, ,ITN4 

=  C. 


n 


=  D 


n 


DPX,  DPXX,  DPG 

=  intermediate  values  used  in  computing  DPI 

DGA(M,LL),DGB(M,LL),  DGC(M,LL)  and  DGD(M,LL) 

3A.   8Bi   3C.      3D. 
=  "5F~  >  te~  »  tf^"  and  «^"  >  respectively,  where 
3Ej   9Ej   3Ej      3Ej 

i=M,  M=l,2, ITN4,  j=LL,  LL=  layer  index, 

representing  the  layer  to  be  differentiated 


TERMS 


Before  the  constants  A.,B. ,•■•,  etc.,  can  be  evaluated, 


certain  elements  used  in  obtaining  these  constants  have  to  be 

defined. 

This  subroutine  defines  the  following  variables: 
F1(G,Y1,Y2)  =  p(h,v,u)  =  2(mh)2+4mh(v-u)+l-8vu+2y 
F2(G,Y1,Y2)  =  q(h,v,u)  =  2mh(2v-l )-(l+8uv-6u) 
X1(I,J,K)    =  part  of  X(I,J,K) 
X2(I,J,K)    =  part  of  X(I,J,K) 
X(I,J,K)     =  XI (I,J,K)  +  EV  X2(I,J,K) 


1+v 


EV 


i+1 


Ei+1  1+vi 
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DMULT 

This  subroutine  performs  matrix  multiplications  of 

1  gv 

D"  (hj)  Xj  D(hj)  or  D"  (hj)  gIiD(hj)  which  is  required  for  the 
evaluation  of  constants  A(M),  B(M),  C(M),  D(M),  DGA(M.LL),  DGB(M,LL), 
DGC(M,LL)  and  DGD(M,LL). 


BESSEL 

This  subprogram  computes  either  JQ(m)  or  J^m)  to  within 
an  absolute  tolerance  of  less  than  1.1  x  10"7.  The  computation  was 
the  following  approximations 

(1 )  Function  JQU) 

For  X  >  3,  Let  Y  =  3/X 


F(x)  =  .79788456  -  .0OCC0077Y  -  .0055274Y2  -  . 0000951 2Y3 
+  .00137237Y4  -  .00072805Y5  +  .00014476Y6 

T(x)  =  X-.78539816-.04166397Y-.00003954Y2+.00262573Y3 
-  . 000541 25Y4-.0029333Y5+. 0001 3558Y6 


J0(x)  =  F(x)  T(x)/^" 


For  X  <  3,  Let  Y  =  X/3 


JQ(x)  =  l.-2.2499997Y2+1.2656208Y4-.3163866Y6+.0444479Y8 
-  .0039444Y10  +  . 00021 Y12 


(2)  Function  J, (x) 


For  X  <  3,  Let  Y  =  3/X 


F(x)  =  . 79788456+. 00000156Y+. 01659667Y2+. 0001  7105Y3 
-  . 0024951 1Y4  +  .00113653Y5  -  .00020033Y6 
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T(x)  =  X  -  2.35619449+.12499612Y+.0000565Y2-.00637879Y3 
+  .00074343Y4  +  .00079824Y5  -  .00029166Y6 


J^x)  =  F(x)  T(x)/A 
For  X  <  3,  Let  Y  =  X/3 


J^x)  =  .5-.56249985Y2+.21093573Y4-.03954289Y6 

+  . 0044331 9Y8  -  .00031 761 Y10  +  . 00001 109Y12 

CALCIN 

This  subroutine  performs  all  deterministic  calculations. 
CALCN2 

This  subroutine  performs  probabilistic  calculations. 
PRORL 

This  subroutine  is  called  by  CALCN2  and  evaluates  the  sums 


2 


and 


n  /3g.  \2 


CURVIT 

A  least  squares  curve-fit  of  laboratory  creep  compliance  or 
relaxation  modulus  data   is  obtained  for  the  power  law  representation 
J(t)   =  A  tn 
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where 


A  =  exo 


n 

Z  1  on J  . 
i=l 

n 

Z    lont. 
1=1 

n 

Z  lonJ. 
i=l 

logti 

n 

Z  (lont.)2 
i  =  l     n 

DET 


and 


n  = 


DET  = 


N 

Z  logt 
i=l 


N 

Z  lont. 
1=1     n 


N 

Z   loqJ. 
i=l 

n 

Z  logJ.lont 
i=l     ' 


DET 


M 

Z  lont. 
i  =  l     n 

N  (lont.)2 
Z 
i=l 


II.   EVALUATION  OF  EQUATIONS 

Evaluations  of  terns  in  the  probabilistic  solution  is  described 
to  provide  familiarity  with  the  structure  of  the  computer  program. 

In  solving 


3A. 

9B. 

3C. 

3D. 

i 
3F- 

,   l 

3E, 

l 

3E. 

and  3E. 

J 

J 

J 

J 
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where 


-»,- 

i 

m                 « 

3Ej 

0 

9B-i 

I 

»     «^ 

0 

*i 

3P, 

J 

_ 

1 

8C, 

3Ej 

*i 

u     — 

Cn 

9D1 

A 

Dn 

rl 

ap 


*['■] 


o 

3Cn 
D 


^TTvTT 


3EJ 


3E. 


rnj  v  xai 


(D 


(2) 


(3) 


In  the  evaluation  of  Equation  (  3  )  nine  cases  must  be 
considered.  The  cases  form  a  tree  structure: 


NS  =  No.  of  layers 
N  =  NS-1 


207 


Case  1         0  =  j  >  N 

This  case  does  not  exist, 


Case  2         1   =  j  <  N 


3P,       8P 


11 


9E,       3E- 


-1    9X1 

ui    aE1  ur 


4(v,-l) 


'1 


N 

rr 

c=2 


4U.-1) 


-1    3X1 
rUl      3L1    Ul 


4(vrl) 


N 

7T 

£=2 

D"]X( 

>v 

4{Vl. 

■1) 

(4) 


where 


Ek      ]  +  Vl 


Ek+1       ]   +  \ 


k  =  1,2,3, ,n-l 


Case  3        2  <  j   <  N 

In  this  case,   the  possible  conditions  are: 


=  2 

N  =  3 

■  2,3 

N  =  4 

=  2,3,4 

N  =  5 

j  =  2,3,4,--,n-2         N  =  n-1 
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Hi  =  r#  °;Vo  rp-1  ^feii ,  fD.HxJ-ipMi 


[ 

U-l   4(vJL-l)J[hj      [ 


-1 3xj 

4(vrl) 


j-i  aLj.^j-i 


4(v..rl) 


3X 


•1 


•1 


°I1A\  ♦  l.iP.i-ixj-idj-i 

4(v.-l)j       EJ  I  4(^,-1)     . 


•I       OA  . 

L  4(Vj-l)     Jj  [i-3^   4(vrl)J 


(5) 


where 


3E, 


9Xi   3Li 


3X.  rU 


3Li    3^  3L 


J" [iT  «1k   "   E^7  5ak]a=1+] 


Case  4        j  <  1 

This  case  does  not  exist. 


(6) 


Case  5        j  =  N  =  1 


3P 
3E 


1    -     II 
7    "     El 


-1   3X1 
D1      3L7D1 

L   4(V1) 


(7) 
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Case  6         1   <  j  »  N 


9P 

9E 


,-i    !*fcl 


i .  r^2  D^Aif[Dj-'  aV  DJ-f|rpj1xJDjirDj'1-i''j-ipj-i 

j  [l-l   4(V1)J[[      4(vj.rl)      ][4(vrl)JL    4(vj.rl) 


-1   3Xi 


4(vr!)     Jj 


i_  n^2  pl\pi 

Ej     4=1     4(v  -1) 


3Xj-l 
i-1 


Bj.fi 


^'Vj 


-1   3Xj 


,4(v,-l) 


(8) 


Case  7        j  <  1 

This  case  does  not  exist. 

Case  8    j  =  NS  =  1 

This  case  exists  for  N  =  0  (semi -infinite  space)  only. 


Case  9        j  =  NS  >  1 


3P. 


U-l       4(v0-l)JL 


-i  5h  D 

4(v..rl) 


] 


ENSU=1     4(v£-1)J[  4(vN-l)    J 


4(Vi)JJ_  4(-v1) 

N0t6:       iSTl   "^V^    =     1      for  NS     =     2 
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(9) 


III.   USER'S  MANUAL 

The  input  requirements  for  the  primary  response  model  are  defined 
in  Table  I  below: 

TABLE  I 

INPUT  REQUIREMENTS 
PRIMARY  RESPONSE  MODEL 

CARD  PARAMETER 

1  ITITLE  (I)  (20A4)  -  Title  and  Job  information; 

columns  5-80  are  printed  out  with 
job  results. 

IFLG,  NS,  IR,  IZ,  WGT,  AR,  PSI  (415,  3E10.5) 

0  ->  Elastic  Analysis. 

IFLG  =  {  1  •*■  Viscoelastic  Analysis  with  lab  creep 
data  to  be  curve-fit  to  power  law 
A  tn. 

IFLG  /  or  1  ^  Viscoelastic  Analysis;  but  coefficient 
"A"  and  exponent  "n"  of  power  law 
Atn  are  to  be  input. 

NS  =  Number  of  Layers. 

IR  =  Number  of  radial  points  where  out- 
put is  desired. 

IZ  =  Number  of  axial  points  where  out- 
put is  desired. 


WGT  =  Axle  Load    ,  Qnly  ^  Qf   the$e 
AR  =  Load  Radius   }  three  may  be  input; 


PSI  =  Tire  Pressure 


however,  any  two  will 
do. 


3  HH  (I)  (I=1,NS-1)(8E10.5)  -  Layer  thicknesses  for  first  n-1  layers- 

4  RR  (I)  (1  =  1, IR)  (8E10.5)  -  Radial  points  where  output  is  desired. 

5  ZZ  (I)  (1=1, IZ)  (8E10.5)  -  Axial  points  where  output  is  desired. 

IFLG  =  0:   • 

6  E(I),C  E(I),  V(I),  1=1,  NS  (8E10.5) 

E(I)   =  Elastic  Modulus   (or  compliance  of  the 
Ith  layer. 
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COVE(I)  =  Coefficient  of  variation  of  elastic 
modulus  (or  compliance  of  Ith  layer. 

V(I)  =  Poisson's  ratio  of  Ith  layer. 

Note:  Either  modulus  or  compliance  may  be  input.  If 
COVE(I)  =  0  for  all  layers;  then  a  deterministic 
solution  is  obtained. 


IFLG  =  1 


6  NT ( J )  (15)  =  Number  of  times  for  which  laboratory 

creep  data  is  input  for  J-layer. 

7  TYME(I.J)  (I=1,NT)(8E10.5)  =  Vector  of  times  for  laboratory 

creep  data. 

8  VALUE  (I, J),  COV(I.J)  (1=1, NT)  (8E10.5) 

VALUE  (I, J)  =  Creep  compliance  at  time  I  for 
layer  J. 

C0V(I,J)  =  Coefficient  of  variation  of  creep 
compliance  at  time  I  for  layer  J 

Note:  Cards  6,  7,  &  e  are  repeated  for  as  many  layers 
as  there  are.  Also,  if  anyone  layer  has  creep 
compliance  data  to  be  input,  then  they  all  must 
have  complainces  input  (not  moluli);  however,  they 
need  not  be  viscoelastic  (i.e.,  all  values  can 
be  the  same  for  elastic  layer). 

9  V(J)  (J=1,NS)  (8E10.5)  =  Poisson's  ratio. 

10  NTO  (15)  =  Number  of  tines  for  which  output 

is  desired. 

11  =  Vector  of  times  for  which  output 

is  desired. 

IFLG  f   0  or  1: 

6    COF(K),  DN(K),  COVE(K)  (K-l.NS)  (8E10.5) 

COF(K)  =  Coeffocient  in  power  law  Atn- 

DN(K)  =  Exponent  in  power  law  Atn • 

COVE(K)  =  Coefficeint  of  variation  of  creep 
compliance  (constant  with  time  for 
each  layer). 

Note:  If  any  layer  is  elastic  input  compliance  =  COF(K) 
and  ND(K)  =  zero.  Any  COVE(K)  f   0  again  signals 
probabilistic  solution. 
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V(J)  (J=1,NS)  (8E10.5)  =  Poisson's  ratio 
NTO  (15) 


T1(I)  (1=1, NTO)  (8E10.5)  = 


Number  of  times  for  which  output 
is  desired. 

Vector  of  times  for  which  output 
is  desired. 


Note: 

LAST  END  (A4)  :  The  above  sequences  of  cards  may  be  repeated  any  number 
of  times  for  different  problems;  however,  the  very 
Last  card  of  job  must  have  END  in  the  first  three  columns 
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FEDERALLY  COORDINATED   PROGRAM  OF  HIGHWAY 
RESEARCH  AND   DEVELOPMENT   (FCP) 


The  Offices  of  Research  and  Development  of  the 
Federal  Highway  Administration  are  responsible 
for  a  broad  program  of  research  with  resources 
including  its  own  staff,  contract  programs,  and  a 
Federal-Aid  program  which  is  conducted  by  or 
through  the  State  highway  departments  and  which 
also  finances  the  National  Cooperative  Highway 
Research  Program  managed  by  the  Transportation 
Research  Board.  The  Federally  Coordinated  Pro- 
gram of  Highway  Research  and  Development 
(FCP)  is  a  carefully  selected  group  of  projects 
aimed  at  urgent,  national  problems,  which  concen- 
trates these  resources  on  these  problems  to  obtain 
timely  solutions.  Virtually  all  of  the  available 
funds  and  staff  resources  are  a  part  of  the  FCP. 
together  with  as  much  of  the  Federal-aid  research 
funds  of  the  States  and  the.  NCHRP  resources  as 
the  States  agree  to  devote  to  these  projects.* 


FCP  Category  Descriptions 

1.  Improved    Highway   Design   and   Opera- 
tion for  Safety 

Safety  R&D  addresses  problems  connected  with 
the  responsibilities  of  the  Federal  Highway 
Administration  under  the  Highway  Safety  Act 
and  includes  investigation  of  appropriate  design 
standards,  roadside  hardware,  signing,  and 
physical  and  scientific  data  for  the  formulation 
of  improved  safety  regulations. 

2.  Reduction    of    Traffic    Congestion    and 
Improved  Operational  Efficiency 

Traffic  R&D  is  concerned  with  increasing  the 
operational  efficiency  of  existing  highways  by 
advancing  technology,  by  improving  designs  for 
existing  as  well  as  new  facilities,  and  by  keep- 
ing the  demand-capacity  relationship  in  better 
balance  through  traffic  management  techniques 
such  as  bus  and  carpool  preferential  treatment, 
motorist  information,  and  rerouting  of  traffic. 


*  The  complete  7-volume  official  statement  of  the  FCP  is 
available  from  the  National  Technical  Information  Service 
(NTIS),  Springfield,  Virginia  22161  (Order  No.  PB  242057, 
price  $45  postpaid).  Single  copies  of  the  introductory 
volume  are  obtainable  without  charge  from  Program 
Analysis  (HRD-2),  Offices  of  Research  and  Development, 
Federal    Highway   Administration,    Washington,   D.C.    20590. 


3.  Environmental  Considerations  in  High- 
way Design,  Location,  Construction,  and 
Operation 

Environmental  R&D  is  directed  toward  identify- 
ing and  evaluating  highway  elements  which 
affect  the  quality  of  the  human  environment. 
The  ultimate  goals  are  reduction  of  adverse  high- 
way and  traffic  impacts,  and  protection  and 
enhancement  of  the  environment. 

4.  Improved  Materials  Utilization  and  Dura- 
bility 

Materials  R&D  is  concerned  with  expanding  the 
knowledge  of  materials  properties  and  technology 
to  fully  utilize  available  naturally  occurring 
materials,  to  develop  extender  or  substitute  ma- 
terials for  materials  in  short  supply,  and  to 
devise  procedures  far  converting  industrial  and 
other  wastes  into  useful  highway  products. 
These  activities  are  all  directed  toward  the  com- 
mon goals  of  lowering  the  cost  of  highway 
construction  and  extending  the  period  of  main- 
tenance-free operation. 

5.  Improved  Design  to  Reduce  Costs,  Extend 
Life  Expectancy,  and  Insure  Structural 
Safety 

Structural  R&D  is  concerned  with  furthering  the 
latest  technological  advances  in  structural  de- 
signs, fabrication  processes,  and  construction 
techniques,  to  provide  safe,  efficient  highways 
at  reasonable  cost. 

6.  Prototype  Development  and  Implementa- 
tion of  Research 

This  category  is  concerned  with  developing  and 
transferring  research  and  technology  into  prac- 
tice, or,  as  it  has  been  commonly  identified, 
"technology  transfer." 

7.  Improved  Technology  for  Highway  Main- 
tenance 

Maintenance  R&D  objectives  include  the  develop- 
ment and  application  of  new  technology  to  im- 
prove management,  to  augment  the  utilization 
of  resources,  and  to  increase  operational  efficiency 
and  safety  in  the  maintenance  of  highway 
facilities. 
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